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Effect of Telegraph Distortion on the Margins of Operation 
of Start-Stop Receivers 


By W. T. REA 


Recent practical and theoretical investig 
tortion on the margins of operation of start 
to the development of improved ms 
have enabled criteria of distortion tolerar 
monitoring receivers and regenerative 1 
application of more convenient and accurat 
mission. This paper describes the causes of dist 
and internally to the receiver and the effects of such d 
margins. Methods of determining the internal dis 
described and some of the more importan 
ing distortion tolerance criteria are discussed 
URING the past decade the proportion of Bell System telegraph 
service operated on a start-stop teletypewriter basis has shown a 
continuous increase. Whereas in 1930 about 65°% of telegraph long- 
distance circuit mileage was manual Morse, the present proportion of 
teletypewriter and teletypesetter service stands at 92°%. The rapid growth 
of teletypewriter switching facilities has been an important factor in this 
development. 

Naturally, this situation has made increasingly important a thorough 
understanding of the factors which affect the performance of start-stop 
receivers. In the present paper, an effort will be made to show some 
relationships between signal distortion and the operating margins of start 
stop receivers. 

A properly designed start-stop telegraph receiver requires only a small 
portion of the time of each signal element to permit a selection to be made; 
ie. to determine whether the signal element in question is marking or 
spacing. The remainder of the signal element gives an operating margin, 
and serves as a reserve to take care of imperfections in the receiver or 
distortions which the telegraph signals may suffer in their passage over lines 
and through repeaters. The greater the signal distortion is, the smaller 
will be the margin which remains in the receiver to overcome the effect of 
such factors as wear of parts, variation of adjustments, or differences in 
speed between transmitter and receiver. 

A consideration of the effects of telegraph distortion on the margins of 
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operation of start-stop receivers may well begin with a brief review of + H 
nature and causes of the various types of distortion commonly experie: 
by telegraph signals. Telegraph distortion is generally considered to |, fu 
divided into three types or components: bias, characteristic distortion, anc m 
fortuitous distortion.! The magnitude of the distortion is expressed in 
cent of a unit pulse. Vi 


THE COMPONENTS OF TELEGRAPH DISTORTION 


Bias, which is the simplest and most common component of distortio: 
may be positive (marking) or negative (spacing). Positive bias appears 
as a uniform lengthening of all marking pulses and an equal uniform shorten- 
ing of all spacing pulses. Conversely, negative bias appears as a uniform 
lengthening of all spacing pulses and an equal uniform shortening of al! 
marking pulses. 

Bias is caused by an improper relation between the levels at which the 
relay or other receiving device responds and the steady-state marking and 
spacing levels of the signal. For example, Fig. 1(B) shows the signals of 
Fig. 1(A) as they might appear as a symmetrical wave on a line. With 
such a wave zero bias will be received when the currents at which the 
receiving relay operates from spacing to marking and from marking to 
spacing are symmetrically located with respect to the average of the steady- 
state marking and spacing currents. That is, zero bias will be received if the 
relay operates from spacing to marking and from marking to spacing at 
B-B, or if the relay operates from spacing to marking at A-A and from 
marking to spacing at C-C, or if the relay operates from spacing to marking 





at C-C and from marking to spacing at A-A. Negative bias will be received 
if the relay operates in both directions at A-A, and positive bias will be 
received if it operates in both directions at C-C. j 


ss 


In Fig. 1(C) is shown an unsymmetrical wave, in which the transient from 
space to mark is more rapid than that from mark to space. In this case, 
positive bias will result when the relay operates in both directions at B-B or 
at C-C, but no bias will result if the relay operates in both directions at 4-1 

In the remaining diagrams of Fig. 1 it is assumed that the relay operates 
in both directions at a level midway between the steady marking and spacing 
levels. Fig. 1(D) shows a wave in which the transients are of such duration 
that the steady-state value is not attained in the shortest pulse length. It 
will be seen that the operation of the relay is delayed less after a short pulse 
than after a long one, and that this is true whether the pulse be marking or 
spacing. This effect is known as negative characteristic distortion, and it 
tends to shorten short pulses and lengthen long pulses. When a series of 
unbiased dots (called telegraph reversals) is transmitted, a steady-stat 
condition is reached, in which the delays become equal on all transitions 
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Hence, the signals are received as sent. When biased reversals are trans- 


mitted, the longer pulses are further lengthened and the shorter pulses are 
further shortened, causing the bias of the received signals to be of greater 
magnitude than that of the transmitted signals. 

Fig. 1(E) shows a wave in which the current overswings the steady-state 
value, and fails to complete the return to steady state within the duration 
of the shortest pulse. It will be seen that the operation of a relay will be 
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Fig. 1—Signal diagrams illustrating causes of distortion. 


delayed more after a short pulse than after a long one, and that this is true 
whether the pulse in question be marking or spacing. ‘This effect is known 
as positive characteristic distortion, and it tends to shorten long pulses and 
lengthen short ones. When unbiased reversals are transmitted, a steady- 
state condition is reached, in which the delays become equal on all transi- 
tions. Hence, the signals are received as sent. When biased reversals are 
transmitted, the shortening of the long pulses and lengthening of the short 
pulses causes the bias of the received signals to be less than that of the trans- 
mitted signals. 
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Fig. 1(F) shows a wave which performs a damped oscillation b« 
settling to a steady state. This type of wave tends to produce a nega 
characteristic effect on certain transitions and a positive character 
effect on others. 

In general, if, on a given transition, the sum of all previous transients js 
such as to delay the operation of the receiving device, positive characterist 
distortion is said to occur. If, on the other hand, the sum of all previous 
transients is such as to advance the operation of the receiving de\ 
negative characteristic distortion is said to occur. 

Bias and characteristic distortion, considered together, are called “sys 
tematic” distortion, because they occur with some regularity, and obey 
certain constant laws. There is another type of distortion that is 
systematic. This is known as fortuitous distortion. It may be caus 
by the effect of various interfering currents on the receiving device. Tig 
1(G) shows a wave upon which interfering currents have been superpos: 
It will be noted that, for a given magnitude of interfering current, the mor 
sloping the wave is in the region of the operating level of the receiving 
device, the greater will be the resulting fortuitous distortion. 

Fortuitous distortion may also occur, in cases of extremely sloping way: 
shape, due to the ‘‘indecision”’ of the receiving device, or, in other words 
due to small variations of its effective operating level from signal to signal 

Fig. 1(H) shows a wave that is affected by interfering currents and i 
which the mark-to-space and space-to-mark transients have different slopes 
in the region of the operating level of the receiving device. The interfering 
current therefore causes fortuitous distortion of different magnitudes o1 
mark-to-space and space-to-mark transitions. It will be shown later that 
distortion of this type affects a start-stop receiver in a particular manner 
which differs from the effect of distortion of the type illustrated in Fig. 1(G 

These, then are the generally-recognized components of telgraph dis 
tortion. More complicated effects ensue when characteristic distortio 
occurs on waves having dissimilar transients in the mark-to-space and 


space-to-mark directions, but a consideration of such phenomena is outside 


the scope of an elementary explanation of telegraph distortion, and is not 
necessary to an understanding of the effects of distortion on the margins of 
operation of start-stop receivers. 
START-STOP DiIsPLACEMENTS 

The basic principles of operation of start-stop receivers have been described 
in previous articles?*. A brief review of these principles will, therefor: 
suffice here. 

The start-stop signal train consists of a start pulse, which is generall; 
spacing, several selective pulses, each of which may be either markiny or 
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spacing, and a stop pulse which is generally marking. The receiving 
mechanism is started by the transition at the beginning of the start pulse, 
and its speed is such that it arrives at the stop position before the end of the 
stop pulse occurs, and remains stopped until the succeeding start transition 


takes place. Thus any speed difference between the transmitter and 
receiver is prevented from cumulating for more than the duration of one 
signal train. 

Since the receiving device starts anew at each start transition, and the 
instants of selection of the selective pulses are spaced in time relative to the 
instant of starting, as shown in Fig. 2(A), the start transition acts as a basic 
reference point to which all other instants of time during the selective cycle 
may be referred. 

The advances and delays of the transitions of the start-stop signal train 
from their normal times of occurrence, relative to the start transition, are 
known as “‘start-stop displacements.” Fig. 2(B) shows the four types of 
displacement that may occur: MB or ‘marking beginning displacement,” 
which is the advance of a space-to-mark transition (beginning of a marking 
pulse) relative to the start transition; SB or “spacing beginning displace- 
ment,” which is the delay of a space-to-mark transition relative to the start 
transition; SE or “spacing end displacement,”’ which is the advance of a 
mark-to-space transition (end of a marking pulse) relative to the start 
transition; and ME or “marking end displacement,” which is the delay of a 
mark-to-space transition relative to the start transition. 

Effect of Bias on Displacement 

Since bias affects all pulses alike, and since in the usual start-stop receiver 
the start transition is mark-to-space, the succeeding mark-to-space transi- 
tions of the signal train are not shifted relative to the start transition. 
Hence the total effect of the bias appears on the space-to-mark transitions. 
Positive bias causes MB displacement alone, as shown in Fig. 2(C). Nega- 
tive bias causes SB displacement alone,-as illustrated in Fig. 2(D). 

The total range through which the selective periods may be shifted, 
relative to the start transition, without producing an incorrect selection is 
known as the orientation range of the receiver. Its limits are read on a 
scale calibrated from 0 to 100 in per cent of a unit pulse-length. Figure 3 
is a graph of teletypewriter orientation range versus input signal bias, for a 
receiver whose range is from 10 to 90 on unbiased signals. Diagrams of this 
type are called ‘‘bias parallelograms.”’ 


Effect of Characteristic Distortion on Displacement 


Characteristic distortion does not affect all pulses of miscellaneous signals 
alike, because, as explained above, the effect on each transition depends 
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upon the signal combinations that have previously been sent over the cir: 
Hence the start transition and the transitions occurring between sele 


pulses are, in general, delayed by varying amounts. All four types 
displacement shown in Fig. 2(B) occur, depending upon whether the tra 
tion in question is mark-to-space or space-to-mark and whether it has lb 
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Fig. 2—Diagrams illustrating start-stop displacements. 
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Fig. 3—The bias parallelogram. 


delayed more or less than the start transition. For example, if a space-to- 
mark transition is delayed less (on an absolute time basis) than the start 


transition, MB displacement occurs; if more, SB displacement. If a mark 


to-space transition is delayed less than the start transition, SE displacement 


occurs, if more, ME displacement. 
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Maximum Displacements Caused by Characteristic Distortion 


The maximum MB displacement will occur when the start transition is 
layed as much as possible and some space-to-mark selective transition is 


delayed as little as possible. This will take place, in the case of negative 

\aracteristic distortion, when as long a marking signal as is possible precedes 
the start transition and a combination of pulses as predominantly marking 
as possible precedes the space-to-mark transition in question. A marking 
signal sufficiently long to permit a steady state to be attained, followed by 
any signal train having the first selective pulse marking satisfies this condi- 
tion, as shown at “‘X”’ in Fig. 4(B), but it will be noted that the MB dis- 
placement extends into the start pulse, where, in the case of a start-stop 
receiver, no selection is made. Hence it will not affect the margin of 
operation of the receiver, provided it is not so large as to prevent the receiver 
from starting. This particular distortion will, however, affect a start-stop 
distortion measuring set‘ or regenerative repeater which is so designed that 
measurements or selections are made during both the selective pulses and 
the start pulse. As far as a start-stop teletypewriter, in which no selection 
occurs during the start pulse, is concerned, the maximum MB displacement 
occurs on the fourth transition of the letter A following as long a marking sig- 
nal as possible, as shown at “Y”’ in Fig. 4(B). This space-to-mark transition, 
being preceded by a spacing pulse of unit length which, in turn, was preceded 
by signals which are predominantly marking, is delayed for a short time, 
whereas the mark-to-space start transition, which was preceded by a long 
marking signal, is delayed for a longer time. Except in the case of unusual 
wave forms, there will be very little difference between the magnitudes of 
the displacements shown at ‘‘X”’ and ‘‘Y’” unless they are both very large, 
since the wave will usually attain steady state during the steady marking 
interval constituted by the first, second, third and fourth selective signal 
intervals. 

In the usual case of positive characteristic distortion, the maximum MB 
displacement will occur when the start transition is preceded by a combina- 
tion of pulses as predominantly spacing as possible, and some space-to-mark 
transition is preceded by the longest spacing signal possible in the start-stop 
code. These conditions are met by repeated, “BLANK” signal trains, 
shown in Fig. 4(E). 

The maximum SB displacement will occur when the start transition is 
delayed as little as possible, and some space-to-mark selective transition is 
delayed as long as possible. This takes place, in the case of negative char- 
acteristic distortion, when a combination of pulses as predominantly spacing 
as possible precedes the start transition and the longest possible spacing 
signal precedes the space-to-mark transition in question. As noted in the 
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preceding paragraph, repeated ‘““BLANK” signals satisfy this description 
Fig. 4(F) illustrates the effect of negative characteristic distortion on this 
signal combination. It will be noted that the resulting SB displacement 
extends into the stop pulse, where usually no selection takes place, and 
hence it would not affect the margin of operation of the start-stop receiver 
except, of course, in the case of a type of receiver, such as a regenerative 
repeater, in which selection of the stop pulse does occur. For the case of a 
receiver which does not select the stop pulse, the maximum SB displacement 
affecting the margin of operation occurs at the second transition of the 
letter ‘“T”’ following repeated “BLANK” signals. Except in the case of 
very large distortions, this displacement will be of nearly the same mag- 
nitude as that illustrated in Fig. 4(F). 

In the usual case of positive characteristic distortion, the maximum SB 
displacement will occur when the start transition is preceded by a long 
marking signal and some space-to-mark selective transition is preceded by a 
combination of pulses as predominantly marking as possible. As noted 
previously, this description is satisfied by a sufficiently long marking signal 
to permit the attainment of steady state, followed by any signal train having 
the first selective pulse marking. Figure 4(C) illustrates the effect of 
positive characteristic distortion on this type of signal. 

The maximum SE displacement will occur when the start transition is 
delayed as much as possible and some mark-to-space selective transition is 
delayed as little as possible. This will take place, in the case of negative 
characteristic distortion, when a long marking signal precedes the start 
transition and a combination of pulses as predominantly spacing as possible 
= precedes the mark-to-space transition in question. This description is 
answered by a long marking signal followed by a “CARRIAGE RETURN” 
signal train, as shown in Fig. 4(H). The SE displacement occurs at the 
end of the fourth selective pulse. 

In the usual case of positive characteristic distortion, the maximum SE 
displacement will occur when the start transition is preceded by a combina- 
tion of pulses which is as predominantly spacing as possible, and some 
mark-to-space selective transition is preceded by the longest possible marking 
signal. This description is satisfied by repeated ““BLANK”’ signal trains 
followed by the letter ‘“‘A,”’ and, as shown in Fig. 4(L), the SE displacement 
occurs at the end of the fourth selective pulse. 

The maximum ME displacement will occur when the start transition is 
delayed as little as possible and some mark-to-space selective transition is 
delayed as much as possible. This will take place, in the case of negative 
characteristic distortion, when the start transition is preceded by a combina- 
tion of pulses which is as predominantly spacing as possible, and some mark- 
to-space selective transition is preceded by, the longest possible marking 
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signal. As noted in the previous paragraph, the letter preceded | 
repeated “BLANK” signal trains satisfies this description, and the 
displacement is as illustrated in Fig. 4(M). 

In the usual case of positive characteristic distortion, the maximum J) 
displacement will occur when the start transition is preceded by a lo: 
marking signal, and some mark-to-space selective transition is preceded by 
combination of pulses as predominantly spacing as possible. As se 
previously, this description is answered by a long marking signal follow« 
by a “CARRIAGE RETURN” signal train. Fig. 4(J) illustrates the M/ 
displacement. 


Effect of Characteristic Distortion on Orientation Limits 


In the usual start-stop system which employs a stop pulse longer than the 
unit selecting pulse, characteristic distortion affects the upper and lowe: 
limits of orientation differently. This effect is due mainly to the longer stop 
pulse, although the fact that the start transition is always mark-to-space 
contributes to the effect. 

In the case of negative characteristic distortion, the displacements (M3 
and SE) which affect the upper end of the orientation range are those in 
which the start transition suffers a long delay and a selective transition 
suffers a short delay. The delay of the start transition can be quite large, 
since it may be preceded by a long marking pulse. Moreover, the delay 
of the selective transition may be very short, since the pulse which precedes 
the transition can be of unit length, and this, in turn, may be preceded by a 
signal of the opposite type which may be of as much as four units duration. 
Hence these displacements, being the difference between a large and a smal! 
delay, are large. 

On the other hand, the displacements, SB and ME, which affect the lower 
end of the range are those in which the start transition suffers only a fairly 
short delay and a selective transition suffers a long delay. The delay of the 
start transition can not be very short for two reasons: first, the start pulse 
cannot be preceded by a steady spacing pulse; and second, what is of more 
importance, the stop pulse is of greater than unit length. The delay of a 
selective transition can be long, as when the transition is preceded by a pulse 
of four or five units in length. (This delay may not be so long as that 
suffered by a start transition which follows a steady-state marking condi 
tion, but it is not much shorter.) Hence the SB and ME displacements, 
being the difference between a long selective transition delay and only a 
fairly short start transition delay, are smaller than the MB and SE dis- 
placements. 

For this reason negative characteristic distortion affects the upper end 
of the range more than it does the lower. 
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In the case of what we have termed “the usual type of positive char- 
acteristic distortion,” the displacements (SB and ME) which affect the lower 
end of the orientation range are those in which the start transition suffers a 
short delay and a selective transition suffers a long delay. The delay of the 
start transition can be quite short, since it may be preceded by a long mark- 
ng signal. Moreover, the delay of the selective transition may be very 
long, since the pulse which precedes the transition can be of unit length and 
this, in turn, may be preceded by a signal of the opposite type which may be 
four or more units in length. Hence these displacements, being the differ- 
ence between a short and a long delay, are large. 

On the other hand, the displacements (MB and SE) which affect the upper 
end of the range are, in this type of distortion, those in which the start 
transition suffers only a fairly long delay and a selective transition suffers a 
short delay. The delay of the start transition cannot be very short for the 
two reasons mentioned previously. The delay of the selective transition 
can be short, as when the transition is preceded by a pulse four to six units in 
length. Hence the MBand SE displacements, being the difference between 
a short selective transition delay and only a fairly long start transition delay, 
are smaller than the SB and ME displacements. 

For this reason positive characteristic distortion of this type affects the 
lower end of the range more than it does the upper. 

In the case of a wave which oscillates, causing positive characteristic 


distortion on some transitions and negative on others, no such general 
statements as are made above are applicable. In practice, cases have been 
observed in which one end of the orientation range was cut and the other 
was actually extended. 

Due to the fact that characteristic distortion delays the start transition by 
different amounts from character to character, it causes the character length 
to vary during continuous automatic transmission. The maximum varia- 
tion in character length is roughly of the same magnitude as the maximum 
displacement affecting the selective pulses. 


Effect of Fortuitous Distortion on Displacement 


Fortuitous distortion causes the start transition to be delayed more or 
less than normal, and has the same effect on the selective transitions. Since 
it is usually equally probable that the maximum fortuitous effects will 
occur on mark-to-space or space-to-mark transitions and will increase or 
decrease their delay, this type of distortion generally produces the four types 
of displacement in equal magnitude, and this magnitude is equal to the 
maximum increase or decrease in the length of pulse. 

An exception to the above statement occurs when the mark-to-space and 
space-to-mark transients give the wave different slopes at the point where the 
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receiving device operates. Then the magnitude of the fortuitous eff 


different on mark-to-space and space-to-mark transitions. If the eff: 
greater on the space-to-mark transitions, MB and SB displacement 
greater than SE and ME. If the opposite, SE and ME are greater | 
MBand SB. Inall cases, however, the orientation range is reduced eq 
at both ends. 

Fortuitous distortion also lengthens and shortens the character si 
does not affect all transitions alike. 


INTERNAL DISTORTION 


Telegraph signal distortion may occur within the start-stop recei 
and it should be expected that the components of distortion will have 
same effect on the margins of operation as the same components exter: 
to the receiver. Consequently, it should be possible to determine 
magnitudes of the various components of internal distortion by their eff: 
on the margins of operation. 

As mentioned previously, the upper end of the orientation range is 
determined by whichever of the displacements MB and SE is the greater 
and the lower end by whichever of the disp'acements SB and ME is the 
greater. To discover the magnitude of the smaller type of displacement it is 
necessary to reduce the larger displacement by distorting the transmitted 
signals. For example, if a receiver has a large internal marking bias, th 
upper limit of orientation is determined by MB displacement, and hence the 
amount of SE displacement caused by internal distortion is concealed 
However, by transmitting signals affected by SB displacement (in other 
words, signals biased to spacing), the total MB displacement is decreased 
until it is less than the internal SE displacement, whose effect on margi! 
can then be found. Thus the internal distortion may be determined by 
observing the effect of external distortion on the margins of operation. 

It is convenient to regard any start-stop receiver as a theoretically periect 
receiver affected by certain types of internal distortion. The interna 
distortion is usually considered to be composed of bias, “skew” (defined 
later) and fortuitous distortion. (The internal characteristic distortion is 
generally included in “internal fortuitous distortion,” since it is usually very 
small, and a fairly elaborate testing procedure is required to separate its 
effects from those of internal fortuitous distortion.) Internal bias 
internal fortuitous distortion are of the same nature as the external effects 
previously described. Skew is said to occur when there exists the type o! 
distortion, previously mentioned, in which the fortuitous effect on space-t: 
mark transitions differs in magnitude from that on mark-to-space trans 
tions. When the former is greater the skew is said to be positive; when t! 
latter, negative. Hence in positive skew, MB and SB displacements tend | 
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be larger; in negative skew, ME and SE dis 
[he magnitude of the skew is defined as the difference between the mag- 
nitudes of the fortuitous effects on space | 
transitions. 


Figure 3 showed the bias parallelogram of 
margin of 10 to 90. Figure 5 shows the bi 
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receiver and illustrates how the components of internal d 


shape of the bias parallelogram. The skewit 
allelograms shown in Fig 5(D) and (C) led t 
for this effect. 
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placements tend to be larger. 
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In telegraph transmission systems skew may be caused by the effect of 


interference on a wave which has different sloy 
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space-to-mark transitions. It may also result from an equivalent electro 
mechanical effect in a start-stop receiver, as will be described later. 


Measurements of Receiver Distortion Tolerance 


In measurements of the distortion tolerance of start-stop receivers there is 
used a distributor which is arranged to transmit signals having any of the 
four types of displacement MB, SB, SE and ME. Positively biased 
signals are transmitted for MB displacement and negatively biased signals 
for SB displacement. The test signals having SE or ME displacement 
are said to be affected by “end distortion.” These differ from any ex 
perienced on transmission circuits in that only the mark-to-space transitions 
of the selective pulses are shifted relative to the start transition, being de- 
layed for ME displacement and advanced for SE displacement, as shown 
in Fig. 2(E) and (F). “End distortion” simulates the mark-to-space displace 
ments produced by characteristic and fortuitous distortion, and it has been 
found in practice that it yields results which enable a receiver’s tolerance to 
these components of distortion to be predicted with a high degree of accuracy. 

When fixed values of displacement are transmitted, the limits of orienta 
tion are measured by means of the range scale of the receiver. Alternately, a 
distributor may be used in which the magnitude of displacement may be 
continuously varied, and this enables measurements of internal distortion 
to be conducted with the orientation fixed, or, indeed, on receivers having 
no means or a limited means of varying the orientation. 


Orientation Settings for Best Tolerance to Test Distortions 


Obviously, the best orientation setting is that which permits the receiver 
to tolerate the greatest amount of any distortion which is expected. If all 
four types of displacement are considered equally likely, the orientation 
should be set at that point at which the minimum tolerance to any type of 
displacement is as large as possible. For example, consider a receiver 
which, with an orientation setting of 49, has the following tolerances to test 
displacements: 


MB 44 
SB 38 
SE 42 
ME 44 


Let the orientation setting be raised 2 per cent, to 51. Then the tol- 
erances are as follows: 


MB 42 
SB 40 
SE 40 


ME 46 
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The shift of orientation has increased the minimum tolerance (spacing 
bias) from 38 to 40. Any further shift would make the tolerance to spacing 
“end distortion” less than the tolerance to spacing bias. This setting is 


; 


called the “center of fortuitous distortion tolerance,’ 
receiver will tolerate the maximum amount of fortuitous distortion. 

If, on the other hand, bias is considered more probable than distortions 
which produce “end distortion” effects, the orientation might be adjusted to 


since at this point the 


the point at which the tolerances to marking and spacing bias are equal. 
For example, suppose the orientation setting of the receiver under con- 
sideration were raised 1 per cent to 52. The tolerances would then be 


MB 41 
SB 41 
SE 39 
ME 47 


This setting is called the 
receiver will tolerate the maximum amount of bias regardless of the sign of 
the bias. 

There is one more setting that is of interest. It is that at which the 
tolerances to marking and spacing “end distortion” are equal. Suppose the 
orientation of the receiver were lowered 4 per cent to 48. The tolerances 
would then be 


center of bias tolerance,” since at this point the 


MB 45 
SB 37 
SE 43 
ME 43 


This setting is called the ‘‘center of end distortion tolerance,” since at this 
point the receiver will tolerate the maximum amount of “end distortion” 


regardless of its sign. 


Calculation of Components of Internal Distortion 


Figure 6 illustrates how the components of internal distortion are deter- 
mined from measurements of distorted signals. Each diagram shows a 
portion of a teletypewriter character consisting of a start pulse, a marking 
selective pulse and a spacing selective pulse. The solid lines show an undis- 
torted signal. The dashed lines show the displacement of a transition due to 
internal bias. The shaded area defines the fortuitous effect which is skew; 
that is, the transition in question may fall anywhere within the shaded area 
during repeated transmission of the signal. The arrows below the figure 
show the extent of the displacement occurring on each transition due to the 
presence of a given displacement of the transmitted signals. The four types 
of displacement are of equal magnitude D. The arrows above the diagram 
designated Lz and L show the lower limits of orientation with, respectively, 
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spacing bias and markin 


“end distortion” (SB and ME displacemer 


Phe arrow l’, and Ll’, show the upper limits or orientation with, res 
tively, marking bias and spacing “end distortion” (MB and SE 
placements). 

Figure 6(A) shows the case of positive internal bias and positive sk: 


Fig. 6(B), positive bias and negative skew; Fig. 6(C), negative bias and pr 
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Fig. 6—Use of distorted test signals in measuring inte-nal distortion. 


tive skew; and Fig. 6(D), negative bias and negative skew. The followin: 
relationships hold, bearing in mind that MB = SB = SE = ME =D 


Fig. Fias Skew Lp Up Lr UE 

(A) 4 + D+s—b 1—b—s—D D i-D 

(B) “ - D—b 1—b-—D D+(—s) 1-—(-—s)-—D 
(C) _ + D+s+(—b) 14+(-b)-s-—D D 1—D 


(D) - D+ (—b) 1+(—b)—D D+(—s) 1—(—s)—D 
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In any figure Lp — Le = s—b 
ind U,—U,z=-—s—) 


\dding and subtracting, we find that: 


, Us, 4 
Internal bias = — — = 


Us a Le l » Lp 


Skew = 


U+L. 
imag ss the center of an orientation range. Hence it may be 


that the internal bias is equal to the difference between the cente 


, , j l [ 
ance to “end distortion” and bias. It will also be noted that any 


} 


of an orientation range. When the test signal displacements detern 


7 U-L 
range limits are equal, the amount of tolerance equals. ~~ + D (assur 
no curvature in the distortion parallelogram). Hence the skew is e 
the difference between the amounts of tolerance to “end distortion’ 


For example, the receiver cited previously has the following character 


Internal bias = 48 — 52 = —4% 
Skew = 43 — 41 = +2% 


Incidentally, this means that internal bias does not reduce the total! 
tolerance of a receiver, but merely shifts the center of bias toleran 
relation to the center of ‘‘end distortion” tolerance. Hence the effe 
internal bias may be compensated for, as far as the bias tolerance 
receiver is concerned, by setting the orientation at the center of bias 
erance. However, internal bias does reduce the minimum “‘end distort 
tolerance of a receiver whose orientation is adjusted to the center of 
tolerance. 


“Switched” Bias 


When biased signals are produced by the action of a biasing curret 


relay driven by a symmetrical wave, and the sign of bias is suddenly revers 


during the transmission of a teletypewriter character, all the succ« 
transitions of that character are affected, not by bias, but by ‘end 
tortion.”” This is shown in Fig. 7, of which (A) shows the original 
signals, (B) shows the signals affected by bias which changes from positi' 
negative at time 7, and (C) shows the effect on the same signals when t 
is changed from negative to positive. 


Signals such as these, in which the sign of bias is changed at interval! 


are said to be affected by ‘“‘switched bias.”’ Since all four types of disp! 
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ment are present in equal magnitude in switched bias signals, the effect o1 
start-stop receiver resembles that of fortuitous distortion. Thus the cent ra 
of switched bias tolerance is the center of fortuitous distortion tolerance a: 
the amount of switched bias tolerance is the amount of fortuitous distorti: ta 
tolerance. This center is also the center of orientation in a receiver having 1 
curvature or symmetrical curvature of the displacement-vs.-orientatio 
limit characteristic. The switched bias tolerance is, of course, one-half t! 
orientation range in a receiver having no curvature of the characteristic. 
In actual field practice, switched bias signals, applied at a central offic: 
are used as a test of tolerance of the teletypewriter at a subscriber station i 
combination with the subscriber loop. They provide a more accurat: 
measure of transmission capabilities than an orientation range measurement 
with undistorted signals from the central office, since not only is the curva 
ture of the distortion parallelogram taken into account, but the character 
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Fig. 7—Switched bias. 
length changes in much the same manner as in signals affected with char 
acteristic or fortuitous distortion. 

The components of internal distortion of a receiver may be estimated from 
bias and switched bias measurements, but they cannot be accurately specified 
thereby. Figure 8 illustrates the difficulty in separating bias and skew by 
means of measurements of the difference between the amounts and centers . 
of tolerance to steady bias and switched bias. Figure 8(A) shows the bias 
and end displacement parallelograms of a receiver having +24 per cent bias 
and +16 per cent skew. The center of tolerance to switched bias is 4 per 
cent above the center of steady bias tolerance and the steady bias tolerance 
is 4 per cent greater than the switched bias tolerance. Figure 8(B) shows 
the parallelograms of a receiver having +4 per cent bias and —4 per cent 
skew. Again, the center of tolerance to switched bias is 4 per cent above the 
center of steady bias tolerance and the steady bias tolerance is 4 per cent 
greater than the switched bias tolerance. 

Of course, the components of internal distortion can be measured by 












rom 
fied 
by 
ters 
bias 
bias 
per 


NCE 


OWS 








START-STOP RECEIVERS 


‘rving both ends of the orientation range with positive and negative bias 


rather than observing the upper end with positive bias and the lower end with 
ative bias. This type of measurement is merely equivalent to using a 
fairly large percentage of bias and zero per cent of end distortion The dis 
vantage of this measurement is that no account is taken of the curvature 
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Fig. 8—Switched bias measurements 


of the end displacement parallelogram, and hence the indicated values of 
tolerance may not be an accurate measure of the receiver’s ability to receive 
distorted signals. 
Internal Fortuitous Distortion 
It is usually considered, in measurements of miscellaneous signals, that the 
difference between the maximum distortion tolerance and 50 per cent (the 
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effects, even though part of it may be due to the effects of internal 

teristic distortion. Hence the internal fortuitous distortion is 
leftined as the difference between 50 and the tolerance to bias or e1 
tortion, whichever of the latter may be the larger. 


For example, in the sample receiver considered on page 220, the 1 


lon 1s: 


50 — 43 = 7 per cent 
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rig. 9—Effect of negative characteristic distortion on bias parallelogran 


Internal Characteristic Distortion 


In practice it is found that the relation between displacement and redu 
tion of margin is sometimes not strictly linear. Especially at large values 
lisplacement, the reduction in margin is often greater than the displacem 
causing it. This effect is due to internal characteristic distortion, w! 

ses an increase in the distortion of shortened pulses. Internal cl 
cteristic distortion, like any other form of characteristic distortion, is cat 
by the failure of some circuit or mechanical element to attain steady stat 


before the occurrence of a succeeding transition. Figure 9 shows an exam] 


















f the bias parallelogram of a recei 
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For the same reason (the curva 


internal characteristic distortion) mea 
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receiver which is, itself, to be used to1 


displacements of approximately the 


which the receiver is to measure 


measure small distortions, we are 
portion of the parallelograms. Her 
bias and skew using small amounts of 


rhe internal fortuitous distortion may 


not affect the shape, but only the 


acteristic. 
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It should be realized that the remoy 
a service receiver. In the case of 


introduction of distortion of a given 


the opposite sign, an 
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since skew is t 
be removed either by reducing the 


he difference between tw 
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The former proc edure will increase the 
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In practice bias tolerance is genera! 
“end distortion” tolerance. The rez 


circuits suffer from some bias (of 1 


1 thus improve the 


the 
r S) 
} lactead : 
be negiected, since it 
austortion-vs irgin 
€ sed tor rece ys 
+] Rica ica! ieee 
le CnNaracteris 5 
an optimum setting o 
ver destined for serv 


} 












splacements 
I 


of skew does not necessarily improve 
is or characteristic distortion the 
nw remove internal distortion of 

performace of the receiver But 
) fortuitous distortion effects, it may 
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tolerance. This is why the orientation of a service receiver is generall 
adjusted to the center of bias tolerance, and small amounts of internal | 

or negative skew are not considered objectionable, since they do not affe 
the tolerance to bias at the center of bias tolerance. By the same token, t! 
presence of positive skew, which indicates a lowered bias tolerance, usua 
calls for a readjustment of the receiver to reduce the fortuitous effect on t 


space-to-mark transitions. As explained above, removing the skew | 


introducing a fortuitous effect on the mark-to-space transitions will not, ot 


course, improve the bias tolerance. 
It is the present practice in the field to specify a minimum bias tolerat 


about 5 per cent greater than the minimum permissible “end distortio: 


} 


tolerance, the orientation being adjusted to the center of bias tolerance fo: 


both measurements. 


Some CAUSES OF INTERNAL DISTORTION 
Up to this point internal distortion has been considered without rega1 
to its probable causes. The more obvious causes will be found to bé 
analogous to those which produce equivalent distortions in telegraph tra1 
mission circuits. 
Bias will result when an element (whether electrical, mechanical, or elt 


tronic) of a receiver possesses dissymmetry toward marking or spacing 


} 
| 


For example, a mechanical element may travel more slowly from spacing to 


marking than from marking to spacing and thus cause spacing bias, or 


range of travel may be divided unequally into marking and spacing port 
thus producing an equivalent effect. 


Characteristic distortion will result when an element (whether electrical! 


mechanical) of a receiver fails to attain a steady state before being acted 


upon by a succeeding transition, or otherwise depends, in its action, upon the 
previous history of the signal train. An example of characteristic distortio1 
is found in the 20-milliampere holding magnet selector when it is equipp 


with a resistive shunt. In this type of selector the armature is actuated by 


cam, which presents it to the pole-face at about the middle of each pulse 
and then disengages it. The armature is then free to release or remail 


operated, according as the received pulse is spacing or marking. The shunt 





that is normally used presents so low an impedance to the magnet windins 
that the motional impedance effect which is produced by the sudden mechani 
cal presentation of the armature to the pole-faces causes a sizeable reductio1 
in the magnet current. In the case of a short marking pulse, the current 
fails to attain steady state before the next mark-to-space transition occurs 
The magnet therefore releases sooner than it does at the end of a lon: 


marking pulse, during which the current has had time to attain steady state 


It will be seen that this is really a characteristic distortion effect, since it is 


due to a failure to reach steady state and depends upon the previous history 
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.e signal train. However, when miscellaneous signals are being received 
effect appears similar to a fortuitous distortion occurring t 
ce selective transitions, and hence it is usua y thought ol as negative 

Fortuitous distortion result when an element is irregular ts actio 

if such action is more irregular on one type of transition tha ¢ 
er, the result will appear as skew. For example, irregular action of the 
eiving clutch affects the selector alike in regard to all selective ~ 
ns, and appears as internal fortuitous distortio1 (nother source of 

ternal fortuitous distortion is the period of indecision that occurs J 

passage of a selective element past a locking member, at which time the 

ce between marking and spacing is largely fortuitous 

\ common cause of skew in teletypewriters may occur in the follo ¢ 

If the armature stops are so adjusted tl for example, the armature 
vel is greater on the narking side than on the spacing side i the arma ¢ 
positive internal bias results. If, now, this bias is compensated for by 
idjusting the armature air-gap and retractive spring tension as to se the 

eiving magnet to operate in a negatively biased manner (rather tl! t 
rrecting the improper armature travel), the armature will be forced to op 

te in a region of the operating wave that is more sloping than the re 

which it releases. Hence, it will operate more irregularly than it releases 
ind thus will be affected by positive skew. 
SELECTOR ACTION 

Over and above the sources of internal distortion which are analogous 
ffect to sources of distortion encountered in telegraph transmiss r 
its, there is another whose action in cau distort S sO 

vious as those just described. This s stort ¢ 

termed ‘‘selector action,” and it depend relation between the 
erating time of a selector element and yf time allowed for said 
ement to act. For the purpose of explaining the effect of time re ns 
vithin the selector on internal distortion, selector mechanisms iv be 
lassified as of three basic types: M, S, and P. 

In a mechanism of type M each selector is initially in the spacing 
nd either remains spacing or operates to marking when subjected to the 
ction of the corresponding received signal element. When it attains the 

marking condition it becomes locked for the duration of the « le 


Early types of start-stop printers having an individual sele 
yr each pulse of the code and employ.ng a separate recciving « 


ire i.lustrative of type M. 


In a device of type S each selector is initially in the marking 


ind either remains marking or operates to spacing when subje 


} 


action of the corresponding received signal element. When it 
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spacing condition it becomes locked and cannot agalr operate to mar 
during that character lhe Siemens-Halske five-selector teleprinter 
r cample oft this type 
} urs } } l " 
In a mechanism of type P, the selector may be in either the mark 
spacing condition 1 " ord g to the tvpe of the previous sigt 
| g 
ment to which it has responded. When subjected to the action of a rece 
puise the selector n neither directio and it remains responsive t 
I 
action of the signal during the entire selecting interva The No. 14 
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Figure 10 (A) illustrates the action of a type VW selector. <A portion 0 


teletypewriter character is shown, consisting of the spacing start puls 


) ) ] . } > ) 
marking first selective puise and a spacing second selective puise | 

’ ) 1 3 9° } . ¥ > 
undistorted signal is shown in solid lines The maximum amounts of ma 


ing and spacing bias that the receiver will tolerate are shown by dashed | 


and are designated 4B and SB. The limiting amounts of marking a 


q 


spacing end displacement are shown by dotted lines and are designated 
and SE. Above the signal train is shown a schematic representation of 1 
action of the selective system. The periods of time 7 are those during w! 


the selector is subject to the action of the received signal, and ¢ is the t 
] 


hat the selector must be subjected to the operative force in order that 
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type of mechanism, the armature general] 


drives a subsidiary se!ect 
member, and the time 7 extends from the instant at which the armatut 
isengaged by its operating cam until the instant when the subsidiary si 
tor becomes locked. As this period is often long in relation to the ma 
releasing time ¢; and the subs‘diary selector operating time fe, holding n 
selectors are often subject to negative internal bias. In those mecha 


+ 


hich the subsidiary selector is flexibly coupled to the magnet armat 


the former’s operation is of type rr It, therefore, may be subject to 
characteristic distortion effect noted in the description of type P operat 
except that the effect, when it occurs in this type of mechanism, affect 


} f ; rand dicen] } } | } . 
the instant of decision tor end displacement and hence resembles ne¢ 
skew rather than internal fortuitous distortion 


An interesting, but somewhat unusual, ettect occurs In any receiver 


1 +] 4 ] ‘ " ] . 
itever type, In wh the lengths of selective period or selector ope: 

time, or both, differ for the various selective pulses, or in which the spa: 

I the selective perio is Is In proper. In a case of this sort, the receiver 


hits an internal bias equal to the difference between the average instant 


decision tor end dispiacem¢ 


decision for bias and the average instant of 
an internal fortuitous distortion equal to the variation of the insta: 
decision having the smaller variation, and a skew equal to the differ 
between the variations of the instant of decision for bias and the insta: 
or end displacement 


CONCLUSIONS 


\ working knowledge of the effect of telegraph distortion on the n 





. ! 1 
art-stop receivers is essential in dealing with a plat 


he use of teletypewriters, regenerative repeaters and start stop 


tortion measuring sets is as widespread as it is in the Bell System. Whe: 


major portion of the communication system operates on a start-stop ba 
desirable that transmission measurements be made on the same basis 
The knowledge of this subject that has been gained in recent years | 


e many improvements in technique both in the field and in 


made possib 
laboratory, and these have led to corresponding improvements in the me 
in telegraph service. The analysis of new start-stop dev 
may now be carried out efficiently and accurately, and this often pern 
the formulation of suggestions leading to improved operation of the devi 
The general level of service excellence has been raised by the setting 

of criteria for the distortion tolerances of station teletypewriters, regenerat 
I 


vided for switching The sources of distorted test signals that are no 


available are useful not only in measuring the tolerances of service receive! 


repeaters and other start-stop devices used in service, including those pi 
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CHAPTER XIII* 


The Mounting and Fabrication of Plated Quartz Crystal Unit 


By R. M. C. GREENIDGE 


, ; 
HIS paper 1s one of a series on piezoelectri quartz plates and de 
primarily with the methods employed in mounting al plat 








operating up to approximately one megacycle for practical 
communication equipment Phe theoretical aspects of mounting crys 
have been covered in ¢ apter VII The LISCUSSIOI s confined to pl 
ing definite noda Ines or points, st as +5 LNG 18 740 xX ( 
rr, DT, MT and NT cuts The mounting of vh-trequency cry 
such as AT and 3] ( Ss Vn) Vibrate i thickness shear moce 
s not included It sho SO noted that t it) je matter 1s tre 
descr plively and that tempt made t oO tne re intr 
details of design or to e€ performance iracterist These matters 
be dealt with fully er paper The designs d methods outlined 
up to date for eat type ol it, the res S ol \ ears of developm« 
on the part ot Bel! system engineers to evolve practi il designs tor comn 
clal manutacture and use. Expanding on the contributions of the ¢ 
ivest itors mentioned b \\ P. Mase ( ipter these ¢ wines 
had, in the ten vears pr to 1939, worked out pra designs and deve 
oped suitable tools and rocesses tor W ( ommercla tion in telep! 
applications. In the last five years, under tl mpetus of war, furt 
improvements have been made in the design and manufacture of ervst 


units, particularly those tor use by the Armed Forces 


The term “Crystal U1 t’’, orig nally adopted by the Bell Svstem to des 


I i t 4 
] ] ] + ‘ syet , > 
nate the complete assembly of a crvsta piate ; mounting and case, 
‘ } Ar + wiarsl a] wuite a oe iy — { art rer no riety 
now een Standardized quite generally in thie rt, replacing a variety « 
} } | , rmerl) ] he : =F 
names by which these devices were formerly calle 1 The basic design fe 


tures of a crystal u 





1. Electrodes, on or near to the crystal surfaces for impressing voltage 


“Quartz Crystal Applications”, W. P. Mason, B.S.T.J., Vol. XXII, Page 191, J 
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2. Supports for holding the crystal plate in its mount 
> ] 

3. A sealed outer 


r case having te necessary 








term 
incorporating the unit electrically and mechanically int 
[wo distinct types of cry stal units have been e f 
e use of press pins or anvils for supporting a 
ite and the other i { on ol e crvst | 
fine wires These « I respect Vel s e | I 
1 Wire Supported discussed latet Le ( 
However, there are fabrication con t t 
h can best be dis« issed at this ] nt 








si { 
i eI ; 
ae aw i and - 
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Fig 13.1 Pre ' 
[Irrespective ot the type ot mounting. rea ath the 
e€ ina crystal unit d 


es | ve Di ont 

ice of such tactors as the precision ol nguilatl ent 
mensions on the fundamental characteristics I the ( i 
ist also be virtually free of impurities or imperfections 

n of a quartz plate i 


I 
plate it must be lay ly 
iterials until the final di ns are reached Depending 
crystal unit, No. 400, No. 600 carl wey 


arborundun 





A. W. Ziegler, Patent 2,275,122 2 
“The Use of X-Rays for Detern t { 
land E, J nstrong, B.S.T.J., t. 1943 
“Raw Quartz, Its Imperfecti 
XXII, Oct. 1943. 
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employed for the final stage grinding. Following this, the plate is t 
oughly cleaned by acid treatments or by the use of solvents and deterg: 


followed by copious washing. It is then etched in commercial hydroflu 


acid to remove all loose particles of quartz that might have remained o1 
surfaces or in the crevices after cleaning. Etching also smooths off 
roughness of the ground surfaces. The effect of this treatment redu 
energy dissipation in the plate itself and increases by many times the « 
ciency of the crystal units. Etching also improves the stability of perfor 
ance of the crystal unit. Standard designs of crystal units require etc! 
of the plates, uniformly on all surfaces, for a period of thirty to forty n 
utes. For units of highest precision and efficiency longer etching peri 
are employed. 

rhe electrodes employed with types of crystal units being described « 
sist of metallic coatings, generally aluminum, silver, or gold deposited o 
the major surfaces of the crystal plate. These coatings are applied by 1 
evaporation process which results in an extremely thin and uniform coat! 
of metal having excellent adherence to the quartz. 

With reference to the mechanical supporting members for the plate 
has been brought out that such supports should be confined as closely 
possible to the nodal points or nodal lines where the motion for all practi 
purposes is zero. It is common practice for the supporting members to |! 
made of metal so that they will serve also as a means of making electric 


connections to the electrode coatings on the surfaces of the plates. 
13.2 PREssuRE Type CRrysTtaL UNITS 
This type of crystal unit was initially developed for use in telephor 


filters. Up to about five years ago it was employed in virtually all cor 
mercial designs of filters. Depending upon the mode of vibration and s 


of the crystal plate the design of the mounting varies. However the pri 





ciples employed for clamping are essentially the same in all cases. Wher 
small longitudinal or face shear plates are involved one pair of pressu: 


pins is used unless two are required for electrical reasons as explained lat: 
For medium size plates of the same type or for face flexure plates, two pai: 
of pins are employed. In the case of large low-frequency longitudi: 
plates double anvils are used instead of pins in order to obtain firmer clam] 
ing of the plate to prevent translation or rotational movement which mig 
cause wear in the electrode surface at the pressure point with resulta 
variations in frequency and resistance. The blocks are usually compos« 
of molded steatite and the springs for exerting the necessary pressure are 
phosphor-bronze. 

Figure 13.1 (B) shows a pressure mounting for holding four crystals whi 


have single coatings on each of their major surfaces. The main requiré 
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must be met for such a mounting are small areas, a 
ent, and adequate pressure to hold the plate in place Some 
ers make slight indentations in the quartz at the point of conta 


ve the mechanical stability of the plate. For crystal plates of 


me-half inch square or smaller, points having an area of 








neter are employed and the pressures used for holding 
te range Irom one to two poun Is For larger plates correspond 
reas of points and increased pressures are employed. The ac 
ent required for ten mil points is of the order of two or three 1 
tained in the mountings shown in Figure 13.1 by using 
leeves for holding the points which are brought into alignment 
fa straight rod and then cemented in place. One of the poi 
hile the other point slides in its sleeve and the pressure requirt 
ed by the spring which presses on the outer end of the sliding px 
I ilanced filter structures it is desirable to use crystal plates wit 
r on each side divided into two equal areas. This reduces to one 
iber of plates that would otherwise be required. In mounting p 
led coating, it is necessary to provide a mounting which n 
t on each side of the plate. Figure 13.1 (A) shows a pre 
ting which accomplishes this. This is the mounting whi 
been used for several years in holding the plates used for the 75-type crn 
nel filters® for the standard terminal common to all broad-band 
systems. The crystal is mounted in the holder in such a way 1 
e two pairs of points clamp the crystal along the nodal line. The rectar 
ensions of the points used for this type mounting for crystals « 
the frequency range from 60 ke up to 120 ke are about 35 mils 


he nodal line and from 10 to 15 mils wide. A 
requirement for such a mounting is that the flat area of 
each side of the plate fallin the same plane. This is accomplis 
i precise milling operation after the points are assembled in the mount 
[he pressure applied to the pair of points is furnished by the flat spr 
l is equalized by the action of the roller centrally located ut 
e springs. The pressure employed is of the order of four to five pout 


pair of points 





[he most commonly used coating for crystal plates held in pressure-tyj 
ngs is aluminum.® Aluminum has been found to be most satisfact: 
is type of unit because its hard surface is more resistant to wear at 

} 1 


of clamping than other metals such as silver, or gold 
i 
Except for a few designs, which are mounted in sealed metal or 
Crystal ( 1 el Filters for Carrier Cable Systems,” C. E. Lane, B.S.7.J., Vol. X\ 


aluminum-coated crystals are similar to those descril 


crystals in paragraphs 13.42 and 13.43 











Lammers, press re-tvpe cl ta | { ‘ ret ri ‘ 
ners Llowever, the entire filter thre reer 


i] 


ealed off after filling with drv ait 


In pressure-type units of this type, variations in freque f the 
f .01°% may be expected if crvstals are transferred from one 
» wal } +} 
nother or relocated in the same mounting Consequently, if high-tre 
recision Is desired, t is necessary to make the tn trequel! 
ith the crystal located in its final positiot Due to the inherent dit ti 
f adjustment, coupled with the close manufacturing tolerance 


nd precision adjustments necessary fot 
lesigns have been virtually discarded in favor of wire supported desis 
However, more recent developments by J. F. Barry on pressure-mount¢ 


type units have brought forth some new ideas which might prov 


vider future application 


13.3. Wire SUPPORTED CRYSTAL UNITS 


This type of mounting is being used extensively on crystal ap} 


ind has superseded the earlier pressure-type units. Various designs 


type of crystal are shown on Fig. 13.2 rhe wire-mounted crystal possesses 
the definite advantage in that after the supporting wires are att ‘ 
he plate, they remain fixed in position throughout the subseque 


facturing process thus facilitating adjustment of the frequi 
frequency-temperature characteristic. Moreover, the supporting wir 

be formed so as to provide a spring mounting for the crystal plat 
protects it from any shocks or vibration it may encounter in shipment 
ise. In the wire-supported design the suspension wire is also employe 

t means of connecting the electrical circuit to the electrode plating 
crystal. By virtue of the solder bond between the wire and the electrode 
this type of unit is free from the possibility of instability in frequer 
formance due to slight changes in position and variations 
ance prevalent in pressure type designs, and for this reason the 


ind frequency precision of wire supported crystals is superior 


From a manufacturing standpoint the wire-s 
vreater number of processing operations than the pressure 


the adjustment operations are considerably easier. Moreover, it is poss 


to realize greater pre ision of freque ney adjust ment by a fact ( 
two or even three The crystals are also more unifort ¢ elle ( 
resistance. Since the mount or cage in which the crystal 1S suspe¢ ( 
omparatively inexpensive, there should in general be little difference 
‘anufacturing costs of the two types ( onsequently, in the wire-su] rte 


crystal units a very appreciable improvement should be gained in per 


+t 
Lit 


ance without increasing the cost of the u 
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13.4 BRICATION OF Wikke Supported Unit 


13.41) Silver Spottin 


Starting with the crystal plate, the first step in manufacture is te 


SHVer Spots to the surfaces of the plate. These spots serve as footings 
vhich the supporting wires are ultimately soldered or sweated They 

placed on the nodal points or along the nodal lines of the plate in order 
detract as little as possible from the intrinsic characteristics of the p 
itself The areas of the silver spots cover the range from about 40 to 


mils in diameter depending upon the amount of solder to be used in atta 


+} { 


ing the wire to the plate Before spotting the plates it is essential 





they be free from any contamination such as grease or organic mater | 


that might affect the fusion of the silver spots into the surface of the quart 


One of the best methods to ensure cleanliness is to boil the plates in aq 


regia, followed by copious rinsing in water. Detergents such as sodi 
meta silicate Are also employed followed by a rinsing The plates n 
finally be boiled in distilled water and carefully dried Throughout 


1¢ plates should be handled with clean tweezers 


subsequent processes tl 


gloved fingers and kept away from any source of contamination. —Prefir 





of the plates at 950° prior to spotting has also been used as a positive va 


of ensuring freedom from any contamination that would affect the fusi: 


of the spots, but this process is not necessary if the first mentioned proc 
is properly controlled 


In spotting, small quantities of a prepared silver paste are placed on 


areas of the plate to which the wires will ultimately be attached. The past 
consists of a compound ol finely divided silver and low melting point glass 
lead borate) thoroughly mixed with a suitable vehicle to facilitate applic 
tion. For spotting purposes it has been found that a paste having a specif 
gravity of between 2.3 and 2.6 gives best results. In use, the materia 
must be constantly agitated or stirred in order to prevent the solid ingre 
ents from settling out This is important, for, unless the concentratior 


silver is maintained around 90 to 95 per cent of the solid matter, it v / 





| t 


not be possible to obtain good wetting of the solder in making the wi 


attachment. 


} 


The placement of the semi-liquid material on the plate 


s accomplishe 
by means of a small stylus, the crystal plate being held in a clamp or \ 

and the stylus guided so as to place the material at the exact location on thi 
plate as desired. A typical tool for doing this work is shown in Figure 13.3 
The point of the stylus should have a slightly rounded end. With t! 


rounded point the tendency of the paste to spread out is minimized at 


consequently the diameter of the spot is substantially the same as that 























) 


tvlus. The round 


led stylus also results in a more uniform dist 
e material. The material is applied to the end of the stylus by spr 


l } 


small amount of the paste on a glass plate from which it is transferre 


e stylus and then deposited on the crystal The material on the glass 
should be wiped off and replaced quite often due to settling and 

{ the mixture, in order to insure uniformly good spots. Generally 
‘ing, anywhere from two to six spots at the most should be possible 
one loading of the transfer plate depending upon the speed of 


itor. 


Regarding the character of the crystal surface, aside from cleanliness, and 


; effect on the ultimate strength of adhesion of the silver spot, experience 





far with all the various cuts of pla 


ctor. Spots have been found to adhere to polished surfaces as we is 


t the same degree ol! strengt 


rround and etched surfaces with abot 


13.412 Firing of Silver Spots 

Following the application of the paste to the plates it is desirable to pre 
‘o) t I } 

Iry the spots in order to remove the low volatile constituents of the vehicle 


rior to firing at high temperature. This may be done in a ventilated over 


over a hot plate at approximately 300°F for about 15 minutes. The 


plates are then placed in a furnace and heated up to between 975 to 10101 


id held at that temperature for a sufficient length of time to obtain good 


thy 


fusion of the spot to the plate. Ordinarily this reaction takes only a fe 


ninutes after the plate has reached the proper temperature. After firing 


he plates are allowed to cool in the furnace to the 
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emoved without danger of the crystal crac king due to cold shock 
| the temperature of the furnace so that the temper 


of the crystal plates does not reich 1063°F, otherwise the crystals 


become electrically twinned and consequently useless. In order to 
shattering of quartz plates due to thermal shocks while heating u 
ooling, fairly long cycles have heretofore been specified. However, 


recent experience has shown that much shorter cycles can be em] 
especially where small crystals are involved. Moreover, there are 

tions that the faster heating, particularly during the last two or thre 
dred degrees temperature rise, results in better spots. During the 











ig. 13.4—Continuous-belt furnace for firing silver spots 


operation the crystal plates may be placed on nichrome wire mesh trays 
Pyrex dishes provided ample provision is made for air to circulate arou: 
the spots. This pre 


I 


as the baking reaction must take place in the presence of oxygen so that 
lead borate will not be reduced to lead, leaving only a partially bond 


aution is essential with the types of pastes employ: 


mixture of silver and lead on the plate. This type of spot when encounter¢ 





lv has a dull appearance after burnishing as compared with the brik 


uned with a good 


Surface obt 


silver spot, and is quite difficult to wet wit 
solder. For firing silver spots, ventilated continuous-belt-type furna 
with open ports are being used with very satisfactory results. Figure 13 
shows a furnace of this type developed by C. J. Christensen for this purpos 

After the firing operation, the spots should be examined to ensure tl 
they are satisfactory 


Besides a visual inspection, it is desirable that 




















| percentage of the plates be used as a nt 


es are attached and pull tested. Satisf 


few seconds a force of at least 


neth of commercial attachments using 6-1 
een three and four pounds. Unsatisfactory 
s stage by stripping off the spots by means 
roxide, and reprocessing in the manne r 
lo {? .* yr P 


\t this pe int of the process the surfaces ot 
er electrodes by the evaporation process pre\ 
igrams per square inch of silver is the we 

oved, which amounts to a thickness of .024 1 
| other special types of crystal units, these co 

e major surfaces of the plates. Howeve 


e silver Is dep sited to some ¢ 


d it is necessarv to remove it Phis pr 


me by lapping the edges of the crysta 


nixture of pumice or a finely divided abrasive 


d water or kerosene in the form of a paste IX 


plate lightly over very fine abrasive cloth is also sat 


s preferable, however, since while it readily remov 


softer than quartz and consequently does not remo 


} } 


late. The use of harder abrasives has a te el 


the quartz unless the operation is performed very 


eaning is completed the plates are washed, dric 


for insulation resistance at 500 volts d.c. to make ce 
material remains between the silver coatings o1 
13.43 Divris f Coating 


lor circuit reasons all but a few types ol crys 


pair of electrodes on each side of the plate Division 


ongitudinal axis is essential on flexure mode cryst 


plate vibrate in flexure. Typical divisions of coat 


rystals shown in Figure 13.2. In the case of the 


ing line is carried around the wire attachments 


of the divided surfaces is connected to one of the 


dividing coating involves the use of a low voltage 


} ] } 


mpressed between the coating to be divided and a 


point of the stylus is brought in contact with the si 


ilong the desired line of division, the silver is burt 


7W._L. Bond, Pat. # 2,248,057 





Oth 
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gap in the plating between 8 and 18 mils wide depending upon the px 
the stylus. Following the burning operation, the plate is immerse 


photographic hypo solution to remove all traces of the burned residue 














Fig. 13.5—Tool for dividing plating (electric stylus). 


which it is carefully washed and rinsed in water and dried. In this proce 
it is important that the plates are not kept in the hypo solution for longé 
than two or three minutes, otherwise discoloration of the silver coating 
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The gap in then tested for presence of metallic particles by grad 


impressing voltages up to 1000 volts a-c across it. If no flashover 
irs the division is satisfactory. If flashover occurs the voltage is 
ntained until the slivers are burned out. The hypo and burning treat 
ts are repeated until a good division is obtained. Figure 13.5 shows an 

electric dividing tool developed for this purpose. In using this method 
been found that the arc at the point of the stylus may cause twinning 





the quartz to a minor extent along the dividing line. This 
ially insignificant although it may be objectionable especially where pré 


ise values of crystal inductance or frequency-temperature performance ar 


f 


required. Where more complicated divisions are necessary, as in the case 
of face flexure and harmonic plates, the electric stylus method is emp! 
iithough methods and tools for performing this operation by other means 


to avoid twinning are being developed. 
13.44 Attachment of Wire Supporting Leads 


Phosphor-bronze wire is employed in wire supported crystal units pr 
marily because of its high tensile strength, and excellent fatigue resistance 
characteristics. Five- and six-mil diameter wires are the most widely used 
sizes, depending upon the mass of the crystal plate, the desired electrical 
performance, and the severity of treatment it is likely to encounter in use 
lo facilitate soldering the wires to the spot on the crystal plate and to the 
crystal support system the phosphor-bronze wire is given a heavy electro 
tinned ‘finish. 59.5-34.5 per cent tin-lead eutectic solder saturated wit! 
approximately 6 per cent silver at 570°F is employed for attaching thes« 
fine wires to the silver spots of the crystals. This solder solidifies at approxi 
mately 360°F with practically no mushy stage. The reason for saturating 
the solder with silver is to discourage migration of the silver in the spot to 
the solder during the soldering operation. Even with this solder it is advis 
able to limit the time for heating of the joint to a minimum. 

One method of attaching the wires to the crystal plate is by means of a 
special machine developed for the purpose. Such a machine is illustrated 
in the photograph on Fig. 13.6. The wire is fed from a spool through the 
head in the movable arm. The head contains a wire guide having a hole 
only slightly larger than the diameter of the wire and a small vise for firmly 
clamping the wire. The crystal plate is clamped in the vise on the hot plate 
which is thermostatically controlled at approximately 240°F. The position 
of the arm carrying the wire is lined up with respect to the crystal plate 


by means of guides so that the wire will be placed exactly on the nodal 





or line of the crystal. In making the attachment, with everything lined yy 
the wire is fed through the guide until it touches the spot on the plate and 


the vise closed. Since the curvature of the wire can never 
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eliminated the distance between the tip of the guide and the crysta 





is kept as small as possible. A small disc of solder is then punched 
press at the left, the little disc remaining in a round slot whose pos 


also lined up with respect to the arm carrying the wire. The moval 





Fig. 13.6—Wire-soldering machine for straight or hooked wires. 


is then rotated until the wire is directly over the solder disc at which px 
it falls into a guide and comes down and spears the disc. The arm is 1 
lifted, picking up the solder at the end of the wire. The solder and 


end of the wire is then wetted with rosin-alcohol flux and the arm rotat: 











it falls into its origi 

fused to the wire and plate by me 
ring iron as shown in the illustrati 
ed on the joint to melt the solder. 


yn is formed around the wire atta 


good wetting of the solder, the spot sl 


bbing the spot on a hard polished n 

t pointed tool of agate, are the best 
Phe hot air blast has now replaced tl 
It nsists of a tube through which air a 


issed over a hot filament and thre 


[he head of the filamen 


mal] 
sillalt 


I 
brought up slowly enough to preheat t! 


flame have been considered, but so far 


ave been made. The advantage of the 





nelting the solder such as a hot radiant wi 


Si 


is adjusted so that 
ist hot enough to melt the solder and comp! 
15seconds time. In using this method with 
to insure that the temperature of the crystal 
plate and that the blast is brought up to the pl 
eat shock of the localized blast may cause 


plates the use of a hot plate may be d 


mentioned is that it can be better controll 


nent of the operator. If an iron is used 
older with the possibi Ly ol displac ing t] 


s already mentioned, the iron must be « 


tip to prevent removal of solder from the 


Considerable maintenance is required t« 


nerating condition. 


The type of attachment described al 


embedded in the solder cone is straicht 


supported crystals. However, it was four 


vibration of the crystal plate caused br 
ilder cone and the wire with resultant fai 


large plates. Recause of this the use 


y for small size plates. In order toe 
I , } 
, ok has been placed at the end of the v 


i to obtain a better anchorage. The hook 


. afl a 
pecial tool affixed to the soldering ma 


r straight wires are otherwise used for 
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of spearing the little solder discs as with the straight wire, the sol 
punched in the shape of a horseshoe and squeezed in place on the ho 
positioned by tool with the hooked wire in place on the spot. Hooked 
attachments will withstand pulls of the order of three to four pounds | 
pulling off. Under severe vibration hooked-wire attachments have the 


tendencies as straight wires towards breaking away of the wires fron 
top of the solder cone forming a sma!l crater in the latter. However 
crater does not progress deeply enough into the solder cone to impair | 


strength or cause failure un ler ordinary conditions. 








Examples of headed phosphor-bronze wire. 


The most recent development for wire supports involves the use of heac: 
phosphor-bronze wires as worked out by A. W. Ziegler. In this procedur: 
individual wire lengths are cut and one end upset in a cold heading t 
which provides a little cone-shaped head with a base of about 22 mils as 
shown in Figure 13.7 for 6-mil diameter wires. The head is carefully pr 
tinned, leaving a small globule of solder at the end. Depending upon t 
size of the crystal plate, globules of 1000, 3000 or 7000 cubic mils of solder 
are used. The attachments are made in a wire soldering tool which attac! 
the wires to both sides of the plate simultaneously. This tool is illustrat: 
in Figure 13.8. The prepared wires are fed into positioning guides, whi 








per REE 










PLATED QUARTZ CRYSTAL UNITS 








heade 
cedure 


1g too 





mils 





on t 


250 BELL SYSTEM TECHNICAL JOURNAL 


are aligned with respect to the crystal plate, which is properly locate 

held between sliding jaws as shown. Prior to the operation the silver 

are burnished and fluxed. The wires in their guides are then slid int 

tact with the plate. The hot blasts are raised from below so as t 

directly at the work. The solder is melted and the attachment com 
in about 10 seconds when the bjasts are withdrawn. Slight press 

maintained on the wires during this operation by springs in the guid 
force the head of the wire to seat on the spot. A small fillet is obt 
around the head making the solder cover an area of about 35 to 40 
diameter. 

Crystal units using headed wire attachments have many advantages 
those made with straight or hooked wires. The puil-off strength is 
uniform and averages slightly better than that of hooked wires, despit 
fact that only a fraction of the amount of solder used with hooked wires 
employed. With this type of attachment the cratering effects encount: 
with previous methods have also been eliminated. The reduced quantit: 
solder on the face of the crystal plate effects a decided improvement ir 
temperature coefficient of the crystal unit as well as in its efficiency 
stability. Although the heading of the wires and the subsequent clea: 
and tinning operations involve more work, the process of making the atta 
ments is simpler and quicker, since the use of individual wires is bett 
adapted to making all the attachments in one operation. Headed wi 
attachments are more uniform in size and shape and give a more workn 
like finish to the job. This type of attachment has now replaced t! 
using straight and hooked wires in virtually all designs of telephone typ 
crystal units using wire supports. While the potential advantages of 
headed wire type of attachment for crystal support wires had been know: 
for many years, the practical exploitation of the idea depended on finding 
commercial means for producing the headed wires. The development of 
suitable machine for this purpose was carried out by the Western Electr 
Company in close collaboration with the Laboratories. Figure 13.9 sho 
such a machine. The fine wire is fed through the lower mechanism to 
die in which it is firmly clamped with a predetermined amount extendi 
above the plate. This part of the wire is then cold-worked by multi 
punches in the head of the machine until a conical shaped head is form: 
in the die cavity. As the individually headed wires are formed they 
cut off to a definite length and expelled as the vise is released and the next 
wire brought into position. Cold heading of the wires is necessary in or 
to retain the elastic properties of the phosphor-bronze springs employed 


the suspensicn. ‘The operation of the tool is simple after the precise alis 


ments of the die and punches have been made. 
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13.46 Mounting the Crystal Plate 


After the suspension wires are affixed to the plate, they are then bent 
serve as springs and to permit soldering into the cages as illustrated 
Fig. 13.2. Two different types of springs are used, one of them invol\ 
one bend and the other two. The direction of the bends and the dista: 
between them have been worked out so that the crystal will be displa 
to about the same extent in all three directions for equal forces. 1] 
cages are of simple construction being made up of mica stampings ; 


it 
l 


metal rods. The assembly of these parts is performed by welding 


AMPLITUDE OF WIRE VIBRATION I 
VERSUS 
NUMBER OF COMPLETE EXCURSIONS FOR 0063 INCH DIA. PH. BR. WIRE 
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Fig. 13.10—Characteristic performance of phosphor-bronze spring wires. 


eyelets, which are staked into the micas, to the rods. In the structures 
shown, the inside micas are provided with rectangular slots which limit the 
sidewise movement of the crystal plate from 25 to 30 mils. The end micas 
are spaced so as to limit the movement of the plate in the lengthwise direc- 
tion by the same amount. Aside from being used as parts of the cage, the 


micas therefore serve as “bumpers” to prevent excessive displacement, and 
possible breakage of the wires or plate if the crystal units are subject to 
extreme vibration or shock. Figure 13.10 is an experimental curve showing 
the minimum number of excursions made by wire-mounted crystal plates 
vibrated at different amplitudes before wire failure occurs for 6.3 mil phos 
phor bronze spring wires with single bends. On the basis of these data, th 
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1osen spacing of 25-30 mils between the crystal plate and the bumper 


ould ensure against any service failure of the unit in this regard. In 
der to center the crystal laterally and longitudinally in the bumper system, 
he plate is assembled first in the cage by means of spacers. The fine wires 
ire then soldered to the vertical rods or ‘‘straights” as they are called, and 
the spacers removed leaving the plate suspended in position. In order not 
to set up any strains in the junctions of the wire to the straight which might 
tend to displace the plate after this operation, the spring wires are usually 
pre-formed to come within about 5 to 20 mils of the straight. The junction 
is made by immersing the intersection of the wire and straight in a ball of 
molten solder. As the wires are withdrawn the ball of molten solder comes 
with them, solidifying in the air and thus joining the fine wire to the straight 
without strain. 

It will be noted from Fig. 13.2 that the wires to the longitudinal crystal 
are equipped with little weights close to the plate. This practice has been 
found desirable on virtually all types of crystal units to alleviate problems of 
wire resonance® which arise in occasional units thereby causing hig! -esist- 
ance as well as a shift in the frequency of the plate. Initially, while these 
effects were noted to some extent in the course of laboratory developments, 
it was not thought that they would be prevalent enough to warrant taking 
precaution to eliminate them by loading the wires, since they can usually 
be corrected by refloating and resoldering the crystal plate thereby changing 
the effective length of the wire. However, it has turned out that in manu- 
facture a large enough percentage of crystals contain resonant wires to 
warrant the use of weights. For low-frequency crystals (up to 200 kc) solder 
balls are placed on the wire at the desired location using a method worked 
out in conjunction with the Western Electric Company. The process is 
performed in somewhat the same manner as that described above for con- 
necting the crystal support wires to the straights, except that the weight 
of the solder deposited and the distance from the plate is more critically 
controlled. For higher-frequency crystals above 200 kc in which more pre- 
cise positioning of the weight is essential, small metal discs are employed. 
They are threaded onto the mounting wire and held in the correct position 
by a definite amount of solder on the back to obtain the desired loading. 
Since the free length of wire must be accurately controlled, the manufac- 
turing aspects of this job have been greatly simplified by the use of headed 
wires in which the variation in height of the solder cones is very small. The 
chart shown in Fig. 13.11 shows the weights of solder balls or discs and the 
position they should take on crystals having frequencies up to about one 
megacycle. It should be noted that the chart covers .0063”” phosphor 


*“Principles of Mounting Quartz Plates,” R. A. Sykes, B.S.T.J., April 1944 
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bronze wire. For any other diameter, d, of phosphor bronze wire, the ney 
distance 
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Earlier, it was menticned that the supporting wires for the plates wer: 
formed with one or two bends. In addition to the function of suspensio: 
these bends also introduce changes in impedance along the wire thus mini- 
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Fig. 13.11—Graph for determining placement of weights on wires for damping vibration 


mizing the possibility of trouble due to wire resonance. The use of a greater 
number of bends in the wire would tend to accomplish the same result as 
that of weights. However, the use of weights is considered more practical 
and has been adopted. As a result of this change, it is possible to employ 
wire supports having only one bend in virtually all crystals. In low- 
frequency crystals (below 2 kc) where the wave-length of the flexural wave 
in the wire is relatively long it is unnecessary to use weights since the wire 
length can be controlled adequately by the termination of the support wire 
at the straight. Depending upon the frequency, the desired length of wire 
is obtained by using either two or three direction bends. 





PLATED QUARTZ CRYSTAL UNITS 


13.47 Housing of Crystal Units 
For the pressure-type units first discussed no provision was made for 
protecting or sealing them other than the hermetically sealed containers in 
which all the other associated components of the filter were enclosed. How 
ever, the wire-supported designs have been worked out so that each unit 


is sealed in its own individual container. Fortunately, the sizes of virtually 


all crystal units are in the range which permits the use of relatively inexpen 
sive radio tube parts for these housings. There are many obvious advan 
tages to the individually sealed unit. After adjustment and sealing it can 
be handled more readily in subsequent assembly operations. It is not 
subject to variations due to changes in ambient humidity and consequently 
does not restrict the assembly of apparatus to conditioned space. It can 
be made up and stored or shipped as an individual unit. It has a higher 
degree of stability. There is one small effect, however, in the case of units 
which are sealed in vacuum. Due to the absence of any gaseous medium 
around the crystal, a slight change in frequency is encountered when the 
tube is evacuated. However, this change is always the same for each 
particular type and size of crystal and can be allowed for in the final adjust- 
ment before sealing. 

Most designs of crystals can be sealed in an atmosphere of dry air although 
better performance results from the use of vacuum. Some crystals must 
be sealed in vacuum for this reason. A decided advantage in favor of 
vacuum-sealed crystals is the elimination of acoustic effects from air 
resonance. 

Both metal and glass tubes are used for housing crystal units. Initially 
it appeared that metal tube radio parts were ideally adapted to crystal use, 
and it was felt that, instead of welding the stem to tube, this sealing opera- 
tion could be done by soldering. However, it was found that while sound 
solder joints could be obtained, extreme precautions were necessary to 
protect the button-type glass seals, through which the leads emerge, during 
the pre-tinning and soldering operations. Even with such precautions, it 
would have been essential to include in every vacuum type tube a means of 
detecting whether or not a leak had developed. For air-filled tubes at 
atmospheric pressure this would not have been necessary since minute leaks 
can be tolerated with little likelihood of the crystal being affected over a 
long period of time. The possibility of welding as is done in the case of 
radio tubes was considered but did not appear justified on the basis of 
equipment cost. Moreover, even with welding there still appeared to be 
problems from leakage and outgassing of the metal since, after the crystal 
is enclosed, the assembly cannot be exposed to high temperature to drive 


off adsorbed gases during the evacuation process. In view of these draw 
’ 
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backs the use of metal tubes has been discarded in favor of glass tubes excep‘ 
in the case of a few special designs. 

The procedure of mounting a crystal unit on a stem and sealing it ir 
glass is much the same as fora radiotube. Figure 13.2 shows crystal unit: 
mounted on stems ready for sealing and also shows units sealed in glass and 
based. The extensions of the straights through the bottom micas are welded 
to the formed wires emerging from the glass seals. In the glass-sealing 
operation care must be taken not to heat up the assembly to the point wher« 
the solder attachments will be melted or even softened enough to permit the 
crystal to change position. To accomplish this it is necessary to use hot, 
sharp-pointed fires localized to the region where the seal will be made. Thx 
use of oxygen-gas flames is virtually essential to accomplish the seal quickly. 
Having the fires strike the bulb at tangency is also desirable. The ordinary 
type of glass-sealing head for use with gas-air fires is not well adapted to 
this work since the rotating pillars require the fires to be held too far away 
from the work thereby necessitating larger flames and consequently more 
heating up of the crystal unit assembly. The screening effect of the pillars 
as they revolve also slows up the work of the fires thus increasing the over-all 
heating of the assembly. A special glass-sealing machine developed for 
sealing crystal units is shown in Fig. 13.12. Immediately following the seal- 
ing operation the glass units should be placed in a suitable annealing box 
or leer where they can cool off very slowly. A large wooden block equipped 
with holes to admit the individual bulbs is convenient. The holes may be 
covered with a cloth to prevent air circulation. 

After the units have cooled they are placed on a vacuum pumping station 
and evacuated. During the first half hour of pumping they are enclosed 
in a heated oven in which the ambient temperature is maintained at about 
240°F. This drives off any traces of moisture that might have entered the 
tube prior to sealing. Following the heating interval, the tubes are pumped 
for another half-hour during which time they will have cooled down to room 
temperature. At this point the pressure in the tubes should be at the 
minimum of which the pump is capable of attaining. This value should be 
at most 20 microns and preferably less. However, with a six. or eight-tube 
station better than 15-20 microns is not likely to be attained unless a liquid 
nitrogen trap is employed in the system for eliminating moisture. 

After the pumping period, vacuum-tyre units are sealed off, with pump 
running, by melting the glass tubulation with a fine-pointed oxygen-gas 
flame as close to the stem as possible. If air is to be admitted, the pump is 
closed off from the system and dry air admitted to the tubes after which 
they are sealed off. After testing the crystal unit to see that it meets its 
requirements, the unit is equipped with a base in the same manner as fol 
lowed for radio vacuum tubes. 
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Fig. 13.12—Glass-sealing machine. 
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13.48 Stabilization of Crystal Units 


Despite the close dimensional tolerances applying to the manufactur: 
the individual crystal plates, the exact frequencies are rarely realized in 
mounted crystals. To bring the crystal to frequency it is therefore ne 
sary to grind off minute layers from the lengthwise or widthwise edges 
the plate depending on the mode of vibration. This adjustment, whi 
causes superficial disruption of the quartz areas affected, results in unsta! 
operation of the unit with respect to frequency and resistance. Unless t 
crystal plate is properly treated after these operations, considerable dri 
in these characteristics will take place, particularly during the initial servi, 
life of the unit To alleviate this condition, the crystal units are first roug! 
adjusted to the approximate frequency and then heat-aged in an oven whi: 
subjects the units to several heating and ceoling cycles between 240°F a: 
75°F. The units are then mounted in their cages as previously describe: 
and fine-adjusted after which they are again aged. This operation als 
tends to drive off any moisture which might be troublesome. With this 
type of accelerated aging the crystals are stabilized to the point wher 
changes in performance can be detected only by the use of the most precis: 
measuring equipment over long periods of time. Crystals so stabilized may 
generally be depended upon, at any one temperature, maintaining thei: 
frequency indefinitely within two or three parts per million provided they 
are not subjected to excess voltage. 


13.49 Cleaning of Crystal Units 


Throughout the manufacturing process it is essential that every precau 
tion be taken to keep the crystal plate and associated parts absolutely fre 
from contamination and dirt. The rigorous cleaning necessary before the 
spotting operation has already been discussed. In all the subsequent opera 
tions care must be taken to prevent the plates coming in contact with sul) 
stances that might tend to cause corrosion. Any particles of foreign matte: 
that may have accumulated on the plate or wires before rough and fin¢ 
adjustments should be carefully washed off. Otherwise the performanc« 
and life of the crystal may be adversely affected. A suitable method for 
cleaning crystal units before sealing consists of washing and rinsing i! 
chemically pure carbon-tetrachloride or other suitable solvent to remov« 
grease, followed by washing and rinsing in hot distilled water at about 
150°F. To facilitate removal of unwanted substances, the parts should b« 
scrubbed gently with a soft brush or agitated in the solution during this 
process. The use of pure alcohol (95%) in addition to carbon-tetrachloricd« 
is also good for this process, but is not essential. The cleansed crystal units 


should be carefully dried out and protected from further contaminatio: 


prior to the sealing operation. 
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13.5 CONCLUSION 


In ending I should like to acknowledge the aid and useful suggestions 
given me by Mr. C. E. Lane and my other associates in preparing this 
article. I should also like to reiterate the fact that the status of the art as 
described was reached after many years of pioneering development by many 
engineers. In some cases the names of individuals associated with specific 


contributions of a major nature have been mentioned. 
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CHAPTER XIV 


Effects of Manufacturing Deviations on Crystal 
Units for Filters 


By A. R. D’HEEDENE 
14.1 THe Errect oF DEVIATIONS IN THE CHARACTERISTICS OF CRYS1 
UNITS ON FILTER PERFORMANCE 


HIS chapter emphasizes primarily the need for close control in the 

manufacture of crystal units for use in filters. The first telephon 
use of crystal units in the commercial manufacture of filters was made by 
the Western Electric Company in about 1936. To make such commercial 
manufacture practical, it was necessary to establish accurate design informa- 
tion and allowable manufacturing tolerances. The quantitative data 
collected for this purpose provided the chief source of material for this 
chapter. While the data is quite extensive, it will be observed that there 
are still some factors which must be treated qualitatively. 

While filter crystal units are like oscillator crystal units in that they must 
have low internal dissipation and a close control of resonant frequency, they 
are different in that many additional characteristics of the filter crystal 
units must also be controlled accurately. Two typical illustrations will 
demonstrate how characteristics other than resonant frequency and Q may 
react on filter performance. 

The first characteristic considered is the slope of the reactance with 
frequency curve in the vicinity of the series resonant frequency. This slope 
is sometimes referred to as the impedance level of the crystal unit. A con 
venient measure is the inductance of the equivalent electrical circuit 
When this inductance departs frem its ncminal value, the performance of 
the filter using the crystal unit may undergo appreciable change. This is 
particularly true of filters in which the schematic contains a lattice or some 
other type of bridge circuit with crystal units contained in all the bridge 
arms. For example, in Fig. 14.1 the solid curve illustrates the transmission 
characteristic obtained from a lattice-type crystal filter, in which both the 
series branches and the diagonal branches contain two balanced crystal 
units. High loss results from a close impedance balance between the 
branches of the lattice. When the inductance of any of the crystal units 
departs from its nominal value, the bridge balance is disturbed and the 
transmission characteristic of the filter is changed. The two dotted curves 
of Fig. 14.1 illustrates the characteristics that result when the inductance 
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values of the crystal unit in either branch depart from their nominal values 
by about one percent. A negative departure in one branch results in about 
the same effect on performance as a positive departure in the other branch. 


The difference between the two curves shown on Fig. 14.1 is that one as- 
sumes a positive departure and the other a negative departure for the in- 
ductance of a branch. 

Due to the close impedance balance which is required for these filters, 
the effect of small departures in resonant frequency will produce rather large 
variations in the transmission characteristic. For example, departures of 
about 10 cycles per second in the crystal units of either branch will produce 
variations in discrimination of about the same type and magnitude as those 
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Fig. 14.1.—The insertion loss characteristic of a crystal band-pass filter as affected 
by deviations in the inductance of the crystal units. 


illustrated in Fig. 14.1 for departures in inductance. On the other hand, 
if the crystal units of both branches exhibit equal departures the entire 
transmission characteristic will be shifted by the frequency departure of the 
crystal units, and there will be no loss in discrimination. 

Another way in which deviations in the properties of crystal units may 
react on filter performance is illustrated by the schematic and curves shown 
in Fig. 14.2. The schematic is the equivalent electrical circuit of a narrow 
band filter, using two balanced quartz crystal units. The filter is designed 
to provide a passed band of about 10 cycles per second with distortion of 
less than 0.2 db. The insertion loss characteristics show that the desired 
transmission can be obtained for various magnitudes of effective resistance 
as long as the resistances in the series and diagonal branches are equal. 
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However, if the effective resistance in one branch is twice as large as that i: 
the other branch a highly distorted characteristic results as shown by th: 
middle curve of Fig. 14.2. 

Both of these illustrations show that filter performance is degraded rapidl) 
if the crystal units of the lattice have characteristics which depart from 
their nominal values by different extents for the two branches. A similar 
effect is produced when the temperature coefficient of resonant frequency 
for the crystal units in one branch differs frem the temperature coefficient 
of the units in the other branch. Deviations occurring in a single unit may 
also affect filter performance. Such deviations include the presence oi 
unwanted resonances of even weak amplitude, inadequate insulation re 
sistance between the metallized coatings or unbalance between the halves 
of plates on which the coating has been divided. 
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Fig. 14.2.—Effect of deviation in the effective resistance of crystal units on the distor 
tion characteristic of a crystal filter. 


The importance of controlling the electrical characteristics of the crystal] 
units is indicated from the above considerations. It is pertinent to correlate 
deviations in the mechanical properties of the crystal unit with the devia 
tions in electrical characteristics. This is the subject of the succeeding 
sections. Consideration is restricted to the plates commonly used in filters, 
that is, X-cut plates, vibrating in extensional or flexural modes, and GT-cut 
plates. 


14.2. Tue Errect oF DEVIATIONS THAT OCCUR IN THE MANUFACTURE OF 
QUARTZ PLATES 


Quartz is an anisotropic material. Accordingly, plates cut from a quartz 
crystal exhibit elastic and piezo-electric properties which depend on the 
orientation of the plates with respect to the principal axes of the crystal. 
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For that reason, any deviation in the orientation of the plates from nominal 
will affect the electrical characteristics of the crystal units. In addition, 
these characteristics are affected by imperfections in the plates due to 
deviations in linear dimensions, to the presence of flaws, or to the condition 
of the surface of the plates. The effects of these deviations differ for various 
cuts of crystal plates, for plates of various chapes ard for the various modes 
of vibration. In the following paragraphs, each type of deviation will be 
considered in turn and data will be presented to show its effect on the charac- 
teristics of crystal units using the various types of plates. 


14.21 DEVIATIONS IN THE ANGLE OF ORIENTATION 


Accurate information is available on the effect of deviation in angle of 
orientation on the characteristics of X-cut plates vibrating in the exten 
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Fig. 14.3.—Equivalent electrical circuit of piezoelectric crystal. 


sional mode. The relation between the electrical characteristics of this 
type of vibration and the properties of the quartz are shown in Fig. 14.3 
This information, with minor changes, is reproduced from a preceding 
publication.'. In Fig. 14.3: /,w and ¢ are the length, width and thickness 
respectively of the plate; A is the dielectric constant; p is the density; 
diz is the piezo-electric constant; and Se is the modulus of compliance 
(inverse of Young’s modulus). All these individual quantities are ex- 
pressed in electrostatic units. The quantities which cepend on the orien- 
tation of the plates are the piezo-electric constant and the modulus of com- 
pliance. The symbols for these quantities usually are primed when they 
are used for a generalized orientation. When unprimed, the symbols desig- 
nate quantities measured along the principal axes. For X-cut plates, devia- 
tions of the plane of the major surface from the YZ plane have relatively 


1 “Electrical Wave Filters Employing Crystals with Normal and Divided Electrodes” 
W. P. Mason and R. A. Sykes, B. S. T. J., April 194), page 222 





| 
| 
| 


264 BELL SYSTEM TECHNICAL JOURNAL 


small effect, while variations in the angle of rotation about the X-axis hav: 
a relatively large effect on these quantities. 

Mason has shown? how the magnitudes of the piezo-electric constant: 
and the moduli of compliance for any angle of rotation may be derived fron 
their magnitudes along the principal axes of quartz. Using these equation 
and the magnitudes for the principal axes tabulated in a recent paper? by 
Mason, dj2 and $92 have been calculated as a function of the angle of rotatio: 
of the plates about the X-axis. In turn, the frequency and inductance con 
stants have been calculated as a function of the angle of rotation, using th« 
relations shown in Fig. 14.3. Figure 14.4 is a plot of the frequency and in 
ductance constants as a function of the angle of rotation for angles between 
about —70° and +70°. It shows how the inductance and resonant 
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ANGLE OF ROTATION 


Fig. 14.4—Frequency and inductance constants of X-cut quartz plates as a function 
of their angle of rotation around the X-axis. 


frequency of these plates will change with deviations in the angle of rotation 
The frequency constant used is the product of the resonant frequency in 
kilocycles and the length of the plate in millimeters. The inductance 
constant used is the inductance per millimeter of thickness of a plate having 
a width-to-length ratio equal to 0.1. 

The change of inductance and resonant frequency with deviations in 
angle of rotation about the X-axis is of particular interest for the angles of 
rotation most commonly used, that is for —18.5° and for +5°. The 
calculations indicate that fr an X-cut plate rotated —18.5°, a devia- 
tion of +1° in the angle of rotation will change the inductance by +1.2% 


2 “Electrical Wave Filters Employing Quartz Crystals as Elements”, W. P. Mason 
B.S. T. J., July 1934, equations on page 451; also in Chapter VI. 
3 “Quartz Crystal Applications,” W. P. Mason, B.S. T. J., July 1943. 
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respectively, and the resonant frequency by +.05% and —.02%, respec- 
tively. Foran X-cut plate rotated +5°, a deviation of +1° will change the 
inductance —0.9%, and +0.6%, respectively, and the resonant frequency 

Deviations in angle of rotation will, in general, affect temperature coeffi- 
cient. The effect is illustrated by Fig. 1.19 of Chapter I,* which, for X-cut 
plates, shows the relation between temperature coefficient and angle of 
rotation. Ths curve shows that the temperature coefficient is practically 
zero for an angle of rotation of +5°. For that reason this particular cut is 
used whenever a low-temperature coefficient is desired. The curve also 
shows that at this point the slope of temperature coefficient as a function of 
angle of rotation is zero. Hence, for the +5° X-cut plate, which is most 
important from the standpoint of temperature coefficient, there will be 
little change due to a deviation in the angle of rotation. 

In GT-cut plates the effect of deviation in the angle of orientation must 
be considered in combination with deviations in linear dimensions. Mason 
shows! that for an angle of rotation of +51 degrees 7.5 minutes and a 
width-to-length ratio of .859, a temperature coefficient close to zero may be 
obtained from —25° C to +75° C. He also has shown that this tempera- 
ture coefficient varies with both the angle of rotation and the width-to- 
length ratio. Because of this, it has been found possible to compensate for 
small deviations in the angle of rotation by adjusting the linear dimensions. 
The net effect of a deviation in angle of rotation, after it has been so com- 
pensated, is to raise (or lower) the temperature region of zero temperature 
coefficient by 11° C for each 10-minute increase (or decrease) in angle of 
rotation. In GT-cut plates, the width dimension directly controls the 
primary resonance. For this reason, it is preferable to adjust temperature 
coefficient by varying the length dimension rather than the width. The 
crystal plates are cut larger than desired. The frequency of resonance and 
the temperature coefficient are then adjusted simultaneously by grinding 
either the width or length dimension as required. Experimental work 
carried on by L. F. Willey of the Laboratories shows that an increase of 1.0 
per cent in the width-to-length ratio results in an increase in temperature 
coefficient of approximately +1.35 parts per million per degree C. In his 
experimental work Willey used the ratio of the secondary to the primary 
frequency as a convenient measure of the width-to-length ratio. The in- 
ductance constant for GT plates, which is about 17 henries per millimeter 
of thickness, will increase by less than 1% for deviations in any of the 
angles of rotation of as much as 30 minutes. However, the inductance may 
depart appreciably from nominal due to the adjustment of width and 
length dimensions. ’ 


4“New Quartz-Crystal Plate, GT, Produces Constant Frequency Over Wide Tem- 
perature Range”, W. P. Mason, Proc. I. R. E., May, 1940, page 220. 
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14.22. DEVIATIONS IN LINEAR DIMENSIONS 

In the case of X-cut plates, the length dimension is used to control th 
location of their resonant frequency. The length is lapped to its final dimen 
sion after all other processes have been completed, so that this dimensior 
will be such as to compensate for the effect of any other deviations that may 
have occurred. The sensitivity of this adjustment depends on the mode of 
vibration. For plates vibrating in their extensional mode, the resonant 
frequency is inversely proportional to the length, as shown in Fig. 14.3, 
while for plates vibrating in their fexural mode, the resonant frequency is 
inversely proportional to the square of the length. The amount of the ad 
justment required depends on the magnitude of the frequency errors that 
may have been introduced due to deviations in the width or in the angular 
orientation of the plates or due to still other causes. The magnitude of 
such frequency errors, in turn, depends to a considerable degree on the angle 
of rotation of the plates. For example, it was shown in Section 14.21 that 
a deviation in the angle of rotation of a +5° plate changes its resonant fre- 
quency more than ten times as much as a similar deviation in a —18.5 
plate. It must be noted that the adjustment of length compensates for 
frequency errors only and that errors in inductance or temperature coefficient 
may be increased by such adjustment. 

Deviations in the thickness dimension principally affect the impedance 
level of the plates. As shown by Fig. 14.3, the inductance is directly and 
the capacity inversely proportional to the thickness. In the case of GT-cut 
plates the thickness dimension is important also because it controls the 
location of the most prominent unwanted resonances, which arise from 
vibrations in thickness flexure. However, plates are designed to avoid 
critical thicknesses and small deviations from the nominal thickness will 
not usually result in plates having unwanted resonances. 

Deviations in the width dimension affect the equivalent electrical charac- 
teristics appreciably. The effect of deviations in width on the frequency 
of X-cut plates vibrating in their extensional mode can be deduced from 
the curves of Fig. 14.5. The curves show that this effect is more pro- 
nounced for larger values of the width-to-length ratio where coupling with 
the width extensional mode becomes appreciable. For a —18.5° plate with 
a width-to-length ratio of 0.8 an increase of 1% in the width dimension 
will decrease the frequency by about .04%. For a +5° plate with a 
width-to-length ratio of 0.4 an increase of 1% in the width dimension 
will also decrease the frequency by about .04%, but for a ratio of 0.6 the 


decrease in frequency will amount to 0.13%. 
Similar information is available for crystals vibrating in width flexure 
from the measurements published by Harrison’ and the calculations pub- 


5 “Piezo-Electric Resonance and Oscillatory Phenomena with Flexural Vibration in 
Quartz Plates”, J. R. Harrison, Proc. I. R. .©., Dec. 1927. 
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lished by Mason.* This information shows that for small ratios of axes 
the resonant frequency will be directly proportional to the width dimension. 
However, as the ratio is increased to 0.5, a change of 1° in the width dimen- 
sion will change the frequency by only 0.5%. 

The effect of the width dimension on the inductance of the plates fre- 
quently is important. Fig. 14.6 illustrates the relation between inductance 
and the ratio of axes. From these curves the effects of deviations in width 
can be deduced. For the two longitudinal plates, the inductance is almost 
inversely proportional to width. For the flexure plate, the decrease of 
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Fig. 14.5.—Frequency constant of the longitudinal mode of X-cut quartz plates as a 
function of their ratio of width to length 


inductance with increase in width is much more rapid. With a ratio of 
axes of 0.6 the inductance decreases about as a square power, while with a 
ratio of 0.1 the decrease is about as the third power. 

The width dimension of the +5° plate has an appreciakle effect on the 
temperature coefficient of the plate. Mason has shown® that while the 
temperature coefficient is zero for a long narrow bar, it increases quite 
rapidly as the width dimension increases, due to coupling between the 
face shear and the longitudinal modes. In the case of an —18.5° plate, 
coupling with other modes is relatively weak. Hence its temperature co- 


§ “Motion of a Bar Vibrating in Flexure Including the Effects of Rotary and Lateral 
Inertia’, W. P. Mason, Jour. Acous. Soc. America, April, 1935, pages 246-249 
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efficient, which is about 25 parts per million per degree C, does not change 
appreciably with changes in width. For a +5° plate vibrating in its 
flexure mode, Fig. 14.7 illustrates measurements made on the variation o! 
temperature coefficient with ratio of axes. For all these X-cut crystals, 
it may be observed that deviations of 1% in the width dimension will not 


10000 f re 
Bee 






| 
| 
ke a ames: + $4 + + 4 














- HENRIES PER MM THICKNESS 


INDUCTANCE 





8.91 


RATIO OF WIDTH TO LENGTH 


Fig. 14.6.—Inductance of the crystal units used in filters as a function of the cuts of the 
plates and their ratio of width to length. 


change the temperature coefficient by more than 5%. Such changes 
are usually negligible. 


14.23 INTERNAL DEFECTS 


Internal defects in the quartz plates may have a large effect on their 
electrical characteristics. These defects vary so widely in type, size and f 
concentration that it is impossible to predict the effects quantitatively. 
General comments regarding the results that may be expected for various 
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defects are described in Chapter IV’. The conclusions drawn there for 
oscillator plates are also applicable to filter plates. These are: (1) Evidence 
that a particular defect is permissible in a given type of plate does not prove 
that a similar defect is permissible in scme other type of plate, and (2) 
proof that a particular defect is permissible in a given type of plate can be 
obtained only by a statistical study. 

Scme qualitative statements can be mace regarding the effect of mechan- 
ical flaws. Cracks result in instability of resonant frequency and effective 
resistance and must ke avoided. The effect of inclusions or chips depends 
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Fig. 14.7.—Temperature coefficient for +5° flexural crystal units as a function of the 
ratio of width to length of the plates. 


on the size of these defects relative to the size of the finished plates and also 
on their location in the plate. 

Twinning in quartz may be either of the optical type (Brazil twin) or of 
the electrical type (Dauphiné twin)’. The effect of these two types of 
twinning on the performance of oscillator crystal units has been described 
thoroughly in Chapter V°. 

When optical twinning is present, the plate will exhibit the same elastic 
properties throughout, but the two portions of the plate will tend to expand 

7“Raw Quartz, its Defects and Inspection”, G. W. Willard, B. S. T. J., October 1943, 
pp. 338-361. 

8 “The Properties of Silica”, R. B. Sosman, Chapter XTT. 


®“Use of the Etch Technique for Determining Orientation and Twinning in Quartz 
Crystals”, G. W. Willard, B. S. T. J., January 1944, pp. 11-51. 
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and contract in opposite phase. Hence there is little change in frequenc) 
constant or temperature coefficient, but there will be a large change in i: 


ductance. Th change in inductance can be estimated roughly by compar 


ing the twinned plate with an untwinned plate in which activity has bee: 
reduced by removing electrical charge from part of the surface. The are: 
of the surface from which this charge is removed would be twice the twinne: 
area and located at about the same position in the plate. 

It is believed that a small amount of electrical twinning is more serious 
than a similar amount of optical twinning, because the twinned areas ar 
of opposite angular sense. Each of the two areas has a different modulus 
of compliance and the effective modulus of the plate has a value intermediate 
between the two different values of modulus. Therefore, the frequenc\ 
constant of the plate will be intermediate between that of the desired «1 
and its electrical twin. For a small amount of twinning, the direction an 
rate of change of frequency can be estimated from the comparison shown 01 
Table I between the standard filter cuts and their electrical twins. 


TABLE I 


Frequency 
Constant—kc. m.m. 


Frequency 


Electrical Twin Constant—ke. m.m 


Filter Plate 


2560 = 3120 
2815 2650 
3280 a 2610 


This verifies the experimentally observed fact that for —18.5° X-cut 
plates, twinning increases the frequency, while for +5° X-cut and GT plates, 
twinning decreases the frequency. Even for small amounts of twinning 
the inductance will increase rapidly for plates of any orientation. Whe: 
the amount of twinning becomes large, the equivalent inductance approaches 
infinity. That is, the crystal will not be set in motion by an applied voltage 

The quantitative effect of twinning (probably electrical) has been meas 
ured on one set of plates by R. M. Jensen. Figure 14.8 includes a 
photograph of the plates used, illustrating the extent of the twinning i 
each. All of the plates are —18.5° X-cut plates, having the dimensions 
30.88 x 10.56 x .86 mm. The tabulation below the photograph compares 
the inductance and resonant frequency measured for each of the plates 
with the one, designated AN-3, which shows the least effect of twinning 
While there is a good correlation between the amount of twinning in th: 
plates and their change in electrical performances, it is not practical t 
estimate accurately the effect of a given amount of twinning. For this 
reason, crystal plates having any twinning should not be used for crystal! 
units for filters. 
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14.24 EtTcHInc 


The surface condition of the quartz plates also has some effect on crystal 
characteristics. This surface condition is determined in large part by the 
apping operation used to obtain final dimensions. As described in Chapter 


K-cut 
lates, 
: Percentage Increase over Values 

nning Measured for AN-3 

: Plate Resonant 
Wher Designation Inductance Frequency 
aches AN-1 +22.12 + .50 
AN-2 +27.20 +.59 
tage AN-3 0 0 
meas AN-4 +3.03 + .04 

; AN-5 +3.12 + .01 

les a AN-6 +276.54 +5.01 
ng ii AN-7 +1.58 — .03 

AN-8 +14.21 +.37 
IS1ONS AN-9 +738.00 +6.63 
pares AN-10 +32.44 + .83 

I Fig. 14.8.—Effect of various degrees of twinning on the performance of —18.5° X-cut 
plates juartz crystal plates. 
ning 

n the XIII", the plates are given a final lap with 400 or 600-mesh carborundum. 
cal ti This, in turn, is followed by an etching bath which removes foreign particles. 


© this A short etch, about eight minutes in 47% hydrofluoric acid, has been found 
rvsta adequate to ensure firm adherence of the metal coating to the quartz. On 


10“The Mounting and Fabrication of Plated Quartz Crystal Units,” R. M. C. Green 
idge, this issue of the B. S. T. J. 
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the other hand, the use of a relatively long etch, 30 minutes or more, i 
desirable when a high Q is desired. The long etch also results in an im 
proved stability of the resonant frequency as a function of current. Thi 
will be discussed in a subsequent paragraph. A disadvantage of a lon; 
etch is the difficulty of controlling the etching process within close toler 
ances. The variations in rate of removing material may be sufficient t: 
affect the uniformity of the linear dimensions of the plates. 

These factors indicate that etching is an important process in preparing 
crystal plates. A close control must be maintained on the strength of th 
acid, the uniformity with which the surfaces of the plates are exposed and 
the duration of the exposure. 


14.3 THe EFFEcTs OF DEVIATIONS DURING FABRICATION OF WIRE- 
SUPPORTED UNIT 


As described in Chapter XIII'®, two types of mountings have been de- 
veloped for supporting crystal plates, the Pressure Type and the Wire 
Supported Type. The wire-supported type of mounting is the more recent 
development and has resulted in crystal units which have a much higher 
degree of stability and can be reproduced within much closer tolerances 
than the units using the pressure type of mounting. Since this chapter is 
concerned chiefly with the problem of obtaining a high degree of precision 
in crystal units, the discussion is restricted to the wire-supported type of 
mounting. 


14.31 SitveR SPotTTING 


For the wire-supported type of mounting the first operation is to bake 
small silver spots on the surface of the crystal plates. In the application 
of these silver spots to the crystal plates three factors are of importance in 
their effect on the characteristics of the plate, namely, the size of the spot, 
its location, and the firing temperature. Since in all crystal designs to 
date the silver spots are applied at or near the nodal line of the crystal plate 
the principal effect of the spots is to increase the stiffness of the plate, so 
slightly increasing the frequency of resonance. Variations of an appreciable 
magnitude in either the amount of silver paste used (that is, the size of the 
spot) or in the location of the spot with respect to the nodal line will change 
the resonant frequency of the plate. Such changes could be corrected 
later, when the plates are adjusted for resonant frequency, »; long as the 
length is increased sufficiently to allow such adjustment. However, if the 
length be increased sufficiently to allow for extreme cases, ave. ige adjusting 
time will be increased materially, while if the allowance is insufficient some 
of the plates may be unusable. For this reason, close control of the size 
and location of the silver spots is well justified. 











im 
Thi 
lon; 
oler 
it ti 


| de- 
Vire 
cent 
cher 
nces 
er is 
sion 


Ee of 


ake 
tion 
‘e in 
pot, 
Ss to 
late 
2, SO 
able 
the 
inge 
cted 
the 
the 
ting 
ome 


size 








MANUFACTURING DEVIATIONS IN CRYSTAL UNITS 


In baking the silver spots, care must be taken to prevent “heat” twinning. 
If the temperature of a quartz plate is raised above the inversion point 
573°C) and then is reduced again, the plate will be electrically twinned.® 
[he firing temperature of the silver paste currently used for the spots is not 
many degrees below this inversion point. Hence, the firing temperature 
may easily become so high as to result in twinned plates. In addition, 
it has been observed that the twinning may occur at a considerably lower 
temperature if the plate is subjected to large thermal stress. For this 
reason, care must be taken to heat the plates uniformly during the baking 
operation. 


14.32 DrIvIsIon oF COATING 


The next operation is to evaporate a coating of silver on the surface of 
the quartz plates. The plates must be thoroughly cleaned before this coat- 
ing is applied in order to ensure firm adherence of the coating. Poor ad- 
herence may cause the coating to peel off the plate, changing all of the 
electrical characteristics of the plate. In many cases the coating must also 
be divided.'!! Two methods are in general use for dividing the silver 
coating on crystal plates, namely, an abrasive method and an electrical 
stylus method. In general, the abrasive method of dividing the coating is 
superior to the electric stylus for all cases requiring a simple straight line 
division, but it has not been found practical for complicated divisions such 
as are desirable for harmonic longitudinal plates and flexure plates. 

In using the abrasive method for dividing the coating only two factors 
are likely to change the characteristics of the crystal plate, these being the 
location and the width of the dividing line. Deviations in the location of 
the dividing line from the lengthwise center line for a longitudinal plate 
will affect the capacity and inductance balance between the two halves of 
the plate. Deviations in the width of a properly centered dividing line 
will cause changes in the inductance of the plate since for a given plate the 
inductance is a function of the ratio of the plated area to the total area of 
the plate. So, for a wide crystal plate deviations in the width of the divid- 
ing line will be negligibie while for narrow plates these deviations can 
cause an appreciable change in the inductance of the plates. 

When the electric stylus is used for dividing the coating, the location and 
the width of f ie dividing line again will affect the performance of the plates. 
In addition, varying amounts of twinning will occur along the division 
lineapparen { due to instantaneous high temperature gradients introduced 
by burning of the silver at the point of contact of the stylus. In measure- 


1 Loc. cit. 
8 Loc. cit. 
11 “Crystal Channel Filters for the Cable Carrier System”, C. E. Lane, B. S. T. J., 
January 1938, pp. 125-136. 
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ments made on a group of —18.5° X-cut crystals on which the coati: 
has been divided carefully with an electric stylus, the increase in the indu 
tance of the plates ranged from 1.4% to 2.6%. Any twinning resulti: 
from the dividing operation will also change the resonant frequency of t! 
plates. 


14.33 SOLDERING OF WIRES TO PLATES 


The next process, that of soldering the supporting wires to the crystal 
plate may have considerable effect on the performance of the unit. The 
deviations which may be introduced depend on the amount of solder used, 
the location of the solder button with respect to the nodal line of the plate, 
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Fig. 14.9.—Change of resonant frequency with temperature of +5° X-cut quartz 
crystal plates. The curves show that when the volume of solder used for joining the 
plates to the supporting wires is appreciable compared to the volume of the plates, the 
frequency-temperature coefficient characteristic is affected by the volume of solder. 


the shape of the solder button, and the possible twinning of the plate during 
the soldering operation. 

The amount of solder used in forming the joint of the wire to the plate 
becomes extremely important when the plate is small. For example, Fig. 
14.9 illustrates the changes in the frequency-temperature characteristic 
resulting from the use of varying amounts of solder for a particular size of 
plate. The units on which these measurements were made used X-cut 
+5° plates of 16mm x6mmx0.5 mm. The types of wire referred to in 
the figure, that is, hooked, straight and headed, were described in Chapter 
XIII". The frequency-temperature characteristic expected on the basis 
of measurements made on crystal units using larger plates is approximated 


10 Loc. cit. 
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losely by the solid curve. This solid curve actually was obtained from 
measurements made on crystal units using plates supported with headed 
wires and using a very small amount of solder. The other curves indicate 
that the temperature coefficient may be increased appreciably due to the 
presence of a larger amount of solder. Further, when the larger amounts 
of solder are used, the characteristics depend on the exact amount of the 
solder, so that the characteristics represented by the dashed curves are 
hard to reproduce. 

The amount of solder used in this operation also affects the Q of the crys- 
tal unit and its resonant frequency. Measurements using several crystal 
plates of relatively small sizes have shown improvements in Q of as much as 
25 per cent when headed wires are used over that obtained with other wires 
using larger amounts of solder. 

Variations in the consistency of the solder joint will, of course, affect the 
adherence of the supporting wire to the plate. A poor joint will result in a 
high effective resistance for the crystal unit and will generally cause in- 
stability both in resistance and in resonance frequency. 

In soldering the supporting wire to the crystal plate two methods have 
been used for melting the solder; namely, the soldering iron, and the hot- 
air blast. With either method, lack of sufficient control can seriously 
change the electrical characteristics of the plate due to twinning. It has 
been observed that this twinning occurs when there is a large temperature 
gradient in the quartz, even at temperatures well below the inversion point. 
Experimental work by G. W. Willard has indicated that it may occur even 
when the temperature of the soldering iron is as low as 300°C. To avoid 
such twinning during the soldering operation, it has been found desirable 
to raise the temperature of the entire plate to just below the melting point 
of the solder. 

Twinning, when it occurs, will affect the crystal plate by causing an in- 
crease in inductance, a change in the resonant frequency, increased effective 


resistance, and a change in the temperature coefficient, as stated previously. 
Also, in crystals with divided plating there will be an inductance unbalance 
between the two halves of the crystal plate set up due to unequal amounts of 
twinning. Several measurements made, using GT plates at 160 kc, showed 
that twinning during the soldering operation decreased the resonant fre- 
quency in a range from 200 to 100 cps and the temperature coefficient of the 
units ranged from 2 to 6 times that of units using untwinned crystal plates. 


14.34 EFFECTS DUE TO WIRE RESONANCE 


As described in Chapter VIII’, the characteristics of crystal units may 
be changed due to vibrations set up in the supporting wires. When any 


12 “Methods of Mounting and Holding Crystals”, R. A. Sykes, B. S. T. J., April 1944. 
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one of the wires is not located exactly on a node of the plate, the plate w 
set the wire into vibration. For certain critical lengths of the wire, it w 
offer consicerable resistance to this motion and there will ke a rapid increas: 
in effective resistance and some change in resonant frequency of the cryst: 
plate. 


The effect of wire vibration can ke cescriked in terms of its electric 
analogy. The vibrating wire, clamped at its far end, may be considered a 
rather special electrical transmission line open-circuited at its far end 
When viewed frcm the crystal plate the impedance changes rapidly with 
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Fig. 14.10.—Effect of wire resonance on the resonant frequency of a crystal unit. 


frequency in a succession of pronounced resonances and anti-resonances. 
In the vicinity of an anti-resonance the electrical equivalent of the vibrating 
wire may be approximated by a coil and condenser in parallel as shown by 
L, and C2 of Fig. 14.10. This acts in series with the mechanical resonance 
of the quartz plate, represented by ZL; and C;. The impedance curves illus- 
trate the effect of the wire resonance on the crystal impedance. R,, the 
equivalent resistance of the crystal plate, is constant for frequencies in the 
vicinity of resonance. X,, the equivalent reactance of the crystal plate, 
increases rapidly as the frequency departs frem resonance. R,and Xo, the 
equivalent resistance and reactance of the wire resonance, are typical of an 
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anti-resonant electrical network. The curve labeled (X; + X:) shows the 
effect of the wire resonance on the response of the crystal plate. It may be 
observed that the apparent resonance has been reduced by a small frequency 
decrement. The amcunt of frequency shift and the increase in effective 
res:stance depend on the Q of the wire resonance, its frequency location 
compared with the resonance of the crystal plate, the mass of the wire rela- 
tive to that of the quartz plate, and the distance frcm the node to the point 
at which the wire is actually fastened to the plate. 

The slope of the frequency-reactance characteristic corresponding to the 
mechanical resonance of the quartz plate is very steep and the effect of the 
wire resonance will be noticed only when an anti-resonant frequency of the 
wire is close to the resonant frequency of the plate. The changes in resonant 
frequency and effective resistance due to wire resonance have been measured 
for some filter crystal units and the measurements are tabulated in Table IT. 


TABLE II 
EFFECT OF WIRE VIBRATIONS ON THE RESISTANCE OF A QUARTZ CRYSTAL PLATE 


oss : “rysta Vistar f Wi Maxi imum 
Mode of Vibration | Resonant Crystal Distance of Wire Maximum | Maximun 


Crystal Type | Mforrysal  [Prequency| Masso | from Nodal Frequency incest i 

+5° X-Cut Flexural 12 kc 51 (N) .060” +2.0 250% 

+5° X-Cut Longitudinal 164 ke .142 (N) .012” +30 640° 
—18° X-Cut Longitudinal 335 ke .075 (M) .002” +90 300% 
—18° X-Cut Longitudinal 552 ke .068 (N) 0.0” +75 1100% 
5th Harmonic 


GT Longitudinal 164 ke .98 (N) .011” +12 370% 
(N) Specified Dimension 
(M) Measured Dimension 


t 


~- 


The relation between the length of a wire and the frequencies at which 
will resonate in flexural modes is expressed by the following equation: 


j=om wd 
8xf 


where v is the velocity of sound in the wire 
d is the diameter of the wire 
lis the length of the wire 
f is the frequency of wire resonance in cycles per second 
m is a number that cepends on the manner in which the ends of the 
wire can move. 

At a particular frequency and for wire of a particular material and diam- 
eter there is a series of critical wire lengths which must be avoided. The 
critical lengths are thcse which cause the wire to present a high impedance 
to the motion of the plate. This high impedance may be considered, from 
the electrical point of view, as corresponding toan anti-resonance of the wire. 


| 
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The critical lengths are defined by the series of numbers m = (n + } 

where » takes the values 1, 2, 3, etc. and apply to successive modes of 
bar clamped at both ends. Beyond the first mode, the critical wire lengt! 
are spaced at equal intervals, corresponding to increments of m each equ 
to mr. 

There is alsoa series of wire lengths which will present minimum impedanc: 
to the motion of the plate. These may be considered as corresponding to a 
resonance of the wire. These minima of impedance are obtained for lengths 
of wire defined by the series of numbers m = (n — 3) mw. They apply to a 
bar which is clamped at one end and, while free to vibrate at the other end 
is constrained to a slope perpendicular to the plate. 

In selecting a desirable length for the supporting wire, it is not essentia| 
that this length be such as to cause the wire to present minimum impedanc« 
to motion of the plate. As a matter of fact, since the wire is of relative], 
low characteristic impedance a small departure from the critical length is 
sufficient to avoid trouble from wire resonance. In order to allow for as 
wide a manufacturing tolerance as possible the supporting wire is usually 
designed to have a length half-way between two successive critical lengths. 
For a 6.3-mil phosphor-bronze wire, the spacing between successive critica] 
lengths ranges from about 58 mils at 100 kc to about 15 mils at 1000 kc. 
Hence, even at 100 kc the length of the supporting wire must be controlled 
within a tolerance of about 20 mils. 

These supporting wires are formed to have definite bends along their 
length and the location of these bends varies slightly from one wire to 
another. In addition, the wires are terminated by solder at both ends. 
Because of these complications it is impractical to meet such close tolerances 
on the effective length of the wires. Furthermore, a wire that does have a 
suitable effective length at room temperatures may exhibit sufficient change 
of properties with variations in temperature so that it becomes of critica! 
length at some other operating temperature. 

Much of the difficulty due to wire resonance is avoided by use of a solder 
ball on the supporting wire, as described in Chapters VIII and XIII. The 
solder ball is located near the quartz plate. Since it serves as a clamp at that 
point, it makes the supporting wire short. By locating and forming the 
solder ball accurately, the length of the supporting wire is controlled within 
a close tolerance. Further, since the wire is shortened by use of the ball 
it is less affected by changes .. temperature. Experience at about 500 kc 
indicates that a tolerance of about 10 mils in locating the solder ball is prac- 
tical and has provided satisfactory operation between —40 C and +85 C. 


14.4 NrrEp For CLEANLINESS AND Low RELATIVE HuMIDITY 


One of the most serious difficulties encountered in manufacturing quartz 
crystal plates is that of assuring sufficient cleanliness. Even minute par- 
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ticles of foreign matter will introduce appreciable changes in crystal per- 
formance. 

Usually, the presence of foreign matter will act to load the crystal and 
will reduce the resonant frequency but there are also instances where the 
added matter tends to stiffen the plate and increase its frequency. The 
latter has been observed to occur as the result of the deposit of a thin film 
of rosin on the surface of the plate. In the presence of foreign matter on 
the surface of the plates, the performance will be unstable with time and 
temperature even after the plate is sealed into a container. Also, erratic 
variations are observed as the plate is shifted from a normal atmosphere to 
a container which is evacuated or filled with dry air. Experience has shown 
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Fig. 14.11—Effect of humidity on the discrimination of channel crystal filters. To 
prevent decrease of discrimination with increase of relative humidity, the crystal units 
must be hermetically sealed. 


that elaborate precautions for insuring cleanliness are justified by the time 
saved in the adjusting processes. 

The need for cleanliness is closely related to the effect of humidity on the 
insulation resistance of crystal units. As used in filters, crystal units must 
provide extremely high impedances at their anti-resonant frequencies. 
These impedances may be as high as 100 megohms. With clean crystal 
plates in relatively dry atmospheres, such insulation resistance can be main- 
tained up to 1000 kc. However, even a trace of salts or other types of con- 
tamination will make the insulation resistance highly sensitive to moisture 
in the adjacent air. While it is relatively difficult to measure insulation 
resistance at high carrier frequencies, the effect of the reduced insulation 
due to moisture is evident on inspection of the discrimination characteristics 
of the filters. For example, Fig. 14.11 illustrates the result of high relative 
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humidity on the transmission characteristic of a typical crystal filter. I: 
may be observed that the discrimination almost disappears for a relativ: 
humidity of 80%. These measurements were mace on a filter co. 

taining well cleaned crystal plates. It will ke found frequently that a 
unsatisfactory discrimination characteristic is produced by considerab|; 
lower values of relative humidity when the plates are not so clean. Ex 
perience has shown that it is impractical to let the relative humidity sur 
rounding the crystal plate exceed 46% for satisfactory filter performance 
When a high degree of accuracy is required, the plates are assembled 
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Fig. 14.12.—Change of resonant frequency of GT-cut plates due to increase of trans 
mitted current. 


in a unit which is either evacuated or filled with air at a relative humidity 
of less than 5%. 


14.5  EFrFect oF CURRENT LEVEL 


Crystal units will undergo change in effective resistance and in frequency 
of resonance as the current transmitted is increased. Some change might 
be expected due to the heating of the plate by the dissipative loss associated 
with the transmission of current. However, the effects are not identical 
with those obtained with a change in ambient temperature. Appreciable 
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changes have been observed even when using GT-cut plates adjusted to zero 
temperature coefficient as shown, for example, in the curves of Fig. 14.12. 
Also, it has been observed that after a plate has been driven hard and the 
transmitted current then reduced, the original resonant frequency is re- 
stored only after a considerable time interval. The data of Fig. 14.12 pro- 
vides a rough correlation between stability and current levels. For ex- 
ample, if the stability desired for a crystal unit using a GT type plate be in 
the order of one part per million, the circuit design should be such as to 
keep the current level in the plate below about 10 microamperes per milli- 
meter of width. 

The parameter used in these paragraphs for measuring current levels is 
the current per unit of width. This appears to be useful as a common basis 
for comparing various plates of any one cut and mode of vibration. Theo- 
retically, in the case of a plate vibrating longitudinally, the current, 7, per 
unit of width, w, is directly proportional to the elongation per unit of length, 
yy, as shown by following equation: 


I/jw=Ky, 


where K is a constant which depends on the cut of plate and mode. 

Figure 14.12 also illustrates the importance of the surface condition of 
the plates. Curve A is the average frequency-current characteristic for a 
group of crystal units using plates etched for twenty minutes in 47% hydro- 
fluoric acid and curve B the average characteristic for a group of crystal 
units using plates etched for over forty minutes but less than ninety minutes. 
Evidently, crystal units using plates which have been etched for a long 
period exhibit a frequency-current characteristic which is appreciably more 
constant than those using plates etched for a shorter period. 





Mathematical Analysis of Random Noise 
By S. O. RICE 


INTRODUCTION 


HIS paper deals with the mathematical analysis of noise obtained by 

passing random noise through physical devices. The random noise 
considered is that which arises from shot effect in vacuum tubes or from 
thermal agitation of electrons in resistors. Our main interest is in the sta- 
tistical properties of such noise and we leave to one side many physica! 
results of which Nyquist’s law may be given as an example.’ 

About half of the work given here is believed to be new, the bulk of the 
new results appearing in Parts III and IV. In order to provide a suitable 
introduction to these results and also to bring out their relation to the work 
of others, this paper is written as an exposition of the subject indicated in 
the title. 

When a broad band of random noise is applied to some physical device, 
such as an electrical network, the statistical properties of the output are 
often of interest. For example, when the noise is due to shot effect, its 
mean and standard deviations are given by Campbell’s theorem (Part I) 
when the physical device is linear. Additional information of this sort 
is given by the (auto) correlation function which is a rough measure of the 
dependence of values of the output separated by a fixed time interval. 

The paper consists of four main parts. The first part is concerned with 
shot effect. The shot effect is important not only in its own right but 
also because it is a typical source of noise. The Fourier series representa- 
tion of a noise current, which is used extensively in the following parts, may 
be obtained from the relatively simple concepts inherent in the shot effect. 

The second part is devoted principally to the fundamental result that the 
power spectrum of a noise current is the Fourier transform oi its correlation 
function. This result is used again and again in Parts III and IV. 

A rather thorough discussion of the statistics of random noise currents 
is given in Part III. Probability distributions associated with the maxima 
of the current and the maxima of its envelope are developed. Formulas 
for the expected number of zeros and maxima per second are given, and a 
start is made towards obtaining the probability distribution of the zeros. 

When a noise voltage or a noise voltage plus a signal is applied to a non- 


1 An account of this field is given by E. B. Moullin, “Spontaneous Fluctuations of 
Voltage,” Oxford (1938). 
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linear device, such as a square-law or linear rectifier, the output will also 
contain noise. The methods which are available for computing the amount 


| of noise and its spectral distribution are discussed in Part IV. 
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SUMMARY OF RESULTS 


Before proceeding to the main body of the paper, we shall state some of 
the principal results. It is hoped that this summary will give the casual 
reader an over-all view of the material covered and at the same time guide 
the reader who is interested in obtaining some particular item of informa- 


| tion to those portions of the paper which may possibly contain it. 


Part I—Shot Effect 


Shot effect noise results from the superposition of a great number of 
disturbances which occur at random. A large class of noise generators 
produce noise in this way. 

Suppose that the arrival of an electron at the anode of the vacuum tube 
at time ¢ = 0 produces an effect F(¢) at some point in the output circuit. 
If the output circuit is such that the effects of the various electrons add 
linearly, the total effect at time ¢ due to all the electrons is 


+2 
I(t) = a, F(t — &) (1.21) 


where the k™ electron arrives at ¢, and the series is assumed to converge. 
Although the terminology suggests that J(¢) is a current, and it will be 
spoken of as a noise current, it may be any quantity expressible in the form 
(1.2-1). 

1. Campbell’s theorem: The average value of J(¢) is 


cal +2 
I(t) = » | F(t) dt (1.2-2) 
and the mean square value of the fluctuation about this average is 


ave. (I(t) — T(H)) = v ‘. F’(t) dt (1.2-3) 
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where v is the average number of electrons arriving per second at the an: 
In this expression the electrons are supposed to arrive independently an: 
random. ve ” dt is the probability that the length of the interval betwee: 
two successive arrivals lies between ¢ and ¢ + dt. 

2. Generalization of Campkell’s theorem. Campbell’s theorem gives 
information about the average value and the standard deviation of the 
probability distribution of /(t). A generalization of the theorem gives 
information about the third and higher order moments. Let 


+00 
I(t) = >> a F(t — &) (1.5-1 


where F(t) and ¢, are of the same nature as thcse in (1.2-1) and ---a 
d2,***@,*** are independent randcm variables all having the sam« 
distribution. Then the n™ semi-invariant of the probability density P(/ 
of J = I(t) is 


—_ +e 


1. = va" [F(t)]" de (1.5-2 


The semi-invariants are defined as the ccefficients in the expansion of the 
characteristic function f(x): 


An . n - 
loge f(u) = » =a (iu) (1.5-3 
where 
° +e . 
f(u) = ave, e = P(I)e'" dI 


The moments may be computed frem the )’s. 
3. As »-—> the probability density P(Z) of the shot effect current ap 
proaches a normal law. The way it is approached is given by 


4 
PIT) wa (x) - Me o (x) 
? -7 (1.6-3) 
Mo 2 ; 
+ [Ms oe (x ) i> Soe 5 KR e” «) | + a 
where the \’s are given by (1.5-2) and 
_ ee I — T (n) P a ; a” ines 
ee 


Since the \’s are of the order of v, ¢ is of the order of y* and the orders of 
respectively. A 


—1/2 2. 


-1 a4 =f a ~3/2 3/2 
a, X\go*, Xo and Ajo” are v oe aol and vy 
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possible use of this result is to determine whether a noise due to random in- 
dependent events occuring at the rate of » per second may be regarded as 
“random noise”’ in the sense of this work. 

4. When /(t), as given by (1.5-1), is analyzed as a Fourier series over an 
interval of length T a set of Fourier coefficients is obtained. By taking 
many different intervals, all of length 7, many sets of coefficients are 
obtained. If » is sufficiently large these coefficients tend to be distributed 
normally and independently. A discussion of this is given in section 1.7 


Part II—Power Spectra and Correlation Functions 


1. Suppose we have a curve, such as an oscillogram of a noise current, 
which extends from? = 0 tot = «. Let this curve be cenoted by /(t). 
The correlation function of J(t) is ¥(r) which is defined as 


. 
Y(r) = Limit - | T(t + 7) dt (2.1-4) 
to T Jp 
where the limit is assumed to exist. This function is closely connected 
with another function, the power spectrum, w(f), of Z(¢). J(t) may be 
regarded as composed of many sinusoidal components. If J(t) were a 
noise current and if it were to flow through a resistance of one ohm the 
average power dissipated by those components whose frequencies lie be- 
tween f and f + df would be w(/) df. 
The relation between w(f) and ¥(r) i 





w(f) = 1[ (7) cos 2nfr dr (2.1-5) 
0 
(7) = [ w(f) cos 2xfr df (2.1-6) 
0 
When /(¢) has no d.c. or periodic components, 
l 
w(f) = Lim pals (2.1-3) 
Tx 7 
where 
S(/) = [ Iie" at. 
The correlation function for 
I(t) = A + C cos (2rfet — ¢) 
is 
V(r) = A’ +5 — cos 2rfot (2.2-3) 





These results are discussed in sections 2.1 to 2.4 inclusive 
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2. So far we have supposed J(¢) to be some definite function for which a 
curve may be drawn. Now consider J(t) to be given by a mathematica! 
expression into which, besides /, a number of parameters enter. w(f) and 
y(r) are now obtained by averaging the integrals over the possible values 
of the parameters. This is discussed in section 2.5. 

3. The correlation function for the shot effect current of (1.2-1) is 


+00 +00 2 
Wilh sae [ F()F(t + 1) dt + |» F(t) a| (2.6-2) 


The distributed portion of the power spectrum is 


wi(f) = 2v| s(f) |? 


where 
+00 


s(f) = F(t)e""" dt (2.6-5) 


The complete power spectrum has in addition to w;:(f) an impulse func- 
tion representing the d.c. component J(¢). 

In the formulas above for the shot effect it was assumed that the expected 
number, v, of electrons per second did not vary with time. A case in which 
v does vary with time is briefly discussed near the end of Section 2.6. 

4, Random telegraph signal. Let J(¢) be equal to either a or —a so that 
it is of the form of a flat top wave, and let the lengths of the tops and bot- 
toms be distributed independently and exponentially. The correlation 
function and power spectrum of J are 


Vr) = ae (2.7-4) 
2 2 
w(f) = 74 (2.7-5) 


where u is the expected number of changes of sign per second. 

Another type of random telegraph signal may be formed as follows: Divide 
the time scale into intervals of equal length 4. In an interval selected at 
random the value of J(¢) is independent of the value in the other intervals 
and is equally likely to be +a or —a. The correlation function of J(é) is 


zero for |7| > h and is 
ate . 083 
a (1 i ) 


for 0 < | r| < hand the power spectrum is 


w(f) = 2h (c=) (2.7-9) 
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5. There are two representations of a random noise current which are 
especially useful. The first one is 


N 
I(t) = >> (an CoS wal + 5, sin wn!) (2.8-1) 
n=l 


where a, and 5, are independent random variables which are distributed 
norma!ly about zero with the standard deviation Vw ff and where 
W, = 2rf,, fn = nf 


The second one is 
N 
I(t) = pa C, COS (wnt _ Gn) (2.8-6) 
n=l 


where ¢, is a random phase angle distributed uniformly over the range 
(0, 2m) and 


Cn = [2w(fn)df}'” 


At an appropriate point in the analysis NV and Af are made to approach 
infinity and zero, respectively, in such a manner that the entire frequency 
band is covered by the summations (which then become integrations). 

6. The normal distribution in several variables and the central limit 
theorem are discussed in sections 2.9 and 2.10. 


Part III—Statistical Properties of Noise Current 


1. The noise current is distributed normally. This has already been 
discussed in section 1.6 for the shot-effect. It is discussed again in section 
3.1 using the concepts introduced in Part II, and the assumption, used 
throughout Part ITI, that the average value of the noise current /(f) is zero. 
The probability that J(¢) lies between I and J + d/ is 

aT —12/2h9 
me (3.1-3) 


where yo is the value of the correlation function, ¥(7), of 7(t) at r = 0 

bo = 40) = | winds (3.1-2) 
w(f) being the power spectrum of /J(¢). yo is the mean square value of 
I(t), ice., the r.m.s. value of I(t) is po. 


The characteristic function (ch. f.) of this distribution is 


ave, e"? = exp — ar (3.1-6) 
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2. The probability that 7(¢) lies between 7; and 7, + dJ, and I(t + 
lies between J. and J> + dJ2 when ¢ is chosen at random is 


f,2 — y’] 0 dT aF exp | oll — Melt 4 Met (3.9 
on 2(¥o0 — ve) | 
where y, is the correlation function ¥(7) of Z(f): 


/ 


Y(7r) = [ wh cos 2nfr df (3.2-3 


The ch. f. for this distribution is 


ave, OTHE o exp | —¥ at + v) — vo (3.2-7 


< 


The expected number of zeros per second of J(f) is 


1 Eve Are ful) a] 3.3411 
T ¥(0) if nay | 3. 
F w( 


assuming convergence of the integrals. The primes denote differentiation 
with respect to 7: 





For an ideal band-pass filter whose pass band extends from f, to f, the ex 
pected number of zeros per second is 


1fp— py? ps 
2| se = fl (3.3-12) 


When f, is zero this becomes 1.155 f, and when f, is very nearly equal to 
fo it approaches fy + fa. 

4. The problem of determining the distribution function for the length 
of the interval between two successive zeros of I(t) seems to be quite diff- 
cult. In section 3.4 some related results are given which lead, in scme 
circumstances, to approximations to the distribution. For example, for 
an ideal narrow band-pass filter the probability that the distance between 
two successive zeros lies between 7 and 7 + dr is approximately 


dr a 


+e — a 
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+ 
om aer + fa)” = 1 
- oF Va [ey a ’ 1 = Pyeey a 
, and f, being the upper and lower cut-off frequencies 
5. In section 3.5 several multiple integrals which éccur in the work cf 
Part III are discussed. 
2-3 6. The distribution of the maxima of J/?) is discussed in section 3.6. The 


expected number of maxima fer second is 


(4) 71/2 Ww ) df 

ey [ fone a 

2-7 |— or) =|s— (3.6-6) 
alls [ rong 

0 


cr 1/2 





For a band-pass filter the expected number of maxima per second is 


5 51/2 
; fo Pee 4 (36-7 
= ———— ‘ (3.0-7) 
= = 3 x « ‘ 
, 1] Pe) he = Ja 
For a low-pass filter where fg = 0 this number is 0.775 f, 
: The expected number of maxima per second lying above the line /(/) = /, 
ition 


is approximately, when /, is large, i 


—s2/s 
1/2 


e '*° & 3{[the expected number of zeros of I persecond]  (3.6-11 


where Yo is the mean square value of /(f). 

For a low-pass filter the probability that a maximum chosen at random 
pe ex from the universe of maxima lies between J and J + d/ is approximately, 
when / is large, 


5 —y2/2 dl é 
12) = 7 plz 3.6-9) 
. 3 vo 
where 
il to 
I 
7 = yli2 
ngth Yo 
diffi- ~ ; aes se ’ a, rary | en 
7. When we pass noise through a relatively narrow band-pass filter one 
¢ > . - . . 
a of the most noticeable features of an cscillogram of the output current is 
or ° . ’ . >» >. . . . 
?, its fluctuating envelope. In sections 3.7 and 3.8 some statistical properties 
ween 


of this envelope, denoted by R or R(t), are derived. 
The probability that the envelope lies between Rand R + dRis 


R 


—¢ "0 gR (3.7-10) | 
Yo 


J 
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where yo is the mean square value of /(/). The probability that R(t) lie: 
between R, and R; + dR, and at the same time R(t + 7) lies betwee: 
R; and Rz + dR, when ¢ is chosen at random is obtained by multiplying 
(3.7-13) by dR, dR;. For an ideal band-pass filter, the expected number 
of maxima of the envelope in one second is 


64110(f, — fa) (3.8-15) 
When R is large, say y > 2.5 where 
R 1/2 
yu yi? Yo =r.m.s, value of (2), 
0 


the probability that a maximum of the envelope, selected at random from 
the universe of such maxima, lies between R and R + @R is approximately 


1.13(y? — 1)e"¥"? a 
vo 
A curve for the corresponding probability density is shown for the range 
0<y<4. Curves which compare the distribution function of the maxima 
of R with other distribution functions of the same type are also given. 
8. In section 3.9 some information is given regarding the statistical 
behavior of the random variable: 


ti+7 

E= / I’(t) dt (3.9-1) 
1 

where ¢, is chosen at random and J(¢) is a noise current with the power 

spectrum w(f) and the correlation function ¥(7). The average value 

mr of E is Tyo and its standard deviation o7 is given by (3.9-9). For a 

relatively narrow band-pass filter 


eh: an 

mr WT(fo — fa) 
when T(fp — fa) > 1. This follows from equation (3.9-10). An ex- 
pression which is believed to approximate the distribution of E is given by 
(3.9-20). 

9. In section 3.10 the distribution of a noise current plus one or more 

sinusoidal] currents is discussed. For example, if J consists of two sine waves 
plus noise: 


I = Pcos pt+ Qcos q+ In, (3.10-20) 
where p and q are incommensurable and the r.m.s. value of the noise cur- 


rent Iy is yo”, the probability density of the envelope R is 


R rf rJo(Rr) Jo Pr)Jo(Or)e*""*? dr (3.10-21) 


where J( ) is a Bessel function. 
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Curves showing the probability density and distribution function of R, 
when Q = 0, for various ratios of P/r.m.s. Iy are given. 

10. In section 3.11 it is pointed out that the representations (2.8~1) 

and (2.8-6) of the noise current as the sum of a great number of sinusoidal 

/ components are not the only ones which may be used in deriving the results 

given in the preceding sections of Part III. The shot effect representation 


+2 
I(t)'= Qo F(t — ti) 
studied in Part I may also be used. 


Part IV—Noise Through Non-Linear Devices 


1. Suppose that the power spectrum of the voltage V applied to the 


square-law device 
I = aV’ (4.1-1) 


is confined to a relatively narrow band. The total low-frequency output 
current J,¢ may be expressed as the sum 


Tet = Tae + Its (4.1-2) 


where Ig. is the d.c. component and J¢; is the variable component. When 
none of the low-frequency band is eliminated (by audio frequency filters) 


2 
le = “ (4.1-6) 


where Ris the envelope of V. If V is of the form 
V = Vy + P cos pt + Q cos gt, (4.14) 


where Vy is a noise voltage whose mean square value is Yo, then 


32 2 
Tae = «(vo ++ ¢) 


x § 2 2 2 2 PO? 
Ig=a)}pYotPyptQyort — (4.1-16) 
2. If instead of a square-law device we have a linear rectifier, 
fo V <0 
I= lav, V>0 (4.2-1) 


the total low-frequency output is 


T 
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When V is a sine wave plus noise, Vy + P cos pt, 


1/2 
7 «(¥) Fi= 


. 2 
ie =“ (P? + dW) (4.2-¢ 
= 


tom 


31; —x) (4.2-3 


where ,/*; is a hypergeometric function and 


P’ _ Ave. sine wave power 
= 7 _ : (4.2-4 
2YWo Ave. noise power 
When x is large 

2 

2 a Yo 1 ‘ 

ly ~ S*}1 - = +. (4.2 
™ 4x 


If V consists of two sine waves plus noise, J. consists of a hypergecmet: 
function of two variables. The equations running frem (4.2-9) to (4.2-1 
are concerned with this case. About the only simple equation is 


A) 


2 a yy 2 ‘ 
lie= —(WWwo +P +0) (4.2-14 
r: 


3. The expressions (4.1-6) and (4.2-2) for 7,¢ in terms of the enveloy: 
R of V, namely 


a R ‘ aR 


nd -—--, 
TT 
are special cases of a more general result 
1 -"~ 
T.¢ = Ao(R) = | F(iu)Jo(uR) du. (4.3-11 
Tc 


In this expression Jo(uR) is a Bessel function. The path of integration ( 
and the function F(iz) are chosen so that the relation between J and V ma\ 
be expressed as 


[= =f F(iu)e’’” du. (4A-1 
2rJc 


A table giving F(iu) and C for a number of common non-linear devices is 


shown in Appendix 4A. 
If this relation is used to study the biased linear rectifier. 


mm 0, V<B 
~ \V-—-B, V>B 
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for the case in which V is Vy + P ces Pt, we find 


By? B+ vs 


Tae ~ "9 be OR 2xP 
ah cate hed (4.3-17) 
eg Se al 

f 7 P? 


when P >| BI, P* > Wo where yo is the mean square value of Vy. 

4. When V is confined to a relatively narrow band and there are no 
audio-frequency filters, the probability density and all the associated sta- 
tistical properties of J,¢ may be obtained by expressing /,¢ as a function 
of the envelope R of V and then using the probability density of R. When 
V is Vw + P ces pt + Q cos et this probability density is given by the in- 
tegral, (3.10-21) (which is the integral containing three Bessel functions 
stated in the above summary of Part III). When V consists of three sine 
waves plus noise there are four Jo’s in the integrand, and so on. Expres- 
sions for R" when R has the above distribution are given by equations 
(3.10-25) and (3.10—27). 

When audio-frequency filters remove part of the low-frequency band the 
statistical preperties, except the mean square value, of the resulting cur- 
rent are hard tocompute. Insection 4.3 it isshown that as the output band 
is chosen narrower and narrower, the statistical properties of the output 
current approach those of a random noise current. 

5. The sections in Part 1V from 4.4 onward are concerned ‘with the 
problem: Given a non-linear device and an input voltage consisting of noise 
alone or of a signal plus noise. What is the power spectrum of the output? 
A survey of the methods available for the solution of this problem is given 
in section 4.4, 

6. When a noise voltage Vy with the power spectrum w(f) is applied to 
the square-law device 


I = av? (4.1-1) 
the power spectrum of the output current J is, when f # 0, 
9 ices 
W(f) = a | w(x)w(f — x) dx (4.5-5) 
x 


where w(—<x) is defined to equal w(x). The power spectrum of J when V 
is either P cos pt + Vy or 


O(1 + k cos pt) cos gt + Vw 


is considered in the portion of section 4.5 containing equations (4.5-16) to 


(4.5-17). 
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7. A method discovered independently by Van Vleck and North shows 
that the correlation function ¥(7) of the output current for an unbiased 
linear rectifier is 


Wr) = Ye + aR, [ -) 354; “| (4.7-6) 
where the input voltage is Vy. The correlation function ¥(r) of Vy is 
denoted by y, and the mean square value of Vy iso. The power spectrum 
W(f) of IJ may be obtained from 


W(f) = 4 | V(r) cos 2xfr dr (4.6-1) 

‘) 
by expanding the hypergeometric function and integrating termwise using 
G.(f) = I W; Cos 2mnfr dr. (4C-1) 


Appendix 4C is devoted to the problem of evaluating the integral for G,(f). 

8. Another method of obtaining the correlation function ¥(7) of J, termed 
the “characteristic function method,” is explained in section 4.8. It is 
illustrated in section 4.9 where formulas for V(r) and W(f) are developed 
when the voltage P cos pt + Vy is applied to a general non-linear device. 

9. Several miscellaneous results are given in section 4.10. The char- 
acteristic function method is used to obtain the correlation function for a 
square-law device. The general formulas of section 4.9 are applied to the 
case of a v™ law rectifier when the input noise spectrum has a normal law 
distribution. Some remarks are also made concerning the audio-frequency 
output of a linear rectifier when the input voltage V is 


Q(1 + rcos pt) cos gt + Vy. 


10. A discussion of the hypergeometric function :F (a; ¢; x), which often 
occurs in problems concerning a sine wave plus noise, is given in 
Appendix 4B. 


PART I 


THE SHOT EFFECT 


The shot effect in vacuum tubes is a typical example of noise. It is due 
to fluctuations in the intensity of the stream of electrons flowing from the 
cathode to the anode. Here we analyze a simplified form of the shot effect. 
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1.1 THe PROBABILITY OF ExActL_y K ELECTRONS ARRIVING AT THE 
ANODE IN TimE T 


The fluctuations in the electron stream are supposed to be random. We 
shall treat this randomness as follows. We count the number of electrons 
flowing in a long interval of time T measured in seconds. Suppose there 
are Ky. Repeating this counting process for many intervals all of length 
T gives a set of numbers K2, Kz --- Ky where M is the total number of 
intervals. The average number », of electrons per second is defined as 


» on kin Sit Bet --- + Ke (1.1-1) 





where we assume that this limit exists. As M is increased with T being 
held fixed some of the K’s will have the same value. In fact, as M increases 
the number of K’s having any particular value will tend to increase. This 
of course is based on the assumption that the electron stream is a steady 
flow upon which random fluctuations are superposed. The probability of 
getting K electrons in a given trial is defined as 
4K) = Lim Number of trials giving exactly K electrons (1.1-2) 
M—2 M 
Of course p(K) also depends upon T. We assume that the random- 
ness of the electron stream is such than the probability that an electron 
will arrive at the anode in the interval (if, 4 + Af) is vAt where Af is 
such that vAt < 1, and that this probability is independent of what has 
happened before time ¢ or will happen after time ¢ + At. 
This assumption is sufficient to determine the expression for p(K) which is 





?(K) = 7)" ” Siig (1.1-3) 


This is the “law of small probabilities” given by Poisson. One method 
of derivation sometimes used can be readily illustrated for the case K = 0. 
Thus, divide the interval, (0, 7) into M intervals each of length At = = 
At is taken so small that vAt is much less than unity. (This is the “small 
probability” that an electron will arrive in the interval At). The prob- 
ability that an electron will not arrive in the first sub-interval is (1 — vAfZ). 
The probability that one will not arrive in either the first or the second 
sub-interval is (1 — vAt)*. The probability that an electron will not arrive 
in any of the M intervals is (1 — vAt)™. Replacing M by 7/At and letting 
At — 0 gives 


p(0) = «”* 
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The expressions for p(1), p(2), --- p(K) may be derived in a somewh 
similar fashion. 


1.2 STATEMENT OF CAMPBELL’S THEOREM 
Suppose that the arrival of an electron at the anode at time ¢ = 0 produces 
an effect F(t) at some point in the output circuit. If the output circuit 
is such that the effects of the various electrons add linearly, the total effect 
at time ¢ due to all the electrons is 
+0 
I(t) = >> F(t — &) (1.2-1 
k=—oo 
where the &‘” electron arrives at ¢, and the series is assumed to converre. 
Campbell’s theorem’ states that the average value of I(t) is 


I(t) =v» / F(t) dt (1.2-2 
—<e 
and the mean square value of the fluctuation about this average is 
+20 
(I(t) — I(i))° = » | F’(t) dt (1.2-3 
— GO 


where v is the average number of electrons arriving per second. 

The statement of the theorem is not precise until we define what we mean 
by “average”. From the form of the equations the reader might be tempted 
to think of a time average; e.g. the value 


ee 
Lim | I(t) dt (1.2-4) 
T Jo 


y ed -) 
However, in the proof of the theorem the average is generally taken over 
a great many intervals of length T with ¢ held constant. The process is 
somewhat similar to that employed in (1.1) and in order to make it clear 
we take the case of /(t) for illustration. We observe /(t) fcr many, say M, 
intervals each of length T where T is large in comparison with the interval! 
over which the effect /(¢) of the arrival of a single electron is appreciable 
Let ,J(t’) be the value of J(¢), t’ seconds after the beginning of the n™ in 
terval. ¢’ is equal to?¢ plus a constant depending upon the beginning time 
of the interval. We put the subscript in front because we wish to reserve 
the usual place for another subscript later on. The value of J(¢’) is then 
defined as 
I(t’) = Limit (Z(t) + l(t’) +--+ + ult’) (3.2-5 
Mo M 
and this limit is assumed to exist. The mean square value of the fluctua- 
tion of J(t’) is defined in much the same way. 


2 Proc. Camb. Phil. Soc. 15 (1909), 117-136, 310-328. Our proof is similar to one give! 
by J. M. Whittaker, Proc. Camb. Phil. Soc. 33 (1937), 451-458. 
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Actually, as the equations (1.2-2) and (1.2-3) of Campbell’s theorem 
show, these averages and all the similar averages encountered later turn 
cut to ke incepencent of the time. When this is true and when the M in- 
tervals in (1.2-5) are taken consecutively the time average (1.2-4) and the 
average (1.2-5) keccme the same. To show this we multiply both sides of 
1.2-5) by dt’ and integrate from 0 to T: 


M 


er ' 
’) = Limit ml (t’) dt’ 
10) = Limit orp 20 I , 

(1.2-6) 


1 M7 
= Limit wd I(t) dt 
m--o MT Jo 
and this is the same as the time average (1.2-4) if the latter limit exists 


1.3. PRooF OF CAMPBELL’S THEOREM 

Consider the case in which exactly A electrons arrive at the anode in an 
interval of length 7. Eefore the interval starts, we think of these K elec- 
trons as fated to arrive in the interval (0, 7) but any particular electron is 
just as likely to arrive at one time as any other time. We shall number 
these fated electrons frcm one to K for purposes of identification but it is to 
be emphasized that the numbering has nothing to do with the order of ar- 
rival. Thus, if ¢, be the time of arrival of electron number &, the probability 
that ¢, lies in the interval (¢,¢ + dt) is dt/T. 

We take T to be very large compared with the range of values of ¢ for 
which F(¢) is appreciably different from zero. In physical applications 
such a range usually exists and we shall call it A even though it is not very 
definite. Then, when exactly A electrons arrive in the interval (0, 7) the 
effect is approximately 


‘ 
I(t) = >> F(t — &) (1.3-1) 
k=1 
the degree of approximation being very good over all of the interval except 
within A of the end points. 

Suppose we examine a large number M of intervals of length T. The 
number having exactly K arrivals will be, to a first approximation M p(K 
where p(K) is given by (1.1-3). Fora fixed value of ¢ and for each interval] 
having K arrivals, Jx(¢) will have a definite value. As M— ~, the average 
value of the Jx(t)’s, obtained by averaging over the intervals, is 


T it T ‘ K 
rah = [tt [MEF ru 9 
0 7 /0 7 k=l 


K T dt 
= —~ F(t — t) 
pa 7 . 


(1.3-2) 
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and ifA <t < T — A, we have effectively 
> +00 
Tah = ; [ F(t) dt (1.3-3) 


If we now average J(t) over all of the M intervals instead of only over 
those having K arrivals, we get,as M— ~, 


I(t) 


a 


>» p(K)I x(t) 


+00 
«eer al | F(t) dt 


ko TK! 
+o 
=v F(t) dt (1.3-4) 


and this proves the first part of the theorem. We have used this rather 
elaborate proof to prove the relatively simple (1.34) in order to illustrate a 
method which may be used to prove more complicated results. Of course, 
(1.3-4) could be established by noting that the integral is the average value 
of the effect produced by one arrival, the average being taken over one 
second, and that v is the average number of arrivals per second. 

In order to prove the second part, (1.2-3) of Campbell’s theorem we first 
compute /?(#) and use 


Um — 10) 


Pi) — 2 1()I(t) + TO 
[?(8) - I(t)’ (1.3-5) 
From the definition (1.3-1) of Zx(é), 


ll 


Txt) = DDD F(t — &)F(t — tn) 


k=1 m=1 


Averaging this over all values of ¢; , fg, «++ ¢x with ¢ held fixed as in (1.3-2 


4 
RO = ee dt f dix Fy — 4) F(t — te) 
> m= 1 a & 


The multiple integral has two different values. If k = mits value is 


[re- n) 


and if k * m its value is 


[ra-w%f F(t = tm) F 
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Counting up the number of terms in the double sum shows that there are K 


them having the first value and K* — K having the second value. Hence, 
ifA4<t< T — Awehave 


I(t) = . a F°(0) dt + eS If. F(t) a] 


Averaging over all the intervals instead of only those having K arrivals 
gives 


> o(K) 1(0) 


K=0 


P(t) 


II 


y [ F’(t) dt + I(t)? 


where the summation with respect to A is performed as in (1.3-4), and after 
summation the value (1.3-4) for /(¢) is used. Comparison with (1.3-5) 
establishes the second part of Campbell’s theorem. 


1.4 Tue DiIstRIBuTION oF /(2) 


When certain conditions are satisfied the proportion of time which J(é) 
spends in the range J, J + dJ is P(J)dI where, as v — ~, the probability 
density P(I) approaches 

= —(I—1)?/202 

wate e€ I (1.4-1) 
where J is the average of J(¢) given by (1.2-2) and the square of the standard 
deviation o7, i.e. the variance of J(¢), is given by (1.2-3). This normal 
distribution is the one which would be expected by virtue of the “‘central 
limit theorem” in probability. This states that, under suitable conditions, 
the distribution of the sum of a large number of random variables tends 
toward a normal distribution whose variance is the sum of the variances 
of the individual variables. Similarly the average of the normal distribu- 
tion is the sum of the averages of the individual variables. 

So far, we have been speaking of the limiting form of the probability 
density P(J). It is possible to write down an explicit expression for P(/), 
which, however, is quite involved. From this expression the limiting form 
may be obtained. We now obtain this expression. In line with the dis- 
cussion given of Campbell’s theorem, we seek the probability density P(/) 
of the values of J(¢) observed at ¢ seconds from the beginning of each of a 
large number, M, of intervals, each of length T. 
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Probability that J(¢) lies in range (J, J + d/) 


io) 


= }° (Probability of exactly A arrivals) x 
Kx 


(Probability that if there are exactly 

K arrivals, x(t) lies in (7, I + d/) 
Denoting the last probability in the summation by Px(J)d/, using notatio: 
introduced earlier, and cancelling out the factor d/ gives 


P(I) = >> p(K)Px(1) (1.4-2 
K=( 
We shall compute Px(/) by the method of “characteristic functions’” from 
the definition 


K 
Ix(t) = >> F(t — t) (1.3-1 
k=1 


of Ix(t). The method will be used in its simplest form: the probability that 
the sum 


X1 + a2 + ees + ae 


of K independent random variables lies between X and X + dX is 
1 in a u 
dX - xf. ” ios I (average value of e *") du (1.4-3 


The average value of e’ “", j.e., the characteristic function of the distribution 
of x, , is obtained by averaging over the values of «,. Although this is the 
simplest form of the method it is also the least general in that the integral 
does not converge for some important cases. The distribution which gives 
a probability of } that x = —1land 3 that a, = +1 is an example of such a 
case. However, we may still use (1.4-3) formally in such cases by employ- 
ing the relation 


+o : 
[ e'’“ du = 2n6(a) (1.4-4) 
wo 

where 6(a) is zero except at a = O where it is infinite and its integral from 
a= —etoa = +e is unity where e > 0. 


When we identify a, with F(¢ — ¢,) we see that the average value cf 


tZeu . 


1S 


T 
x | exp luk — 4) dt 
0 


3 The essentials of this method are due to Laplace. A few remarks on its history are 
given by E. C. Molina, Bull. Amer. Math. Soc., 36 (1930), pp. 369-392. An account o 
the method may be found in any one of several texts on probability theory. We men.ion 
“Random Variables and Probability Distributions,” by H. Cramér, Camb. Tract in 
Math. and Math. Phys. No. 36 (1937), Chap. [V. Also “Introduction to Mathematica 
Probability,” by J. V. Uspensky, McGraw-Hill (1937), pages 240, 264, and 271-278. 








in 


ps 









4-3 


ution 
is the 
tegral 
gives 
uch a 
iploy- 


4-4) 


from 


lue of 


ry are 
unt of 
en ion 
act in 
natica 


MATHEMATICAL ANALYSIS OF RANDOM NOISE 


4 


\il of the K characteristic functions are the same and hence, from (1.4-3), 


P,(I)dT is 
1 oi,» 1? K 
dI — / ” eae Gi exp linF(t — dr) du 
2r x 7 0 


Although in deriving this relation we have taken K > 0, it also holds for 
K = 0 (provided we use (1.4~4)). In this case Po’) = 6(J), because J = 0 
when no electrons arrive 

Inserting our expression for Px(J) and the expression (1.1-3) for p(K 
in (1.4-2) and performing the summation gives 


PUI) = : exp (itu — oT 
LT J-x 


T 
+} | exp [iuF'(t — dr) du = (1.4-5 
0 


The first exponential may be simplified scmewhat. Using 


. 


T 
yl = » | dr 
0 


permits us to write 


aT 


T 
—vT + » | exp liuF(t — r)| dr = » | (exp [tu (t — r)] — 1) dr 
0 0 


Suppose that A < ¢ < T — A where A is the range discussed in connection 
with equation (1.3-1). Taking | F(¢ — r)| = 0 for|¢ — +r] > A then 
enables us to write the last expression as 


+0 ' 
»| [e"? — 1) de (1.4-6 


oe 


Placing this in (1.4—5) yields the required expression for P(/): 


PUD) = ~ | exp (-it + » | [ers — 1} dt) du (1.4-7) 


An idea of the conditions under which the normal law (1.4-1) is ap- 
proached may be obtained from (1.4-7) by expanding (1.4-6) in powers of 


* 6 . . 3 . “ 
u and determining when the terms involving u° and higher powers of 1 





may be neglected. This is taken up for a slightly more general form of 
current in section 1.6. 
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1.5 EXTENSION OF CAMPBELL’S THEOREM 
In section 1.2 we have stated Campbell’s theorem. Here we shall give 
an extension of it. In place of the expression (1.21) for the J(¢) of the s! 
effect we shall deal with the current 


+00 


I(t) = >> a F(t — tk) (1.5-1 
k==—o 
where F(t) is the same sort of function as before and where --+ a1, de, 
a, , **+ are independent random variables all having the same distributio: 


It is assumed that all of the moments a” exist, and that the events occur at 
random 

The extension states that the mth semi-invariant of the probability density 
P(I) of I, where IJ is given by (1.5~1), is 


+00 
—— var [ [F(e)]" dt (1.5-2) 


where vy is the expected number of events per second. The semi-invariants 
of a distribution are defined as the coefficients in the expansion 


N 
i Na os N — 
loge (ave. e™“) = D> — (iu)” + o(u™) (1.5-3) 
n=] n!\ 
i.e. as the coefficients in the expansion of the logarithm of the characteristi: 
function. The X’s are related to the moments of the distribution. Thus if 
m,, m2, +++ denote the first, second --- moments about zero we have 


N 
avee“ =1+ >> m (iu)” + o(u™) 
n=l : 


By combining this relation with the one defining the \’s it may be shown that 


I =m, = Ai 
I? = m2 = Ao + Arm 
I? = mz = X3 + 2d + Ame 


It follows that A, = I and Xx = ave. (J — I)’. Hence (1.5~-2) yields the 
original statement of Campbell’s theorem when we set » equal to one and 
two and also take all the a’s to be unity. 

The extension follows almost at once from the generalization of expression 
(1.4-7) for the probability density P(J). By proceeding as in section 1.4 
and identifying a, with a,/(¢ — t) we see that 


aed 


T 
ave.e " = 7 g(a) da | exp [iuaF(t — t,)] dt, 
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where g(a) is the probability density function for the a’s. It turns out that 
the probability density P(J) of J as defined by (1.5-1) is 


1 [- ’ [- 
i sone c _ ( 
P(I) ae exp ( ilu + y g(a) da 


[ [ewer 1} dt) du (1.5-4) 
The logarithm of the characteristic function of P(/) is, from (1.5—4), 


+00 $20 
»[ q(a) da [ [ever — 1) dt 


> (iu)” a n es ~n 
= oo # q(a) daa F'(t) dt 


fut, n!\ re) 


Comparison with the series (1.5-3) defining the semi-invariants gives the 
extension of Campbell’s theorem stated by (1.5-2). 

Other extensions of Campbell’s theorem may be made. For example, 
suppose in the expression (1.5-1) for /(é) that 4,, te, +--+ tx, +++ while still 
random variables, are no longer necessarily distributed accerding to the 
laws assumed above. Suppose now that the probability density p(x) is 
given where x is the interval between two successive events: 


le i+ nm (1.5-5) 


lbotm=h + at x 


ll 


ts 
and soon. For the case treated above 
p(x) = ve”. (1.5-6) 


We assume that the expected number of events per second is still v. 
Also we take the special, but important, case for which 


F(t) = 0, t<0 (1.5-7) 
F(t) =e, t> 0. 


For a very long interval extending from ¢ = 4; tot = T + t, inside of which 
there are exactly K events we have, if ¢ is not near the ends of the interval, 


T(t) 


Il 


a\F(t — th) + aF(t —ti— mm) +-°- 
+ AxsiF (t — ty —_ x1 ccc a= Xx) 


ay,F(t’) + aoF(t! — x1) +++ + axyiF(! — m1 — +++ — xe) 
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Pt) = ai’) + a! — mm) + ++ + aka’ — m+) — 2x) 
+ 2aiaeF (t') Fit’ = x1) + eee Zaid cyl (t') F(t’ = Zeros = 


+ 2ara3F(t' — x)F(t) — xy — x2) fees + 


where t’ =¢—%,. If we integrate J°(t) over the entire interval 0 <t’ < 7 


and drop the primes we get approximately 
rT 9 5 9 
| Pod = Git + ade) 

0 


+ 2ara2ge(xi) + 2aiaxp(a1 + x2) + +++ + 2aaxsie(ar + ++: + 


+ 2araxe (x2) a cee a cee + 2aKaKkiy VK 


where 
+00 
g(2) = / F(i)F(t — x) dx 


When we divide both sides by T and consider K and T to be very large, 
K aj vee + Pov 
T a aaa ¢(0) = va?y(0) 


” 


i, P 
pias ae ¢(x1) + a243 9(%2) + +++ + ax @x41¢(Xx)| = T average Qy Arai ol x; 


= va? | ¢(x) p(x) dx 


1 K-1 

T [ai as y(m + x2) + > -| = T aVe. Ay Qrag o(XK + Xk41) 
; = wt [dn [dee p(m)p(ade(m + ) 
. 0 0 


Consequently 


if «a 
I(t) = Lim 4 | I°(t) di 
T Jo 


T—2 


ll 


va?e(0) + 2a" | [ p(x)e(x) dx 
f 


+ [ ax [ dx2 p(x) p(x) ¢(1 + x2) + <> 
0 0 


For our special exponential form (1.5-7) for F(), 


—ar 


¢(x) = 


2a 
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d the multiple integrals occurring in the expression for /?(¢) may be written 
terms of powers of 


q -| p(x)e"™ dx (1.5-8) 


Thus 


and since 
wand i FU) dt = va/a 


we have 


Equations (1.5-8) and (1.5-9) give us an extension of Campbell’s theorem 
subject to the restrictions discussed in connection with equations (1.5-5) 
and (1.5-7). Other generalizations have been made‘ but we shall leave the 
subject here. The reader may find it interesting to verify that (1.5-9) 
gives the correct answer when /(x) is given by (1.5-6), and also to investi- 
gate the case when the events are spaced equally. 


1.6 APPROACH OF DISTRIBUTION OF J TO A NORMAL LAw 


In section 1.5 we saw that the probability density P(7) of the noise current 
I may be expressed formally as 


+20 oo 
P(I) + =~ | exp| ith + >> (iu)"d, | du (1.6-1) 


n=1 


where A, is the #th semi-invariant given by (1.5-2). By setting 


9 
< 


A =o 


pati h i=! nae 


Go o 


4See E. N. Rowland, Proc. Camb. Phil. Soc. 32 (1936), 580-597. He extends the 
theorem to the case where there are two functions instead of a single one, which we here 
denote by /(/). According to a review in the Zentralblatt fiir Math., 19, p. 224, Khint- 
chine in the vsull. Acad. Sci. URSS, sir. Math. Nr. 3 (1938), 313-322, has continued and 
made precise the earlier work of Rowland. 
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expanding 
o 
exp >, (iu)"),/n! 
n=3 


as a power series in u, integrating termwise using 
| +90 at 
° nl ° to 1 ) 
~— (iuo)" exp | —inox — du = (—)"o ¢'"’(x), 
2r - 00 2 


n) 


1 d” 
a (xX) = — “aa n€ 
V 2x dx 


" . ‘ . ° ° 1/2 
and finally coilecting terms according to their order in powers of » 


et ta 


r} —_ 
» Zlve 


- ; — ~4 P he —5 2 ~7 itt 
PU) ~ oo '(x) - = y* a) +] 5 yo” (x) 4. 36 ea |+ vee 


3! 


(1.6-3 


The first term is O(v"”), the second term is O(v" 7. and the term within 
brackets is O(v°*). This is Edgeworth’s series.’ The first term gives the 
normal distribution and the remaining terms show how this distribution is 
approached as y—> x. 


1.7 THe FourteER CoMPpoNENTs OF [(¢) 


In some analytical work noise current is represented as 


‘ 

ay 2rnt . 2ant i‘ 

I(t) = = + > (@, cos —— + 8, sin —— (1.7-1 
2 nel 7 T 

where at a suitable place in the work T and JN are allowed to become infinite 

The coefficients a, and 6, ,1 <n < N, are regarded as independent random 

variables distributed about zero according to a normal law. 

It appears that the association of (1.7—1) with a sequence of disturbances 

. - 16 . 
occurring at random goes back many years. Rayleigh’ and Gouy suggested 
that black-body radiation and white light might both be regarded as se- 
quences of irregularly distributed impulses. 

Einstein’ and von Laue have discussed the normal distribution of the 
coefficients in (1.7-1) when it is used to represent black-body radiation, this 
radiation being the resultant produced by a great many independent os- 

5 See, for example, pp. 86-87, in ‘Random Variables and Probability Distributions” 
by H. Cramér, Cambridge Tract No. 36 (1937). 

6 Phil. Mag. Ser. 5, Vol. 27 (1889) pp. 460-469 

7 A. Einstein and L. Hopf, Ann. d. Physik 33 (1910) pp. 1095-1115. 

M. V. Laue, Ann. d. Physik 47 (1915) pp. 853-878. 
A. Einstein, Ann. d. Physik 47 (1915) pp. 879-885. 
M. V. Laue, Ann. d. Physik 48 (1915) pp. 668-680. 
I am indebted to Prof. Goudsmit for these references. 
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cillators. Some argument arose as to whether the coefficients in (1.7-1 
were Statistically independent or not. It was finally decided that they 
are independent. 

Che shot effect current has been represented in this way by Schottky.” 
[he Fourier series representation has been discussed by H. Nyquist’ and 
also by Goudsmit and Weiss. Remarks made by A. Schuster'’ are equiv- 
alent to the statement that a, and 6, are distributed normally. 

In view of this wealth of information on the subject it may appear super- 
fluous to say anything about it. However, for the sake of completeness, 
we shall outline the thoughts which lead to (1.7-1). 

In line with our usual approach to the shot effect, we suppose that exactly 
K electrons arrive during the interval (0, 7), so that the noise current for 
the interval is 
K 


I(t) = > F(t — &) (1.7-2) 


k=] 


The coefficients in the Fourier series expansion of /x(t) over the interval 
(0, 7) are dnx and ba, where 


. a sf dant 
Qux — tnx > | F(t — &) exp| —2 —~. | dt 
T =i 40 7 
2 ? [- | dan | 
oer F(t) exp | —i —— (t + tk) | dt 
Te Le I 


to 


Il 


| 
| 
x 
is) 
! 
6 
3 
& 
= 
~ 
—_ 
! 
ww 


In this expression 


2ri, 
= — 
7 1 
(1.7-4) 
a dae ee 
Re" = C, — iS, = | F(t)" dt 
T dw 
In the earlier sections the arrival times ¢; , fe, --- ¢« were regarded as A 


independent random variable each distributed uniformly over the interval 
(O, T). Hence the 6;,’s may be regarded as random variables distributed 
uniformly over the interval 0 to 2r. 

Incidentally, it will be noted that in (1.7-3) there occurs the sum of A 
randomly oriented unit vectors. When A becomes very large, as it does 

5 Ann. d. Physik, 57 (1918) pp. 541-507. 

® Unpublished Memorandum, ‘Fluctuations in Vacuum Tube Noise and the Like,” 
March 17, 1932. 


10 Investigation of Hidden Periodicities, Terrestrial Magnetism, 3 (1898), pp. 13-41 
See especially propositions 1 and 2 on page 26 of Schuster’s paper 
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when vy > ©, it is known that the real and imaginary parts of this sum ar 
random variables, which tend to become independent and normally < 
tributed about zero. This suggests the manner in which the normal « 
tribution of the coefficients arises. Averaging over the 6,’s in (1.7-3) giv 
when » > 0 


Gnx = bx = 0 (1.7-5 
Some further algebra gives 


~ = 2 
2 


Qnk = bar = 


ax On = G@:xnQmk = bik bmx = 0 


where m * mand n,m > 0. 
So far, we have been considering the case of exactly A arrivals in our 


interval of length 7. Now we pass to the general case of any number of 


arrivals by making use of formulas analogous to 


™ ed . 
an = D1 p(K)ank (1.7-7 
k=0 
as has been done in section 1.3. Thus, for x > 0, 
d, = b, = 0 
22 ¢ 9 9 
a,=b, =" Ri =o, (1.7-8 
2 
Gab, = Gadn = 6,6, = 0, n~m 


In the second line we have used g, to denote the standard deviation of a, 
r . Ss ° 

and 6,. We may put the expression for «, in a somewhat different form 

by writing 


Se = 7 nf, Af = 1 (1.7-9) 


where f, is the frequency of the mth component. Using (1.7-4), 


| +00 2 
on = 2vAf [ Fe?" dt (1.7-10) 
| °o | 


Thus, o» is proportional to »/T. 

The probability density function P(a; , «-+ ay, b;, «++ by) for the 2N co 
efficients, a), +++ aw, b:, +++ by may be derived in much the same fashion 
as was the probability density of the noise current in section 1.4. Here V 
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s arbitrary but fixed. The expression analogous to (1.4-5) is the 2N fold 
integral 


P(a,,-++, by) = (2) | du, --- [ doy (1.7-11) 


exp [—i(aim + +--+ + byty) — »T + vTE) 


where 


2x N 
E= =f dé exp E +» (tnCn + 0,5,) cos nO + (tC, — u,S,) sin nt | 
T 40 n=l 
(1.7-12) 


in which C, — iS, is defined as the Fourier transform (1.7-4) of F(t). 

The next step is to show that (1.7-11) approaches a norma! law in 2N di- 
mensions as y—> ©. ‘This appears to be quite involved. It will be noted 
that the integrand in the integral defining E is composed of N factors of the 
iorm 


exp [ip, cos (n8 — Yp)] 
= Jo(pn) + 2i cos (m8 — Yn)Ji(pn) — 2 cos (2n6 — 2p,)J2(pn) + °° 


where 
2 2 2 2 2 2 2 2 2 
Pn = (un + MC. + Sia) = wT on(Un + Un). 


As v becomes large, it turns out that the integral (1.7-11) for the prob- 
ability density obtains most of its contributions from small values of u and ». 
By substituting the product of the Bessel function series in.the integral for 
E and integrating we find 


E= II Ido.) +A+B+C 


where A is the sum of products such as 
—2i cos (Wirt — We — WOSi(px)Ji(pt)Ji(pere) times N — 3 Jo’s 


inwhichO << k<land2<k+1< N. Similarly Bis the sum of products 
of the form 


—2i cos (Wox ~~ We) J 1( pox) J2( px) times N — 2 Jo's 


C consists of terms which give fourth and higher powers in wand v. There 
are roughly V’/4 terms of form A and N/2 terms of form B. 
Expanding the Bessel functions, neglecting all powers above the third and 
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proceeding as in section 1.4, will give us the normal distribution plus the firs! 
correction term. It is rather a messy affair. An idea of how it looks ma 
be obtained by taking the special case in which F(¢) is an even function ci 
and neglecting terms of type B. Then 


~(r2+y2)/ 
(r2+y2)/2 


P(a,, ees dy, by, »++ by) - (1 + n) II‘ 


2 
n=1 Qra nr 


(1.712 


where 


dn _ b, 
Xn —=- , Vn weeny) 
on on 


n = (207)? Do [xee(aeae — ye ve) + 2 vertyevel (17-13 
ke 


and the summation extends over 2 << k +1 < Nwithk < 1. 

It is seen that if JT and WN are held constant, the correction term 7 ap- 
proaches zero as v becomes very large. A very rough idea of the magnitude 
of n may be obtained by assuming that unity is a representative value of the 
x’s and y’s. Further assuming that there are N* terms in the summation 
each one of which may be positive or negative suggests that magnitude of 
the sum is of the order of V. Hence we might expect to find that 7 is of 
the order of N(2vT)~"”. 


PART II 
POWER SPECTRA AND CORRELATION FUNCTIONS 


2.0 INTRODUCTION 


A theory for analyzing functions of time, ¢, which do not die down and 
which remain finite as ¢ approaches infinity has gradually been developed 
over the last sixty years. A few words of its history together with an ex- 
tensive bibliography are given by N. Wiener in his paper on “Generalized 
Harmonic Analysis’”’."’ G. Gouy, Lord Rayleigh and A. Schuster were led 
to study this problem in their investigations of such things as white light 
and noise. Schuster’ invented the “periodogram” method of analysis which 
has as its object the discovery of any periodicities hidden in a continuous 
curve representing meteorological or economic data. 

" Acta Math., Vol. 55, pp. 117-258 (1930). See also “Harmonic Analysis of Irregular 
Motion,” Jour. Math. and Phys. 5 (1926) pp. 99-189. 

2 The periodogram was first introduced by Schuster in reference 10 cited in Section 
1.7. He later modified its definition in the Trans. Camb. Phil. Soc. 18 (1909), pp. 107- 
135, and still later redefined it in ‘‘The Periodogram and its Optical Analogy,” Proc. Roy. 
Soc., London, Ser. A, 77 (1906) pp. 136-140. In its final form the periodogram is equiva- 


lent to 4w(f), where w(f) is the power spectrum defined in Section 2.1, plotted as a func- 
tion of the period T = (2xf)7. 
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The correlation function, which turns cut to be a very useful tool, appar- 
ently was introduced by G. I. Taylor.’ Recent!y it has been used by quite 
. few writers’ in the mathematical theory of turbulence. 

In section 2.1 the power spectrum and correlation function of a specifi 
function, such as one given by a curve extending tof = «, are defined by 
equations (2.1-3) and (2.1-4) respectively. That they are related by the 
Fourier inversion formulae (2.1-5) and (2.1-6) is merely stated; the dis- 
cussion of the method of proof being delayed until sections 2.3 and 2.4. In 
section 2.3 a discussion based on Fourier series is given and in section 2.4 a 
parallel treatment starting with Parseval’s integral theorem is set forth. 
The results as given in section 2.1 have to be supplemented when the func- 
tion being analyzed contains a d.c. or periodic components. This is taken 
up in section 2.2. 

The first four sections deal with the analysis of a specific function of ¢. 
However, most of the applications are made to functions which behave as 
though they are more or less random in character. In the mathematical 
analysis this randomness is introduced by assuming the function of ¢ to be 
also a function of suitable parameters, and then letting these parameters be 
random variables. This question is taken up in section 2.5. In section 2.6 
the results of 2.5 are applied to determine the average power spectrum and 
the average correlation function of the shot effect current. The same thing 
is done in 2.7 for a flat top wave, the tops (and bottoms) being of random 
length. The case in which the intervals are of equal length but the sign 
of the wave is random is also discussed in 2.7. The representation of the 
noise current as a trigonometrical series with random variabie coefficients 
is taken up in 2.8. The last two sections 2.9 and 2.10 are devoted to prob- 
ability theory. The normal law and the central limit theorem, respectively, 
are discussed, 


2.1 Some RESULTS OF GENERALIZED HARMONIC ANALYSIS 


We shall first state the results which we need, and then show that they are 
plausible by methods which are heuristic rather than rigorous. Suppose 
that J(¢) is one of the functions mentioned above. We may think of it as 
teing specified by a curve eatending from’ = —x tot = x, IJ(t) may 
be regarded as composed of a great number of sinusoidal components whose 
frequencies range from O0to+«. It does not necessarily have to be a noise 
current, but if we think of it as such, then, in flowing through a resistance of 
one ohm it will dissipate a certain average amount of power, say p watts. 

3 Diffusion by Continuous Movements, Proc. Lond. Math. Soc., Ser. 2, 20, pp. 196 
212 (1920). 


4 See the text ““Modern Developments in Fluid Dynamics” edited by S. Goldstei 
Oxford (1938). 
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That portion of p arising from the components having frequencies betwe: 
fand f + df will be denoted by w(f)df, and consequently 


p= [ wna (21-1 


Since w(f) is the spectrum of the average power we shall call it the “‘powe: 
spectrum” of J(t). It has the dimensions of energy and on this account is 
frequently called the ‘“energy-frequency spectrum” of /(¢). A mathematica 
formulation of this discussion leads to a clear cut definition of w/(f). 

Let &(t) be a function of ¢, which is zero outside the interval 0 <¢ < Tand 
is equal to /(t) inside the interval. Its spectrum S(f) is given by 


T 
sin = [ moet at (2.1-2 
0 
The spectrum of the power, w(f), is defined as 
/ 2 
w(f) = Limit 2| sf| (2.1-3 
T-2 Zz. 


where we consider only values of f > 0 and assume that this limit exists. 
This is substantially the definition of w(f) given by J. R. Carson” and is 
useful when J(t) has no periodic terms and no d.c. component. In the 
latter case (2.1-3) must either be supplemented by additional definitions or 
else a somewhat different method of approach used. These questions will 
be discussed in section 2.2. 

The correlation function ¥(r) of I(¢) is defined by the limit 


1? 
¥(r) = Limit rl T(t)I(t + 1) dt (2.1-4 


which is assumed to exist. (7) is closely related to the correlation coeffi- 
cients used in statistical theory to measure the correlation of two random 
variables. In the present case the value of /(¢) at time / is one variable and 
its value at a different time ¢ + 7 is the other variable. 

The spectrum of the power w(f) and the correlation function ¥(7) are 
related by the equations 


wf) = a ¥(r) cos 2xfr dr (2.1-5) 


¥(r) = I w(f) cos 2nfr df (2.1-6) 


46 “The Statistical Fnergy-Frequency Spectrum of Random Disturbances,” B.S.7.J 
Vol. 10, pp. 374-381 (1931). 
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—~ 
w 


[It is seen that ¥(r) is an even function of r and that 
¥(0) = p (2.1-7) 


When either ¥(r) or w(f) is known the other may be obtained provided the 
corresponding integral converges. 


2.2 PoweR SPECTRUM FOR D.C. AND PERIODIC COMPONENTS 


As mentioned in section 2.1, when J(¢) has a d.c. or a periodic com- 
ponent the limit in the definition (2.1-3) for w(f) does not exist for f equal 
to zero or to the frequency of the periodic component. Perhaps the most 
satisfactory method of overcoming this difficulty, from the mathematical 
point of view, is to deal with the integral of the power spectrum.” 


ho 


f 
I w(g) dg (2.2-1) 


instead of with w(f) itself. 
The definition (2.1-4) for (7) still holds. If, for example, 
I(t) = A +C cos (2xfot — ¢) (2.2-2) 
¥(r) as given by (2.1-4) is 


2 


¥(r) = A’ + cos 2nfor (2.2-3) 
The inversion formulas (2.1-5) and (2.1-6) give 


f . in af 
[ we dg =? [ yo) 2 as 
0 7 JO T 


© f 
| cos 2nfr i{ [ we) ds 
0 0 


16 This is done by Wiener," loc. cit., and by G. W. Kenrick, “The Analysis of Irregular 
Motions with Applications to the Energy Frequency Spectrum of Static and of Telegraph 
Signals,” Phil. Mag., Ser. 7, Vol. 7, pp. 176-196 (Jan. 1929). Kenrick appears to be one 
of the first to apply, to noise problems, the correlation function method of computing the 
power spectrum (one of his problems is discussed in Sec. 2.7). He bases his work on re- 
sults due to Wiener. Khintchine, in “Korrelationstheorie der station iren stochastischen 
Prozesse,” Math. Annalen, 109 (1934), pp. 604-615, proves the following theorem: A 
necessary and sufficient condition that a function R(¢) may be the correlation function of 
a continuous, stationary, stochastic process is that R(t) may be expressed as 


(2.2-4) 


Il 


¥(r) 


+2 
R(t) = [ cos tx dF(x) 
2 
where F(x) is a certain distribution function. This expression for R(t) is essentially the 
second of equations (2.2-4). Khintchine’s work has been extended by H. Cran:r, “On 
the theory of stationary random processes,” Ann. of Mut’, Ser. 2, Vol. 41 (194)), pp. 
215-230. However, Khintchine and Cram$r appear to be interested pri narily in ques- 
tions of existence, representation, etc., and do not stress the concept of the power spectrum. 
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where the last integral is to be regarded as a Stieltjes’ integral. When t}, 
expression (2.2-3) for ¥(7) is placed in the first formula of (2.2—4) we get 


[ ( A? : when 0 <f < fo 
w(g) dg = “re, 
“0 (s é 4 + 2 ’ “ f > fo 


When this expression is used in the second formula of (2.2-4), the increments 
of the differential are seen to be A’ at f = Oand C’/2at f = fo. The re- 
sulting expression for ¥(7) agrees with the original one. 

Here we desire to use a less rigorous, but more convenient, method of 
dealing with periodic components. By examining the integral of w(f) as 
given by (2.2-5) we are led to write 





“—-_- 
to 
bo 


te 


' Cc 
w(f) = 24*6(f) + = 8(f — fo) (2.2-6 
where 6(x) is an even unit impulse function so that if e > 0 
ne 1 € ? " x 
5(x) dx = ; 5(x) dx = 3 (2.2-7) 
0 2 € 
and 4(«) = 0 except at x = O, where it is infinite. This enables us to use 


the simpler inversion formuias of section 2.1. The second of these, (2.1-6), 
is immediately seen to give the correct expression for ¥(r). The first one, 
(2.1-5), gives the correct expression for w(f) provided we interpret the in- 
tegrals as follows: 


| cos 2afr dr = 48(f) 

Qa (2.2-8) 
| cos 2nfor cos 2afr dr = 45(f — fo) 

0 


It is not hard to show that these are in agreement with the fundamental! 
interpretation 


+0 =. 
[ el ay = / elt dy = 8(f) (2.2-9) 


which in its turn follows from a formal application of the Fourier integra! 
formula and 


+00 +0 ‘ 
| xine ay = [ i(fye or df =1 (2.2-10 


“20 


We must remember that fo > 0 and f > 0 in (2.2-8) so that 6(f + fo) = 0 
for f > 0. 
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The definition (2.1-3) for w(f) gives the continuous part of the power 
spectrum. In order to get the part due to’ the d.c. and periodic com- 
ponents, which is exemplified by the expression (2.2-6) for w(f) involving 
the 6 functions, we must supplement (2.1-3) by adding terms of the type 


2A73(f) + © caf = fo) = = | Limit 1? ane, 


v [Limit ? te) L “ay — fo) 


(2.2-11) 


The correctness of this expression may be verified by calculating S(f) for 
the J(¢) of this section given by (2.2-2), and actually carrying out the limiting 
process. 

2.3. Discussion OF RESULTS OF SECTION ONE—FOURIER SERIES 

The fact that the spectrum of power w(f) and the correlation function 
y(r) are related by Fourier inversion formulas is closely connected with 
Parseval’s theorems for Fourier series and integrals. In this section we shall 
not use Parseval’s theorems explicitly. We start with Fourier’s series and 
use the concept of each component dissipating its share of energy inde- 
pendently of the behavior of the other components. 

Let that portion of J(t) which lies in the interval 0 < ¢ < T be expanded 
in the Fourier series 


. 
I(t) = 5 + > (« n COS aunt +b,s in 27) (2.3-1) 


n=l 


where 


Z ie 
a, = af T(t) « ‘s dt 


(2.3-2) 
t ? 
a ; I I(t) sin ie dt 
Then for the interval —r <¢ < T — 7, 
9 ? 
I(t+7r) = = 5 + = («. cos = 7) + $b, sin = a ?) (2.3-3) 
n=l 


Multiplying the series for 7(¢) and 7(t + +r) together and integrating with 
respect tod gives, after some reduction, 


T 
rl I(0)I(t + 1) dt 
; _ (2.3-4) 
-%+ 5) a’, + 6%) cos _ ™, +0(7) 


n=l 
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where the last term represents correction terms which must be added be 
cause the series (2.3-3) does not represent /(/ + 7) in the interval (T — +, 7 
when 7 > 0, or in the interval (0, — r) if 7 < 0. 

If /(¢) were a current and if it were to flow through one ohm for the ir 
terval (0, T), each component would dissipate a certain average amount 
of power. The average power dissipated by the component of frequenc\ 
fn = n/T cycles per second would be, from the Fourier series and elementary 
principles, 


*(a’, + b%.) watts, n #0 
. (2.3-5) 
; watts, n= 0 


The band width associated with the mth component is the difference in 
frequency between the » + 1 th and mth components: 


Imi Sf. = a 7 Ps 


Hence if the average power in the band f, f + df is defined as w(f)df, the 
average power in the band fnui — fy is 


w(fn) Inst — fa) = () ; 


and, from (2.3-5), this is given by 


v(z)4 =i. +8), #0 


7a. = 
1 : (2.3-6) 
= # o 
w(0) + q? % = (0 


When the coefficients in (2.34) are replaced by their expressions in terms 
of w(f) we get 


ew tl 

| I() I(t + 2) dt + o( 7) 
= n Qunr 
d ‘ (3) sk 


[Gores 
-[ 


(2.3-7) 


w(f) cos 2xfr df 
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here we have assumed T so large and w(f) of such a nature that the summa- 

on may be replaced by integration. 

If J remains finite, then as T — © with 7 held fixed, the correction term 
on the left becomes negligibly small and we have, upon using the definitions 

2.1-4) for the correlation function ¥(7), the second of the fundamental 
inversion formulas (2.1-6). The first inversion formula may be obtained 
from this at once by using Fourier’s double integral for w(f). 

Incidentally, the relation (2.3-6) between w(f) and the coefficients a, and 
6, is in agreement with the definition (2.1-3) for w(f) as a limit involving 
| S(f) | *. From the expressions (2.3-2) for a, and 6,, the spectrum S(/,) 
given by (2.1-2) is 


SU.) = 5 (ay — ib,) 


Then, from (2.1-3) w(f,) is given by the limit, as T— ~, of 


2 


oo 2 2T 72 2 
71 SU) | = T 4 (a, + b;,) 


(a, + 03) 


mis 


og n\ . 
and this is the expression for w (2) given by (2.3-6). 


2.4 DISCUSSION OF RESULTS OF SECTION ONE—PARSEVAL’S THEOREM 


The use of Parseval’s theorem” enables us to derive the results of section 
2.1 more directly than the method of the preceding section. This theorem 
states that 


[ Fi(f)Fi(f) df = [ G,(t)G2(—2) dt (2.4-1) 
where F;, G; and F2, G2 are Fourier mates related by 
+00 
F(/) = CeO" de 
oe (2.4-2) 
G(t) = F(fye"*"* df 


It may be proved in a formal manner by replacing the F; on the left of 
(2.4-1) by its expression as an integral involving G,(t). Interchanging the 


17 E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals, Oxford (1937). 








318 BELL SYSTEM TECHNICAL JOURNAL 


order of integration and using the second of (2.4-2) to replace F2 by Ge giv: 
the right hand side. 

We now set Gi(¢t) and Go(t) equal to zero except for intervals of length 7 
These intervals and the corresponding values of G; and Ge are 


G,(t) = I(t), ie ee (2.4-3 
G2(t) 


ll 


I(-t+7r), 7 -T<t<r 
From (2.4 3) it follows that Fi(f) is the spectrum S(f) of Z(¢) given by equ: 
tion (2.1-2). Since /(¢) is real it follows from the first of equations (2.42 
that 

S(-f) = S*(f), (2.4-4 
where the star denotes conjugate complex, and hence that | S(f) | ? is an 
even function of f. 

The first of equations (2.42) also gives 


F.(f) = | I(—t + re 7°" de 


eT ~ 
= | (te? dt (2.4-: 
0 
-_ sinter” 
When these G’s and F’s are placed in (2.4-1) we obtain 
+-00 T-t 
| S(f) |°e?*"" df = | I(t)T(t + 1) dt (2.4-6 
- 0 


where we have made use of the fact that G2(--/) is zero except in the inter- 
val —r <¢ < T — rand have assumed 7 > 0. If + < O the limits of 
integration on the right would be —7 and T. 

Since | S(f) |’ is an even function of f we may write (2.4-6) as 


T 2 Co) 2 | 7 12 
; | T(t)I(t + r) dt + 0(7) * - a cos 2efrdf (24-7) 


If we now define the correlation function (7) as the limit, as T— «, of the 


left hand side and define w(f) as the function 


_. Sis 
w(f) = Limit ‘ f>0 (2.1-3) 


T-2 


we obtain the second, (2.1-6), of the fundamental inversion formulas. As 
before, the first may be obtained from Fourier’s integral theorem. 
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In order to obtain the interpretation of w(f)df as the average power dis 
sipated in one ohm by those components of J(t) which lie in the band f, 
f+ df, we set 7 = Oin (2.4-7): 


1 T @ 

Limit — [ I(t) dt = | wf) df (2.4-8 
roe T do J 

rhe expression on the left is certainly the total average power which would 

be dissipated in one ohm and the right hand side represents a summation 

over all frequencies extending from 0 to x. It is natural therefore to in 

terpret w(f)df as the power due to the components in f, f + df. 

The preceding sections have dealt with the power spectrum w(/) and corre- 
lation function ¥(r) of a very general type of function. It willbe noted 
that a knowledge of w(f) does not enable us to determine the original func- 
tion. In obtaining w(f), as may be seen from the definition (2.1—3) or from 
(2.3-6), the information carried by the phase angles of the various compo 
nents of 7(¢) has been dropped out. In fact, as we may see from the Fourier 
series representation (2.3-1) of 7(¢) and from (2.3-6), it is possible to obtain 
an infinite number of different functions all of which have the same w* / 
and hence the same y(7). All we have to do is to assign different sets of 
values to the phase angles of the various components, thereby keeping 
a, + 6% constant. 


2.5 HARMONIC ANALYSIS FOR RANDOM FUNCTIONS 


In many applications of the theory discussed in the foregoing sections 
I(t) is a function of ¢ which has a certain amount of randomness associated 
with it. For example /(¢) may be a curve representing the price of wheat 
over a long period of years, a component of air velocity behind a grid placed 
in a wind tunnel, or, of primary interest here, a noise current. 

In some mathematical work this randomness is introduced by considering 
I(t) to involve a number of parameters, and then taking the parameters to 
be random variables. Thus, in the shot effect the arrival times /; , fe , -- + tx 
of the electrons were taken to be the parameters and each was assumed to be 
uniformly distributed over an interval (0, 7). 

For any particular set of values of the parameters, /(¢) has a definite power 
spectrum w(f) and correlation function ¥(r). However, now the principal! 
interest is not in these particular functions, but in functions which give the 
average values of w(f) and ¥(r) for fixed f and 7. These functions are ob 
tained by averaging w(f) and ¥(7) over the ranges of the parameters, using, 
of course, the distribution functions of the parameters. 

By averaging both sides of the appropriate equations in sections 2.1 and 
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2.2 it is seen that our fundamental inversion formulae (2.1-5) and (2.1-6 
are unchanged. Thus, 


wf) 


ll 


4 | ¥(r) cos 2xfr dr (2.5-1) 
0 


V(r) I w(f) cos 2xfr df (2.5-2) 


where the bars indicate averages taken over the parameters with f or 7 held 
constant. 


The definitions of ® and ¥ appearing in these equations are likewise ob- 


tained from (2.1-3) and (2.1-4) 


= 
salitie - ask (2.5-3) 
and 
7 Dw titans 
¥(r) = Limit | I(t)I(t + 7) dt (2.5-4) 
T-0 


The values of ¢ and ¢ are held fixed while averaging over the parameters. 
In (2.5-3) S(f) is regarded as a function of the parameters obtained from 
I(t) by 


S(f) = Te" de (2.1-2) 


Similar expressions may be obtained for the average power spectrum for 
d.c. and periodic components. All we need to do is to average the ex- 
pression (2.2~11) 

Sometimes the average vaiue of the product J(t)7(¢ + 7) in the definition 
(2.5-4) of P(r) is independent of the time T. This enables us to perform 
the integration at once and obtain 


V(r) = (OT + 7) (2.5-5) 


This introduces a considerable simplification and it appears that the simplest 
method of computing wf) for an J(¢) of this sort is first to compute ¥(r), and 
then use the inversion formula (2.5-1). 


2.6 First EXAMPLE—THE Snot EFFECT 


We first compute the average on the right of (2.5-5). By using the 
method of averaging employed many times in part I, we have 


I(t + 7) = p> P(K) Tx(OT x(t + 7) (2.6-1) 
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ere p(K) is the probability of exactly A electrons arriving in the inter- 


Val (0, T), 


a eh" ... 
p(K) = “zr ° , (1.1-3) 
ind 
I(t) = 2 F(t — &) (1.3-1) 
=i 


Multiplying 7x«(¢) and Jx(¢ + 7) together and averaging f1, f2, +++ tx 
over their ranges gives 


. - - £ ” dt T dtr - \ ry 
IOQIkt +7 =D] wees | AFG - WF +7 —t) 
k=1 m=1 40 7 0 7 


This is similar to the expression for J2(¢) which was used in section 1.3 to 

prove Campbell’s theorem and may be treated in much the same way. 

Thus, if ¢and/ + 7 lie between A and T — A, the expression above becomes 
? 


> 


=i FOF + 7) dt + Sm : | F(t) a] 


When this is placed in (2.6-1) and the summation performed we obtain 
an expression independent of T. Consequently we may use (2.5-5) and get 


+00 
y(r) = » | F()F(t + 7) di + I(t)’ (2.6-2) 
ep 
where we have used the expression for the average current 
. — 
I(t) = » | F(t) dt (1.3-4) 
— 20 
In order to compute #(/) frem (7) it is convenient to make use of the 
fact that ¥(r) is always an even function of 7 and hence (2.5-1) may also 
be written as 


wif) = 2 : ¥(r) cos 2xfr dr (2.6-3) 


Then 


+20 


w(f) = w | dt F(t) a dr F(t + 1) cos 2nfr 


+2 
ts 2 / I(t). cos 2nfr dr 
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+00 +2 
= 2» Real Part of [ dt F(t)e"""' / dt’ F(t)" 
x — 30 
9 inion ’ 
aa TH [ eo"! dr 
: x 
= 2v|s(f)|* + 27(8 &(f) (2.6-4 
In going from the first equation to the second we have written (’ = / 4 


and have considered cos 2xf7 to be the real part of the corresponding e\ 
ponential. In going from the second equation to the third we have set 


s(f) = | F(t)e?*"* dt (2.6-5 
and have used 


$20 
| et dt = 8(f) (2.2-9 
~ 2 

The term in @(f) involving 6(f) represents the average power which woul 
be dissipated by the d.c. component of J(¢) in flowing through one ohm 
It is in agreement with the concept that the average power in the bani 
0 < f < «,¢€ > 0 but very small, is 


| w(f) dt = 2 | 5(f) df 
0 0 (2.6-6 


= I(t) 

The expression (2.6-4) for @(f) may also be obtained from the definitio: 
(2.5-3) for ®(f)plus the additional term due to the d.c. component ob 
tained by averaging the expressions (2.2-11). We leave this as an exercise 
for the reader. He will find it interesting to study the steps in Carson’s: 
paper leading up to equation (8). Carson’s R(w) is related to our @(f) by 


w(f) = 2nR(w) 


and his f(iw) is equal to our s/f). 
Integrating both sides of (2.6-4) with respect to f from 0 to ~ and using 


P= q w(f) df 


. 


gives the result 


~ 
te 
| 
=! 


= 2» | | s(f) df (2.6-7 


18 Loc. cit. 
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[his may be obtained immediately from Campbell’s theorem by applying 
Parseval’s theorem. 

As an example of the use of these formulas we derive the power spectrum 
of the voltage across a resistance R when a current consisting of a great num- 
ber of very short pulses per second flows through R. Let F(t — f,) be the 
voltage produced by the pulse occurring at time 4. Then 

F(t) = Ret) 
where ¢(¢) is the current in the pulse. We confine our interest to relatively 
low frequencies such that we may make the approximation 


+20 


sf) = Re(t)e"™* dt 


to} 


= R | ¢(t) dt = Rq 


where q is the charge carried through the resistance by one pulse. From 
(2.6-4) it follows that for these low frequencies the continuous portion of 
the power spectrum for the voltage is constant and equal to 


w(f) = IwR*g’ = 21R’q (2.6-8) 


where I = vg is the average current flowing through R. This result is often 
used in connection with the shot effect in diodes. 

In the study of the shot effect it was assumed that the probability of an 
event (electron arriving at the anode) happening in d/ was vdt where v is the 
expected number of events per second. This probability is independent of 
the time ¢. Sometimes we wish to introduce dependency on time.” As an 
example, consider a long interval extending from 0 to T. Let the prob- 
ability of an event happening in ¢, ¢ + dt be Kp(t)dt where A is the average 
number of events during T and p(t) is a given function of ¢ such that 


[ p(t) dt =1 


For the shot effect p(t) = 1/T. 

What is the probability that exactly A events happen in J? As in the 
case of the shot effect, section 1.1, we may divide (0, 7) into N intervals 
each of length At so that NAt = T. The probability of no event happening 


in the first A? is 
1— Kp (5) At 


18 A careful discussion of this subject is given by Hurwitz and Kac in “Statistical 
Analysis of Certain Types of Random Functions.”” I understand that this paper will 
soon appear in the Annals of Math. Statistics 
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The product of N such probabilities is, as V — ©, At—0; 


ep| -* | p(t) a| = «* 
0 


This is the probability that exactly 0 events happenin T. In the same wa 
we are led to the expression 


ms ” Mea (2.6-9 


for the probability that exactly K events happen in T. 

When we consider many intervals (0, 7) we obtain many values of A and 
also many values of J measured ¢ seconds from the beginning of each interva! 
These values of J define the distribution of J at time¢. By proceeding as i: 
section 1.4 we find that the probability density of J is 


+00 4 
P(I,t) = bd | du exp| —iul + k | p(x) (e™""-? — 1) ax| 
Qn 0 


The corresponding average and variance is 
T 
i= x | p(x) F(t — x) dx 
0 


— if 
(Q-TIP = ral p(x) F(t — x) dx (2 6-10) 
0 


If S(f) is given by (2.1-2) and s(f) by (2.6-5) (assuming the duration of 
F(t) short in comparison with T) the average value of | S( f) |’? may be ob- 
tained by putting (1.3-1) in (2.1-2) to get 


Sx(f) = s(f) py Pade 


‘ — : 
Expressing Sx(f) Sx(f), where the star denctes conjugate complex, as a 
double sum and averaging over the /,’s, using p(t), and then averaging over 
the K’s gives 


{2 


IS) P = KI sf) ak + Kk l p(xe Prt ax! | (2.6-11) 


This may be used to compute the power spectrum from (2.5-3) provided 
p(x) is not periodic. If p(x) is periodic then the method of section 2.2 
should be used at the harmonic frequencies. If the fluctuations of p(f) are 
slow in comparison with the fluctuations of F(t) the second term within the 


brackets of (2.6-11) may generally be neglected since there are no values of 
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f which make both it and s(f) large at the same time. On the other hand, 


if both p(t) and F(t) fluctuate at about the same rate this term must be 


considered. 
2.7 SECOND EXAMPLE—RANDOM TELEGRAPH SIGNAL'® 


Let I(t) be equal to either a or —a so that it is of the form of a flat top 
wave. Let the intervals between changes of sign, i.e. the lengths of the 
tops and bottoms, be distributed exponentially. We are led to this dis- 
tribution by assuming that, if on the average there are u changes of sign per 
second, the probability of a change of sign in ¢, ¢ + dt is udt and is independ- 
ent of what happens outside the interval ¢, 4 + dt. From the same sort of 
reasoning as employed in section 1.1 for the shot effect we see that the 
probability of obtaining exactly A changes of sign in the interval (0, T) is 

K 
(K) = en (2.7-1) 

We consider the average value of the product 7(t)J(¢+ 7). This product 
isa’ if the two /’s are of the same sign and is —a’ if they are of opposite sign. 
In the first case there are an even number, including zero, of changes of sign 
in the interval (¢, 4 + 7), and in the second case there are an odd number of 
changes of sign. Thus 


Average value of 7(t)7(¢ + 7) (2.7-2) 


2 *“s* 
=a X probability of an even number of 
changes of sign in#,¢ + 7 


— a X probability of an odd number of 
changes of sign int, ¢ + 7 


The length of the interval under consideration is | f+ 7 —¢| = | 7 | seconds. 
Since, by assumption, the probability of a change of sign in an elementary 
interval of length Ad is independent of what happens outside that interval, 
it follows that the same is true of any interval irrespective of when it starts. 
Hence the probabilities in (2.7—2) are independent of ¢ and may be obtained 
from (2.7-1) by setting T = |7 |. Then (2.7-2) becomes, assuming r > 0 
for the moment, 


(I(t + 7) 


a'[p(0) + p(2) + p(4) + -*°] 
— a[p(1) + p(3) + A(5) + ***] 


sa cent _ ut (ur) ae) | 


1! 2! 
= gte™" (2.7-3) 


6 Kenrick, cited in Section 2.2. 
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From (2.5-5), this gives the correlation function for J(t) 
Wr) = ae" (2.7-4 


The corresponding power spectrum is, from (2.5-1), 


(f) = 4a [ e™" cos 2xfr dr 
0 





Correlation functions and power spectra of this type occur quite fre 
quently. In particular, they are of use in the study of turbulence in hydro 
dynamics. We may also obtain them from our shot effect expressions if we 
disregard the d.c. component. All we have to do is te assume that the 
effect F(t) of an electron arriving at the anode at time ¢ = 0 is zero for 
¢ < 0, and that F(t) decays exponentially with time after jumping to its 
maximum value at ¢ = 0. This may be verified by substituting the value 


F(t) = 2a 4/Be™ t>0 (2.7- 


for F(t) in the expressions (2.6-2) and (2.6-4) (after using 2.6-5) for the 
correlation function and energy spectrum of the shot effect. 

The power spectrum of the current flowing through an inductance and a 
resistance in series in response to a very wide band thermal noise voltage is 
also of the form (2.7-5). 

Incidentally, this gives us an example of two quite different J(¢)’s, one a 
flat top wave and the other a shot effect current, which have the same corre- 
lation functions and power spectra, aside from the d.c. component. 

There is another type of random telegraph signal which is interesting to 
analyze. The time scale is divided into intervals of equal length hk. In an 
interval selected at random the value of /(¢) is independent of the values in 
the other intervals, and is equally likely to be +a or —a. We could con- 
struct such a wave by flipping a penny. If heads turned up we would set 
I(t) =ain0<t<h. If tails were obtained we would set 7(#) = —a in 
this interval. Flipping again would give either +a or —a for the second 
interval h < t < 2h, andsoon. This gives us one wave. A great man\ 
waves may be constructed in this way and we denote averages over these 
waves, with ¢ held constant, by bars. 

We ask for the average value of J(#)](¢ + 7), assuming 7 > 0. First 
we note that if r > A the currents correspond to different intervals for a!| 
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values of ¢t. Since the values in these intervals are independent we have 
Qitt +7) = Tt) t+ 7) = 0 
for all values of ¢ when 7 > h. 

To obtain the average when 7 < / we consider ¢ to lic in the first interval] 
0 <t <h. Since all the intervals are the same frem a statistical point 
of view we lose no generality in doing this. Ift+ 7 < h,ie.,t << h — 7, 
both currents lie in the first interval and 


T@l(t + 7) = a 
Ift > h — + the current /(¢ + 1) corresponds to the second interval and 
hence the average value is zero. 


We now return to (2.5-4). The integral there extends from 0 to T. 
When 7 > &, the integrand is zero and hence 


¥(r) = 0, r>h ($3.7) 


When 7 < 4, our investigation of the interval 0 < ¢ < enables us to write 
down the portion of the integral extending from 0 to /: 


I I(t)I(t + r) dt =| a’ dt + | O dt 
0 0 h—r 


=a (hk — 7) 


Over the interval of integration (0, 7) we have T/h such intervals each 
contributing the same amount. Hence, from (2.5-—4), 


Wr) = Limit %-— (h — 1) 
T-20 T h 
(2.7-8) 
=a'(1-1), 0O<sr<h 
h 
The power spectrum of this type of telegraph wave is thus 
9 . T 
w(f) = 42° | (1 = ‘) cos 2rfr dr 
; 0 h 
(2.7-9) 


a sin rfh\ 
= 2); Je 
a ( mh ) 


This is seen to have the same general behavior as #(f) for the first type 
of telegraph signal given by (2.7-5), when we relate the average number, 
u, of changes of sign per second to the interval length / by wh = 1. 
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2.8 REPRESENTATION OF NOISE CURRENT 


In section 1.8 the Fourier series representation of the shot effect curren 
was discussed. This suggests the representation* 


I(t) = ys (Gn COS wat + by sin wt) (2.8-1 


where 


On = 2tfn ; Ja a nAf (2.8-2) 


a, and b, are taken to be independent random variables which are distributed 
normally about zero with the standard deviation V wl fndf- w(f) is the 
power spectrum of the noise current, i.e., w(f) df is the average power which 
would be dissipated by those components of [(/) which lie in the frequency 
range f, f + df if they were to flow through a resistance of one ohm. 

The expression for the standard deviation of a, and b, is obtained when 
we notice that Af is the width of the frequency band associated with the nth 
component. Hence w(f,)Af is the average energy which would be dissi- 
pated if the current 


An COS Wal + b, Sin wyl 


were to flow through a resistance of one ohm, this average being taken over 
all possible values of a, and 6,. Thus 


2 2 a . mY ® 2 “2 PH . 
W(frAf = Gn COS Wal + 24,5, COS Wal SiN wl + bpSin wt =a,=b, (2.8-3 


The last two steps follow from the independence of a, and b, and the identity 
of their distributions. It will be observed that w(/), as used with the repre- 
sentation (2.8-1), is the same sort of average as was denoted in section 2.5 
by w(f). However, w(f) is often given to us in order to specify the spectrum 
of a given noise. 

For example, suppose we are interested in the output of a certain filter 
when a source of thermal noise is applied to the input. Let A(f) be the 
absolute value of the ratio of the output current to the input current when a 
steady sinusoidal voltage of frequency f is applied to the input. Then 


w(f) = cA) (2.8-4) 


* As mentioned in section 1.7 this sort of representation was used by Einstein and 
Hopf for radiation. Shottky (1918) used (2.8-1), apparently without explicitly taking 
the coefficients to be normally distributed. Nyquist (1932) derived the normal distribu 
tion from the shot effect. 
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li W is the average power dissipated in one ohm by /(f), 


T Cs) 
W = Limit =f P(t) dt = [ w(f) df 
T+0 1 0 0 


c [ A’(f) df 


which is an equation to determine c when W and A(f) are known. 

In using the representation (2.8-1) to investigate the statistical properties 
of I(t) we first find the corresponding statistical properties of the summation 
on the right when the a’s and 0’s are regarded as random variables distrib- 
uted as mentioned above and / is regarded as fixed. In general, the time 
{ disappears in this procedure just as it did in (2.8-3). We then let V—- « 
and Af —> 0so that the summations may be replaced by integrations. Fi- 
nally, the frequency range is extended to cover all frequencies from 0 to ~. 

The usual way of looking at the representation (2.8-1) is to suppose that 
we have an oscillogram of /(/) extending from’ = Otot = ~. This oscil- 
logram may be cut up into strips of length 7. A Fourier analysis of I(t) 
for each strip will give a set of coefficients. These coefficients will vary 
from strip to strip. Our representation (7Af = 1) assumes that this varia- 
tion is governed by a normal distribution. Our process for finding sta- 
tistical properties by regarding the a’s and 6’s as random variables while ¢ 
is kept fixed corresponds to examining the noise current at a great many 
instants. Corresponding to each strip there is an instant, and this instant 
occurs at ¢ (this is the / in (2.8-1)) seconds from the beginning of the strip. 
This is somewhat like examining the noise current at a great number of 
instants selected at random. 

Although (2.8-1) is the representation which is suggested by the shot 
effect and similar phenomena, it is not the only representation, nor is it 
always the most convenient. Another representation which leads to the 
same results when the limits are taken is” 


N 
I(t) = D> en cos (wat — ¢n) (2.8-6) 
n=l 


where ¢; , ¢2, °** ¢w are angles distributed at random over the range (0, 27) 
and 


Cn = [2w(frAfl?, wa = 20fn, fn = nf (2.8-7) 


19 This representation has often been used by W. R. Bennett in unpublished memoranda 
written in the 1930's. 
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In this representation J(¢) is regarded as the sum of a number of sinusoi: 
components with fixed amplitudes but randem phase angles. 

That the two different representations (2.8-1) and (2.8-6) of I(t) lead 
to the same statistical properties is a consequence of the fact that they ar 
always used in such a way that the “central limit theorem*” may be used 
in both cases. 

This theorem states that under certain general conditicns, the distribu- 
tion of the sum of V random vectors approaches a normal law (it may be 
normal in several dimensions**) as N—> «. In fact from this theorem it 
appears that a representation such as 

N 
I(t) = Do (an cos wnt + by sin wn) (2.8-6 
n=l 


where a, and b, are independent random variables which take only the values 
we 11/2 _ ° ° ° ° 

+ [w(f,)Af]”, the probability of each value being 3, will lead in the limit 

to the same statistical properties of 7(¢) as do (2.8-1) and (2.8-6). 


2.9 THe NORMAL DISTRIBUTION IN SEVERAL VARIABLES”" 


Consider a random vector r in A dimensions. The distribution of this 
vector may be specified by stating the distribution of the A components, 
x1, %2,°°' Xx, 0fr. ris said to be normally distributed when the prob 
ability density function of the x’s is of the form 

. -/2 1—1/2 7 
(27) *"| M |"? exp [—3x’M' x] (2.9-1 


where the exponent is a quadratic form in the x’s. The square matrix V/ 
is composed of the second moments of the x’s. 


Pi 6° mix] 


M=| - . (2.9-2) 
Mix °° uxx | 
where the second moments are defined by 
wu = x ; ty = aye , etc. (2.9-3 
| M | represents the determinant of M and 2’ is the row matrix 
x’ = [vy sey es Xx] (2.9-4 


x is the column matrix obtained by transposing x’. 


* See section 2.10. 

** See section 2.9. 

20H. Cramér, “Random Variables and Probability Distributions,” Chap. X., Cambridg 
Tract No. 36 (1937). 
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[he exponent in the expression (2.9-1) for the probability density may 
be written out by using 


J —. “ Mis hla @ 
Ms = } be Xr Xe (2.9-5) 


where M,, is the cofactor of u,, in M. 

Sometimes there are linear relations between the x’s so that the random 
vector r is restricted to a space of less than A dimensions. In this case the 
appropriate form for the density function may be obtained by considering 
a sequence of A-dimensiona] distributions which approach the one being 
investigated. 

If r, and r2 are two normally distributed random vectors their sum 7; + 72 
is also normally distributed. It follows that the sum of any number of 
normally distributed random vectors is normally distributed. 

The characteristic function of the normal distribution is 


K K 
dais 1 
t2z)7 Tizeret oe. izKnT f 
Ses KK = exp| — : > Dy bre Sr Ze (2.9-6) 


< r=l s=l 
2.10 CENTRAL Limit TimeoREM 

The central limit theorem in probability states that the distribution of the 
sum of N independent random vectors r; + re + --+ + ry approaches a 
normal law as VN — « when the distributions of 7 , r2 , «++ ry satisfy certain 
general conditions.’ 

As an example we take the case in which 7, , rz , «+ - are two-dimensional 

21 ‘ pe 6 . 

vectors’ , the components of r, being x, and v,. Without loss of generality 
we assume that 


The components of the resultant vector are 
X= Xt tees + ay 
Fant nat «++ + we 


and, since 7; , r2, *** are independent vectors, the second moments of the 
resultant are 


(2.10-1) 


wo = M=atagtees tay 
wn = Y= vit y2 + wes + yx (2.10-2) 
we = XV = ayy + xeye + +++ + ayy 


7 Incidentally, von Laue (see references in section 1.7) used this theorem in discussing 
the normal distribution of the coefficients in a Fourier series used to represent black-body 
radiation. He ascribed it to Markoff. ; 

*1 This case is discussed by J. V. Uspensky, “Introduction to Mathematical Probabil 
ity’, McGraw-Hill (1937) Chap. XV. 
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Apparently there are several types of conditions which are sufficient ' 
ensure that the distribution of the resultant approaches a normal law. One 
a ta ° 21 
sufficient condition is that 


2 
~3/2 — 
Kil i | Xn ° —0 
n=l 
. (2.10-3) 
—3/2 pa 
ya” DY [yn > 0 
n=l 
The central limit theorem tells us that the distribution of the random 
vector (X, Y) approaches a normal law as N — «. The second moments 
of this distribution are given by (2.10-2). When we know the second mo- 


ments of a normal distribution we may write down the probability density 
function at once. Thus from section 2.9 


Mf = @ = Mm a | M ry] Mo2 _ 
Mio M22 — (412 Miu 


| M| = uuu — nie 
x’ = [X, Y] 
a’Mo* x = | M|'(ueX* — 2weXVY + wu Y’) 
The probability density is therefore 
(wii Heo it uaz) ; exp |=" x* — Bu r S aon XY | (2.10-3 
2x 2(ui1 wee — M12) 


Incidentally, the second moments are related to the standard deviations 
a1 , 02 of X, Y and to the correlation coefficient rt of X and Y by 


2 
2 


iu = 01, Meo = O2, Mi2 = 70102 (2.10-4 


and the probability density takes the standard form 


(Q— 7)? 1 x? ‘ XY nm =) (2.10-5 
_> Tr ‘ — 10-5 
210102 - 21 — r’) 0} . 0102 02 


*1 This is used by Uspensky, loc. cit. Another condition analogous to the Lindeberg 
condition is given by Cramer,?° loc. cit. 


(To be concluded) 











ent | 


10-3) 


andom 
yments 
id mo- 
lensit) 


2.10-3 


‘jations 


2.10-4 


2.10-5 


indeberg 








Abstracts of Technical Articles by Bell System Authors 


Crossbar Toll Switching System... L. G. ABRAHAM, A. J. Buscu, AND 
I. F. SHIPLEY. A new crossbar toll switching system was placed in service 
in Philadelphia in August 1943. Important improvements offered by this 
system include: 

1. Transmission objectives are met more readily and substantial econo- 
mies are obtained in outside plant and in repeater equipment. 

2. Extended use of toll dialing results in operating economies and im- 
proved service to subscribers. Calls over circuits still on a ringdown basis 
are also handled more expeditiously and with operating economies. 

3. Flexibility due to the use of sender and markers to control establishing 
connections will be useful in meeting future toll switching requirements. 


Low-Frequency Quartz-Crystal Cuts having Low Temperature Coefficients? 
W. P. Mason AND R. A. Sykes. This paper discusses low-frequency, low- 
temperature-coefficient crystals which are suitable for use in filters and 
oscillators in the frequency range from 4 to 100 kilocycles. Two new cuts, 
the MT and NT, are described. These are related to the +5-degree X-cut 
crystal, which is the quartz crystal having the lowest temperature coefficient 
for any orientation of a bar cut from the natural crystal. When the width 
of the +5-degree X-cut crystal approaches in value the length, the motion 
has a shear component, and this introduces a negative temperature coefficient 
which causes the temperature coefficient of the crystal to become increas- 
ingly negative as the ratio of width to length increases. 

The MT crystal has its length along nearly the same axis as the +5-degree 
X-cut crystal, but its major surface is rotated by 35 to 50 degrees around the 
length axis. This results in giving the shear component a zero or positive 
temperature coefficient and results in a crystal with a uniformly low tempera- 
ture coefficient independent of the width-length ratio. A slightly higher 
rotation about the length axis results in a crystal which has a low-tempera- 
ture coefficient when vibrating in flexure and this crystal has been called the 
NT crystal. The NT crystal can be used in a frequency range from 4 to 50 
kilocycles, while the MT is useful from 50 kilocycles to 500 kilocycles. 

A special oscillator circuit is described which can drive a high-impedance 
NT flexure crystal. This circuit, together with the NT crystal, has been 
used to control the mean frequency of the Western Electric frequency- 
modulated radio transmitter. 

1 Elec. Engg., Transactions Section, June, 1944, 

*Proc. I, R. E., April 1944, 
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Electronics; Todey and Tomorrow.8 JoHN Mitts. John Mills, scientis: 
teacher, author, telephone and radio engineer, was first introduced 
the electron while still a fledgling under the tutelage of that eminent 
scientist R. A. Millikan at Chicago University. That was before Millika 
had, to quote the author, “first put salt on its tail.’’ Electronics, as a 
special branch of physics and engineering, has come to adulthood under the 
eyes of this observant author. From intimate association with this and 
allied fields John Mills writes knowingly. It is a book intended for the i: 
telligent layman rather than for the expert. 

“Electronics; Today and Tomorrow”’ likewise contains much of interest 
concerning the electronics of yesterday but, as one would expect, very little 
about the electronics of tomorrow because, as the author points out, much oi 
this ‘should await the victorious conclusion of the present conflict.” Th: 
book generally presents an interesting introduction to many things electron 
and throughout is interspersed with examples of present day techniques 
which employ electronic devices. All line drawings, diagrams and photo 
graphs have been omitted. 

The author begins his latest book with a brief capitulation of familiar 
engineering applications such as long-distance telephony, broadcasting, 
sound motion pictures and television which have been achieved through the 
use of electronic devices. He goes on to an historical account of some 
underlying discoveries and a discussion of atomic structure. An introduc 
tion to static electricity with the classical concepts of positive and negative 
charge is followed by establishment of the ideas of charges in motion, electri- 
cal current, discharges in gases, X-rays and their generation. 

The remainder of the book is divided into Part I—Electron Tubes, and 
Part If—Electronic Devices. In Part I the author follows the development 
of the art in nearly chronological order from diodes through the modern 
multi-electrode tubes giving uses and applications of each and explaining the 
purpose of the grids introduced. Part II discusses more complicated 
structures such as cathode-ray tubes, kinescopes, iconoscopes, electron 
microscopes, kodatrons, magnetrons, rumbatrons, klystrons, etc. and elabo- 
rates upon their practical applications. ‘This is, in fact, a wondrous field of 
developments. Finally there is a chapter on cyclotrons which the author 
says “‘comes into this book because it requires for its operation a powerfu! 
oscillator—or oscillator plus powerful electron amplifier—which will suppl) 
a high voltage at a frequency of megacycles.”” While many industrial ap 
plications, such perhaps as tin plating, might be included on the same basis, 
nevertheless the discussion of the cyclotron does give an opportunity to 
explain the concepts of modern physics more completely than was under 
taken earlier in the book and is very interesting reading. 


’Published by D. Van Nostrand Company, Inc., New York, N. Y., 1944. 




















enti 
ed 
ninent 
lika 
» aSa 
ler the 
is and 
the i 


iterest 
y little 
uch of 

The 
‘tron 
niques 
photo 


imiliar 
asting, 
gh the 
| some 
roduc 

gative 
slectri- 


‘Ss, and 
pment 
nodern 
ing the 
licated 
lectrou 
elabo 
field of 
author 
ywerfu! 
suppl) 
‘ial ap 


e basis, 
nity to 
under 








ABSTRACTS OF TECHNICAL ARTICLES 

Impedance Concept in Wave Guides.’ S. A. ScHELKUNOFF. The im- 
pedance concept is the foundation of engineering transmission theory. If 
wave guides are to be fully utilized as transmission systems or parts thereof, 
their properties must be expressed in terms of appropriately chosen im- 
pedances or else a new transmission theory must be developed. The gradual 
extension of the concept has necessitated a broader point of view without 
vhich an exploitation of its full potentialities would be impossible. 

In the course of various private discussions it has been noted that there 
exists some uneasiness with regard to the applicability of the concept at very 
high frequencies. In part this may be attributed to relative unfamiliarity 
vith the wave guide phenomena and in part to the evolution of the concept 
itself. Some particular aspects of the concept have to be sacrificed in the 
process of generalization and although these aspects may be logically un- 
important, they frequently become psychological obstacles to understanding 
in the early stages of the development. For this reason several sections of 
this paper are devoted to a general discussion of the impedance concept be- 
fore more specific applications are given; then by way of illustration it is 
proved that an infinitely thin perfectly conducting iris between two different 
wave guides behaves as if between the admittances of its faces there existed 
an ideal transformer. This theorem is a generalization of another theorem 
to the effect that when the two wave guides are alike, the iris behaves as a 
shunt reactor. Actual calculation of the admittances and the transformer 
ratio depends on the solution of an appropriate boundary value problem. 

More generally, wave guide discontinuities are representable by 7-net- 
works. In some special cases these networks lack series branches, and in 
other cases the shunt branch. 


Theory of Cathode Sputtering in Low Voltage Gaseous Discharges.® 
CHARLES Harp Townes. To determine the amount of sputtering in a 
glow discharge three functions must be known: the number of ions of a given 
energy striking the cathode, the amount of cathode material released froin 
the cathode by each ion, and the fraction of material released from the 
cathode which diffuses away. Expressions derived for these allow deter- 
mination of the dependence of total rate of sputtering on the geometry of 
the discharge, pressure, cathode fall, current, and constants of the gas and 
cathode surface. The result is most accurate for very low voltage, high 
pressure discharges. Comparison with experimental! data shows quanti- 
tative agreement under these conditions. 


‘Quarterly of A pplied Mathematics, April 1944 
5Phys. Rev., June 1 and 15, 1944. 
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